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1 Related Applications 

5 This application claims priority to provisional application no. 60/168,754 filed on December 
6, 1999, titled, "An E-Commerce Infrastructure for Value Chains", the contents of which 
are herein incorporated by reference. This application also claims priority to provisional 
application no. 60/194,880, titled, "Method and System to Mediate Commerce", filed on 
April 6, 2000, the contents of which are herein incorporated by reference. 



io 2 Field of the Invention 

The invention relates to a method and system for discovery of trades between parties. In 
particular, the invention is a system which allows buyers to define their preferences and sellers 
to define their capabilities, then determines which trading points maximize the utility of the 
buyer. The system suggests trades by exploiting the flexibilities and tradeoffs encoded by both 
15 parties, thus providing win-win trades. A second level of optimization ranks the trades with 
all suppliers, allowing the buyer to rapidly determine the best alternatives. The system allows 
for rich negotiation spaces and supports continuous, discrete, and range or interval decision 
factors. 



3 Background of the Invention 

20 The present invention relates to methods of automatic exploration and exploitation of the 
flexibilities possessed by negotiating parties to uncover improved win-win agreements. The 
invention describes computationally efficient mechanisms that are applicable whether there 
are one or many selling parties. The precise number and types of negotiating dimensions are 
irrelevant as long as they are numerical. Thus the present invention applies equally to the 

25 optimal determination of terms in the purchase of a commodity or an arbitrarily complex 
artifact. 

The past 5-10 years have seen remarkable growth in software tools to help firms with 
enterprise- wide planning (ERP software) and supply chain management (SCM software). 
While these tools do a wonderful job at integrating disparate data sources within and be- 
30 tween firms, the opportunity exists for significant further cost reductions. 
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The same time period has also seen a tremendous rise in the widespread use of the internet 
by both consumers and businesses. Forecasters are predicting that within a few years e- 
commerce between businesses (B2B) and between consumers and businesses (B2C) will grow 
to in excess of a trillion dollars per year in annual revenues. 
5 Electronic markets have proliferated over the last few years with the advent of B2C 

(business-to-consumer) and B2B (business-to-business) electronic commerce. Such market 
places have yielded significant cost savings by lowering the transaction costs between buy- 
ers and sellers. Buyers have also profited through increased competition between suppliers. 
However, electronic markets have hurt suppliers, since the zero-sum negotiation over price has 

10 been at their expense. The present invention describes a tool whereby cost savings for both 
parties are derived from the discovery of win-win trades. Fundamentally, the system works 
by allowing trading parties to describe their desired trade across multiple dimensions and to 
express their flexibility around this ideal trade. Through an algorithmic exploration of their 
flexibilities, the present invention can discover trades that are near the ideal trades of both 

15 parties, enabling both to win. 

The adoption of B2B and B2C electronic commerce was facilitated by the migration of 
catalogues online. This familiar method of presentation ameliorated the significant cultural 
change to electronic trade. For the foreseeable future, electronic commerce will be dominated 
by online catalogs. At present, online catalogues are direct translations of their hardcopy 

20 counterparts where the attributes of a product are described and a price quoted. Inevitably 
however, online catalogs will become more expressive. Catalog entries will be able to represent 
price breaks for large quantity orders, lot sizes, etc. Thus it is important that any software (like 
the present invention) that uncovers mutually beneficial trading scenarios is able to operate 
with such catalogs. Consequently, in the present invention there is an asymmetry between 

25 buyer (usually a human) and seller (usually an online catalog). 

One of the reasons catalogs have come to dominate electronic commerce is that the types 
of goods that can be represented in catalogs are simple. Whether the product is pens or paper 
clips, different vendor's offers differ little from each other (a pen is a pen is a pen), and a quick 
scan of a catalog gives a buyer enough information to make an informed purchase. These types 

30 of goods are low margin and inexpensive. In contrast, the vast amount of purchasing between 
businesses involves materials which are directly connected with business operations - car parts, 
turbines, etc. Such direct goods are the future of electronic commerce. Unlike present-day 
engines, any truly useful procurement tool must be able to support direct materials with 
complex attributes and complex inter- relationships between its components. 

35 Electronic commerce offers unprecedented opportunities for more informed decision-making 

for both buyers and sellers. The past few decades have seen the widespread adoption of enter- 
prise resource planning (ERP) systems, to the point that now almost every major company 
has some form of ERP software. ERP functions as the digital nervous system of a company, 
transmitting and logging information between the company's many different business func- 

40 tions. ERP software keeps track of inventory, monitors the state of purchase orders, signals 
when a company should reorder direct and indirect materials, and a myriad of other functions. 
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Consequently, ERP databases are a rich source of information to optimize a company's oper- 
ations. Yet today this information is rarely used to make more informed buying and selling 
decisions. The present invention can utilize such information sources to optimize a company's 
interactions with suppliers and customers. 
5 One important manner in which this optimization can occur is through an analysis of all 

cost factors. Current buying and selling practices often focus on limited goals, e.g., minimize 
the total purchase price. Myopic purchasing strategies often result in higher total cost of 
ownership when all cost factors relevant to a product in its lifetime of use are included. These 
other cost factors can be significant. Why save the money in taking delivery two days late 

10 if the receiving docks will be full at that time and an additional shift needs to be hired to 
clear the docks? Why order the cheaper drill bit if it is much more expensive to replace when 
it breaks? The present invention improves trades by minimizing the total cost of ownership 
of a product, yielding significant savings to its users. Many total cost factors are difficult to 
quantify - e.g. what is the cost of dealing with a unionized versus a non-unionized supplier? 

15 Consequently, the present invention supports qualitative (best guesses and intuition) as well 
quantitative factors. 

All companies are situated in a supply or value chain. At each step in the chain, a 
company purchases from its suppliers, transforms these inputs, and sells the output to its 
customers. The termination of the supply chain is the sale of the final product to the end 

20 consumer. Since the only influx of external capital comes from the end consumer, companies 
have realized that they compete not only as individuals but also as entire supply chains. 
As result, software products have recently become available which attempt to streamline the 
operations of links within the entire supply chain. This software, variously called supply chain 
optimization (SCO) or advanced planning and optimization (APO), operates on the basis of 

25 forecasted demand at various points within the supply chain. Based on these predictions, plans 
are generated telling companies how much to produce and how to schedule their operations. 
SCO systems are a valuable source of intra-company information - data the present invention 
capitalizes on. Because SCO software relies on forecasted demand, it is only as helpful as 
the forecast is accurate, and, unfortunately, in many cases demand is very difficult to predict. 

30 How can the software know that laundry detergent will go on special at grocery stores in the 
Northeast in 7 weeks? As a result of the volatility in demand and the many other unpredictable 
perturbations that plague supply chains, companies keep significant buffers in the form of 
inventories. In addition to planning, businesses must also be able to adapt to unplanned 
effects. Such adaptation requires flexibility and a means to exploit that flexibility. The present 

35 invention exploits the flexibility of trading parties to streamline the operations of supply chains 
by smoothing the boundaries between trading parties. 

The present invention is therefore a system to allow trading parties to express trading 
desires and constraints across many and varied different factors. These trading preferences 
are informed by many different data sources to optimize for a company's internal operations 

40 and its connections to it's supply chain through an analysis including total cost factors. The 
flexibility expressed by all trading parties is exploited to locate win-win opportunities for all 
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parties if they exist. 

4 Summary of the Invention 

We describe the present invention in its application to facilitating trade between buyers and 
sellers, but note that the mechanisms described are much more general. We can easily imagine, 
5 for example, using the present invention to match individuals (with the desires and skills) to 
projects. 

The inspiration for the present invention comes from utility theory developed by economists 
since the 1960's. Since we are interested in multiple dimensions of negotiation, we draw from 
the multi-attribute utility theory literature. 1 Utility is an abstract concept which has been 
10 formalized in various ways. For the present purposes utility, u, is a number between 0 and 1 
representing a party's willingness to trade. Larger values indicate a greater willingness. 

4.1 The Negotiation Space 

In any negotiation the parties must come to agreement on the factors requiring negotiation. 
We call these factors dimensions or variables. As an example, when purchasing a car, the buyer 

15 may be concerned with price, time of delivery, and color. Each factor price, time, and color 
is a dimension. Most dimensions can be classed as one of three types: continuous, discrete, 
or range/interval. A continuous dimension is one like price for which the buyer's utility varies 
smoothly across that dimension. The buyer's utility at $23 001.00 is almost the same as the 
utility at $23 000. Color is a discrete dimension. Since the car may only be available in 

20 black, white, and silver, the domain of this dimension is the finite set of values {black, white, 
silver}. Moreover, the buyer's utility may be quite different for the three colors. The third 
class of dimensions is called interval dimensions. An interval dimension arises often in B2B 
negotiations. If a machined part is built to some tolerance (e.g., the inner diameter of a screw 
is between 24.5 and 25.5 mm), the range of variability in the dimension is specified as an 

25 interval. In the language of statistical quality control, a certain percentage of the machined 
parts will fall in this range. These three broad classes of variables capture almost all the types 
of attributes relevant to B2B negotiation. 

The present invention operates over any number of continuous, discrete, and range or 
interval variables. We call the negotiation space X and any point in the negotiation space 

30 (x, r) G X. It is important to recognize that the single trading point (x, >r, r) may have 
multiple components, e.g., price = $23 000, time of delivery = 3 weeks, color = black. 

In the present invention, the space of negotiation is agreed upon by all parties involved 
prior to the commencement of any negotiation. We can, however, imagine more dynamic 
situations in which dimensions are introduced and discarded over time. 
1 For a good introduction to multi-attribute utility theory see [1]. 
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4.2 The Buyer's Utility Function 

A party defines it's utility function over this space so that every (x,x-, r) is assigned a util- 
ity number indicating the party's willingness to trade. We indicate the utility function as 
u((x, x, r)). A great deal of work has been done on the appropriate form for utility func- 
5 tions. In the present invention, we take a simple form for the utility function for two reasons. 
First, we would like the form of the utility to be conducive to rapid computation. Sec- 
ond we would like the utility to be simple enough to be easily understood by and elicited 
from users of the invention. With no loss in generality, we write the utility function as 
m((x, x, r)) = exp(— d((x,>f, r))) where rf(x) is interpreted as the distance of trading point 

10 (x, r) from the most preferred trade. 

So that we can operate against seller catalogs, only the buyer needs to define a utility 
function. Across the continuous dimensions, the buyer's utility is defined by specifying the 
most preferred (or ideal) continuous dimensions and the manner in which utility drops off as 
we move away from this ideal. For the discrete dimensions, the utility is specified in tabular 

15 form since there are a finite number of alternatives. Again, the buyer must specify it's ideal 
discrete values and how utility decays away from those values. In section 6.1 we describe how 
this is accomplished. The range dimensions contribute to utility similarly; the buyer specifies 
an ideal range and the utility decays for ranges other than the ideal according to their distance 
from the ideal. 

20 The utility function can also express tradeoffs between variables, e.g., I may take delivery 

in 5 weeks if the price drops to $20 000, or I may accept the white car if I can take delivery in 
2 weeks. The tradeoffs may be between pairs of continuous dimensions (as in the first case), 
between pairs of discrete variables, or between continuous and discrete variables (as in the 
second case). 

25 4.2.1 Normalization and Weighting 

When utility is defined over different types of variables, it is important to normalize the 
contributions of each variable so that the buyer can weight the importance of the various 
contributions to utility. This is a difficult problem. How should a buyer's color preferences 
be normalized so that they can be traded off against time of delivery? The present invention 

30 solves this problem by requiring that the average distance of any negotiation variable from 
its ideal value is the same for all dimensions. Since the buyer is more interested in those 
regions of the negotiation space where the utility is high, the average is weighted by utility. 
This procedure defines a manner in which to define a baseline where all dimensions contribute 
equally. Given this baseline, the buyer can then weight the various contributions and obtain 

35 useful results. 

4.2.2 Utility Elicitation 

Since utility is fundamental to the present invention, its elicitation from the buyer is important. 
Utility may be defined using any of a number of sources: 
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1. graphical user interfaces associated with the invention 

2. standard benchmark criteria applicable to the domain 

3. formal methodologies like the analytical hierarchical process [2], or discrete choice anal- 
ysis [3] 

5 4. inferred through models 

We expand briefly upon method 4. As discussed in the background section, it is important 
to buyers to minimize their total cost of ownership. If we have a function representing these 
costs as a function of the negotiation variables, and perhaps other factors, this function can 
be used to infer a utility function which will act to minimize the total costs. Later we describe 
10 how this can be accomplished. 

4.3 A Supplier's Capabilities 

As noted previously, suppliers are treated differently from buyers so that the tool can operate 
against catalog information with no human intervention required on the part of the seller. In 
fact, we do not require sellers to define a utility at all. 

15 A supplier cannot offer deals at all points within the negotiation space X, e.g., he certainly 

can't offer the black car tomorrow for free! A capability then represents the ability of a supplier 
to deliver and defines a subspace of X. It can include such things as price discounts on large 
volume orders, variation in delivery time as a function of price, etc. Since these relationships 
are already specified by businesses in terms of simple rules like "the price per unit is $10.00 

20 if 1 to 999 units are ordered and $9.50 per unit if 1000 or more units are ordered", suppliers' 
capabilities are represented in the present invention by piecewise linear functions. 

4.4 Negotiation Constraints 

Both parties may have constraints which must be satisfied in order for them to trade. For 
example, the buyer may not buy the car unless he gets it within 6 weeks, or he may not 

25 purchase the car if it is available only in white. These are examples of continuous and dis- 
crete constraints, respectively. A continuous constraint sets a requirement on the continuous 
variables. In the present invention, continuous constraints must be either linear or quadratic. 
Discrete constraints involve discrete variables. A discrete constraint can be expressed as a list 
of the allowed (or disallowed) combinations of the discrete variables for which the trade will 

30 be acceptable. For example, if the buyer would accept either the black or the silver car, the 
constraint would list both these colors as viable. It is important to note that both continuous 
and discrete constraints may involve one or more variables. We can also express constraints 
involving both types of variables by allowing the continuous constraints to differ depending 
on the discrete variables. 
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4.5 Utility Optimization 



With the major components of the invention in place, we describe how the overall system 
works. As a procurement tool for the buyer, there are two levels of optimization. First, for 
any given supplier we maximize the buyer's utility, subject to the supplier's capabilities to 
5 find that trade which makes the buyer as happy as possible. Since we are optimizing within 
a supplier's capabilities, the supplier has expressed a willingness to complete the trade at 
whatever point is determined to be optimal. The tool then optimizes across suppliers to rank 
them according to utility at the optimal point. A graphical user interface allows a buyer to 
investigate the trades suggested by the tool by sorting according to any dimension or by the 
10 overall utility. 

Utility, while a useful concept in assessing an overall score, may be of limited use to a buyer 
due to its abstract meaning. Consequently, we can also apply the total cost of ownership 
function to the results to rank order the suggested trades according to their various cost 
components. Recall that for any trade x € X, the total cost of ownership function returns 
15 the various cost contributions. This additional information aids the buyer in his purchasing 
decision. The utility number for each trade is still useful because the total cost of purchase 
function includes only those cost factors which can be quantified, whereas the utility also 
includes "softer" qualitative factors. 

4.5.1 Aggregation 

20 In addition to optimizing against one supplier at a time, the present invention can also be used 
to optimize against an arbitrary aggregation of suppliers. This is important if, for example, 
no single seller can supply the large volume requested by a buyer. In this mode of operation, 
the buyer specifies sets of suppliers participating in the aggregation and the dimensions over 
which aggregation can occur, and the tool finds the optimal combination in which to distribute 

25 the volume dimension over the allowed suppliers. 

4.6 An e-Commerce Infrastructure for Value Chains 

This patent application also describes an integrated solution for B2C and B2B e-commerce 
that would be built on top of ERP and SCM software and which would provide a number of 
compelling benefits to companies. Amongst the benefits are: 

30 • multidimensional markets which allow consumers to implicitly define their preferences 

over many criteria. This allows both consumers to express what it is they really value, 
allows companies to position themselves clearly in the space of value, and allows for 
efficient matches between trading partners 

• optional anonymity of market participants and their trading desires when that is appro- 
35 priate 2 

2 Anonymity is not, however, a pre-requisite of the proposed invention. 
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• explicit pricing of the flexibility possessed by the consumer and all businesses in the 
supply chain which allows for more robust operation of the entire supply chain. This 
concept is very different from other types of markets (e.g. auctions, reverse auctions, 
exchanges) where transactions are specified exactly. The flexibility introduced by any 

5 party, whether consumer or supplier, is propagated and exploited through the entire 

supply chain. 

• capture and quantification of true consumer demand leading to improved forecasting 
and product development by suppliers 

• automated markets that integrate supply chain networks through coordination across 
10 and within company boundaries. Coupling of the automated markets with local (i.e. 

at the company level) optimization tools fed by real-time company data allows for opti- 
mization and cost savings across the entire supply chain. 

It should be recognized that supply chains may be very different in the near future. Current 
supply chains are based on physical objects made valuable through a sequence of transforma- 

15 tions resulting in a product purchased by an end consumer. With the move to an information 
economy the supply chain of the near future may not involve physical goods at all. In par- 
ticular the entire supply chain may consist of value adding operations converting raw data to 
consumer-desired information. Such supply chains will have the same coordination problems 
current ones do. Our proposed solution applies equally well to these future supply chains and 

20 by supply chain we mean this more general notion. 

5 Brief Description of the Figures 

Figure 1 shows an architecture for the invention. 

Figure 2 shows a schematic of a buyer-specific capability with examples indicating potential 
input. 

25 Figure 3 shows a schematic of a supplier-specific preference with examples indicating potential 
input. 

6 Detailed Description 
6.1 Theory 

In this section we outline the mathematical foundations of the optimization process in sufficient 
30 detail to allow for computer implementation. 

6.1.1 The Negotiation Space 

In Table 1 we define the parameters which collectively define the space of negotiation X. For 
each of the n c continuous variables, we specify an allowed range over which that continuous 
dimension may vary as X{ 6 Xi — \x^,Xi[, where x is the n c - vector of lower continuous bounds 
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n c 


number of continuous dimensions 


n d 


number of discrete dimensions 


n r 


number of range dimensions 


X 


n c - vector of values for continuous dimensions 


X 


rid- vector of values for continuous dimensions 


Xi 


value of ith continuous variable 


Xi 


value of ith discrete variable 



Table 1: Definition of the negotiation search space. 



and x is the n c -vector of upper continuous bounds. Each discrete variable assumes a value 
from within its domain X{ G XV Without loss of generality, we label the domain of discrete 
variable i by T>i = [1, • • • , d{] where di > 0 is an integer giving the number of possible values 
discrete variable x^ may assume. 
5 With these definitions, we define the space of negotiation by the tensor product X = 

X\ ® • - ■ ® X Uc <g) T>i <g> • • ■ <g> T> nd . Range variables are treated separately and not negotiated 
over. 

6.1.2 The Utility Function 

The utility function is a mapping from X into the interval [0,1]. As indicated earlier we 
10 assume the utility to have the form u(k,x) = exp[— e?(x, x)] where d(x,x) is interpreted as a 
distance. In what follows we will assume that in its simplest form the distance function has 
the form 

d(x,x,r) = (x-m(^))*C- 1 (x)(x-/i(^)) +Z(x) + R(r;x). 

Each contribution to the distance function is positive. We consider each contribution to the 

15 distance in turn, beginning with the range variable contribution R(r; x). 

First, we note that the range distance depends on the setting of the discrete variables. 
This allows the buyer to express different preferences for the range variables depending on 
discrete factors. The total range distance is summed up over all possible range variables 
so that R(r;x) = ^27=i Ri( r i'i The vector r indicates the preferred values for all range 

20 variables. If range variable i is specified as the interval r« = (Zii^i) (where r 2 - > r^) then r is 
an n r -vector of such tuples. The distance contribution, from one range variable will depend 
on the application. If the range variables are meant to represent the tolerances on machined 
parts where issues of statistical quality control are important, then the distance between two 
ranges might be related to the overlap between Gaussian distributions. If is interpreted 

25 as a Gaussian having mean (r^ + fi) /2 and standard deviation proportional to r» — then 
an appropriate range distance is given in Appendix A. Other choices for the range distance 
function are certainly possible. 

The continuous distance is quadratic and determined by the positive semidefmite n c x n c 
matrix C _1 . We have allowed this matrix to vary with the setting of the discrete variables and 

30 indicated this explicitly through C _1 (x) . 3 The n c - vector fi may also depend on x and indicates 

3 Where no confusion will arise we eliminate the dependence of C _1 on x for notational simplicity. 



the point at which the utility is maximal - p is thus identified with the ideal value for the 
continuous variables. The precise quadratic form is convenient, but, using recent developments 
in interior point methods, other convex functions are also computationally tractable [4]. 

The discrete distance is determined through the function Z{x) which maps the discrete 
space T>i <g> • ■ • <g> V nd onto the positive real line [0, oo] . In keeping with the assumption that 
distance is a function of only pairs of components Xi,Xj, we assume the discrete distance has 

theform4 

Each contribution Z i:j is a table consisting of didj entries, where ^(xj, x,-) can be interpreted 
as the distance if discrete dimension i has value x; conditioned on discrete dimension j having 
value xj. The diagonal terms Z { offer an unconditional distance. The most preferred value 
for the ith discrete dimension is that for which Zi(^i) = 0. 

Rather than require the user to enter the distances explicitly, there are numerous ways in 
which the distances can be generated automatically based upon a buyer's ranking of preferred 
values. Further details can be found in Appendix B. 
Weighting of Dimensions 

In many cases it is important for simple modifications of the distance function to re-weight 
the contributions to the total distance. If w c is an n c -vector of weights for the continuous 
dimensions, we can accomplish this by letting C" 1 = W C C _1 W C where W c is the diagonal 
matrix W c = diag(w c ). 5 In a similar way we modify the discrete distance to Z W] ij(xi,Xj) = 
Wd-iWd-jZiji^Kj) where w d is the ^-vector of weights for the discrete variables and w d . yi 
is its ith component. The range contribution is also modified so that R w -i{ri) = w r ;iRi{ri) 
where w r is the n r -vector of weights for the range variables and w r -,i is its ith component. For 
convenience the weights are normalized so that (l*w c ) 2 + (l*w d ) 2 + l*w r = 1. With little 
additional complexity the dimension weights can be made dependent on the setting of the 
discrete variables but we will assume throughout that the weights are constant. 
With these modifications, the total distance function becomes 

d(x,*) = ( X -p(x)) t C- 1 {>c)(x-ii{x)) + Z w (x) + R w (r). (1) 

where Z w (x) = £?=i WdA w * t iZii?<i)+Y%Uw w^Z^^, Xj)} and i* w (r) = E£i u>r;iRi(n) 
Assigning weighting factors is useful only if the relevant contributions have been previously 
normalized so that they are all roughly the same magnitude. This serves as the baseline for 
which all weights are equal. The question immediately arises as to what criteria to use to 
weight the distance contributions. 

We shall determine scaling factors, Q d > 0 and Q r > 0, so that the average distances per 
dimension of the discrete, range, and continuous contributions are equal, where by average we 
mean a utility- weighted average over the entire space of possible trades. This weighting places 
4 Later we shall generalize this distance to include weighting of dimensions. 

5 M = [mi,j] = diag(v) is the diagonal matrix formed by setting mi,i = Vi and rrnj = 0 for j ^ i. 
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more emphasis on the better trades 

If d c , da, and d r are the continuous, discrete, and range contributions to the total distance, 
then after multiplication by the scaling factors d = d c + Qddd + Qrd r - The scaling factors are 
determined through the utility weighted average distances defined by 

S> T, r fy <*" Qfd r exp[-Q r <U - Q d d d - d c ] = S„ exp[-Q d d d ]i: r d r exp[-q r rf T ] / v du exp[-dd 

JV du exp[-Q r d r - Qddd - tic] r S* expt-Qd^] exp[-Q r d r ] J v du exp[-d c ] 

S„ Jy du Q d d d exp[-Q r d,. - Q d d d - d c ] = d d exp[-Q d d d ] exp[-Q r d r ] f v du exp[-d e ] 

E„ E, Sv d " exp[-Q.d r - Qrfdi - d c ] Wd expt-Qdd.i] £ r exp[-Q,«i,] JV du exp[-d c ] 

_ E> Sr Jv du d * exp[-Q F d r - Qjdj - d c ] = expE-Qjdj] S% exp[-Q r d r ] f v du d c exp[-d c ] 
~ E v Er JV du «P[-9rir - Qd<*d - d c ] expI-Qidrf] E r exp[-Q r d,.] f v du exp[-d c ] 

A few comments on the above equations are in order. First, J2x indicates the repeated sum 
10 YLx • " " YLn over an possible discrete trades. ]T) r indicates a sum over all the range variables 
and the integral over volume V indicates integration over the continuous trading volume of 
interest. Finally, we have not included a scaling factor Q c on the continuous distance, since 
this can be made equal to 1 if we reinterpret Q r as Q r /Q c and Q d as Qd/Qc- 6 Each of the 
averages is an explicit function Qd and Q r . 
15 The requirement on equal average contributions determines the two unknowns Q r and 

Q d through the equations: (d r )/n r = (d c )/n c and {d d )/n d = (d c )/n c . These two nonlinear 
equations are coupled in terms of Q r and Qd and must be solved simultaneously for Q r and 
Qd- Further details are found in Appendix C. 

6.1.3 Constraint Specification 

20 Buyers and sellers may express constraints over both continuous and discrete variables. 
Continuous Constraints 

For simplicity (and because additional expressiveness is rarely required) we assume that the 
buyer's constraints over the continuous variables are linear. 7 This allows a buyer to express a 
constraint, e.g., the time of delivery must be within 10 days or I will not trade, i.e., t < 10. 
25 We allow for both inequality and equality constraints which can be expressed as G^x = gf^ 
and G 2 6) x < g^ 6) . If there are cf* equality constraints then G^ has rows. Similarly, G 2 6) 
has rows if there are c£* inequality constraints. We allow the constraints to depend on the 
setting of the discrete variables, and to be explicit we often write G^(^), g^ (*c), G^(x), 
and g<V). 

30 Discrete Constraints 

We use a standard methodology to represent and process constraints over discrete variables [5]. 
Abstractly, a constraint over a (perhaps proper) 8 subset of the discrete variables is represented 
as a list of all the allowed combinations the variables may assume. An example representation 
of a pair of discrete constraints is given in Table 2. There are two solutions to this set of 
35 constraints: >c x — 1, ^2 = 2, ^3 = 3 and x\ = 3, *€ 2 = 2, >e 3 = 1. We indicate these solutions 
as {{k u 1), (x 2 ,2)}, (3*3, 3)}} and {(^1, 3), 0*2, 2)}, <X3, 1)}} respectively. Each solution where 

6 We singled the continuous variables out, since there will always be continuous variables in any trading 
scenario. 

7 With little additional complexity we can also handle quadratic constraints, see [4]. 
8 A proper subset of a set is the set itself. 
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2 


1 


3 


2 


3 


1 


3 


1 


2 


3 


2 


1 



>*2 

2 



(a) 



Table 2: An example set of constraints involving 3 variables where the domains of all variables 
are V x = £>2 = £>3 = {!, 2 > 3 I- Constraint (a) requires that the values assumed by >c u x 2 , and 
x 3 are all different from each other, and constraint (b) requires that the value assumed by x 2 
is even. See text for the solution to both constraints. 

all the variables have been identified with specific values from their domains is called a labelling. 

Computationally efficient representations are used to ensure that only feasible combinations 
of variables are ever processed. Numerous third-party libraries offer constraint programming 
functionality. 9 

6.1.4 Utility and Total Cost of Ownership 

The buyer's utility function and associated constraints may be difficult for many users to 
define. In this section we show how models of the buyer's business can be used to define 
utility in a natural manner. 

We imagine a function which provides an estimate of the total cost of ownership for any 
given purchase. Cost contributions to this function might include piece part costs, freight costs, 
setup costs, quality assurance costs, repair costs, etc. It is important to include all quantifiable 
costs associated with the lifetime of use of the purchased product because it is this function 
we will be minimizing. Significant savings may be obtained by taking a longer-term view 
of the purchase. Revenues (negative costs) generated from the purchase are also included 
in the function so that the function represents some measure of profitability associated with 
the purchase. We write the total cost of ownership function as C 0 (x, r;0). We explicitly 
indicate the dependence on the negotiated trade parameters x, x, and r, as well as other 
factors p. The other factors might include forecasted demand, current inventory levels, etc. 
These factors will vary over time, and they can be extracted from the buyer's ERP and supply 
chain management systems (SCM) in real-time just before the purchase to ensure continuous 
real-time optimization. See section 6.2.1 for further details. 

Minimization of C 0 (x,*,r;/9) defines an ideal trade dependent on current conditions: 
XoptOS), *opt(fi), r opt (0). If desired, these can be used to define /i = x opt (0, r = r opt (£) 
and the desired ideal discrete configuration x opt QS) (having distance contribution Z = 0). 
Moreover, the tradeoffs between continuous dimensions around this minimum can be obtained 

9 See for example www.mozart-oz.org or www.ilog.com. 
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through calculation of the Hessian matrix H = [hij] where the i,j matrix element is given by 
= d 2 C Q (x,*r,r;ff) | 

We then identify C _1 with H. In this way, little trading flexibility is obtained in directions 
where total cost of ownership rises rapidly, while significant flexibility is obtained in directions 
where total cost of ownership increases slowly. 

In summary, a total cost of ownership model defines both the most preferred trade param- 
eters and the flexibility possessed around the preferred trade. The model pulls dynamically 
from real-time data sources to provide the most up-to-date optimization based on total costs 
of ownership and other important qualitative factors the buyer may wish to describe in the 
utility function. The same function and its constituent costs may also be used to help analyze 
proposed trades from suppliers. 

6.1.5 Supplier Capabilities 

As discussed in the introduction, suppliers represent their capabilities through a specification 
of the subspace of X in which they will trade. A supplier's capabilities must specify the 
allowed continuous, discrete, and range variables. The allowed range variables are expressed 
as the pairs {r_j,rj), one for each range variable. For example, if a supplier produces 25mm 
inner diameter screws to within a tolerance of 0.5 mm, then the range variable is simply 
(24.5,25.5). These are compared with the buyer's ideal range and contribute to the distance 
function through the H w (>c) function. 

Capabilities over continuous and discrete variables are more complex. 

Continuous Capabilities 

Continuous capabilities are viewed naturally as responses to a buyer's request. Thus we 
distinguish between a buyer's requested continuous vector x^ and a seller's response x.( s \ A 
vector- valued function, f (x< & ) , , x ) returns the response based on the buyer's request and 
also, perhaps, other previously defined supplier responses. Component fa of f defines the ith 
continuous variable, i.e. x^ — /i(x^,x( s )). 

Currently, suppliers are used to quoting price discounts for large volume orders and these 
price discounts are expressed as piecewise linear functions. Consequently, we restrict to 
have the following form (where we distinguish between the functions depending on the buyer 
and seller variables): 

^^f^^ + fi^^ (2) 

k 

An example of how this may be used to define a supplier response is the following: We 
assume three continuous dimensions - price, volume, and time of delivery and indicate these 
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as [xi,X2,X3] = \p,v,i\. Then a response may be formed as 

V W = f VtV (v®) 
*W = =/*,«(««) 
5 P (s) = /p(« (,) ,* W ) = fpA* {8) ) + /p,*(* W )- 

The /„,„ function returns the volume a supplier will fulfill as a function of what the buyer 
asked for. If the supplier can deliver any volume, this will be the identity function. If the 
supplier delivers only in certain lot sizes, this function may have a staircase shape, etc. The 
ft, v function indicates the time it will take a supplier to deliver a certain volume. So, for 

10 example, if larger shipments require longer transportation, then this dependence is given by 
this function. Finally, we turn to the price determination. In this example the price depends 
on the quantity u (s) being shipped and the f PjV might represent price discounts for large volume 
orders. There is also an incremental price contribution based on the time of delivery. If faster 
delivery is more expensive, this is reflected in f p j- 

15 For a given setting of the discrete variables, each fi^{x k s \x^) and f^ k \x[ b \ te^is 

a one-dimensional piecewise linear function. Consequently, the functions can be specified 
by listing the breakpoints. If f-° k \xk\ ^ s) ) has k\'J. breakpoints, then we list these as 
{a4 a) (&ifc = l),x[ s) (b\ s) k = 2),--- ,4 S) ( & ifc = K< fbY° and function values at these break- 
points are = = 2)),--- = *#))}■ s ™y> 

20 ffl, is a piecewise linear function defined by the nfj. breakpoints {xf\b\ b) k = l),x[ b) (b^l = 

2), ■ • ■ , 4 6) ipfl = ^d function values at these breakpoints {/$ (zf (&g = 1)) , /$ (4 &) (&g 

2)), • • • , ffl (x k b \bfl = nfl))}. The breakpoints are indexed by the integers and b^. 

An interval is specified by assigning a value G - 1] and b 2 G I 1 '*?* ~~ ^ SO 

that 11 

25 x%\b#) < 4 S) < 4 S) (&$ + 1) Vi and 4 & >(&$) < xf < 4 &) (&g + D Vi. (3) 

Within each interval the functions are linear, so we have 

k 

where ci(vW, {&$}) = c$ (x^ , ) + In the above equations, the inter- 

cepts and slopes are given explicitly by 



4"cff + i)/g(4"(M) - 4'W)/ff + ')) 



10 The breakpoints are assumed to be in increasing order. 

11 Note that many choices for {b^l} or {bf\} would be inconsistent with these constraints. A simple way to 
fix this is to define a union of all breakpoints (s) and (&), and have all f t ,k have the same breakpoints. 
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,.w.>, _ /^4'^ + i))-/g(4'W)) 

respectively. An analogous result holds for the (&■$■) and (&•$<)• 
5 To eliminate any cyclic dependence on a^ s) we must impose an ordering on xf^ so that 

can only depend on where j < i. Consequently, we can write 

4" = «<*<", {*& • • • {$» + E"^^"' W + E^'K^?- 



Written as a matrix equation, the above becomes 

(I - M«(*« {6}S»)xW = c(xW, {6g}) + M<»({&S2})x< 6 > 

10 where c(x<«>, {&«},{$}) = [c^'M&SmO - {$}, {&$})]*, 

x«=[*« ••■ xW] t ,xW = [arW •■■ *£]', and 



mM(*« {$}) = 



0 o 



15 In most cases x( s ) will depend only on a subset of the variables in x w . If x (s) depends on 
n' < n of the x^ variables, then is an n x n' matrix. In the example given everything 
depending only upon the volume the buyer requested. 

Since (x^) is lower triangular and can be inverted in time 0(n), we can rapidly 
express as 

20 ' ' (4) 

as long as the b\*l are chosen to also satisfy x%\b\*l) < < x^\b\ 8 ^ + 1). These constraints 
will be used in section 6.1.6 which formulates the optimization problem. 

We also allow a supplier to express additional linear constraints so that, for example, 
he may represent that he does not deliver on Sunday. Thus the supplier may define the 
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matrices g[ s) (x), g£ s) (x), and the vectors g[ s) (x), g 2 s) (x) such that g[ s) x^ = g[ s) and 
G 2 s) x( s )x < g 2 s) . g[ s \>c) and G 2 s) (» have c[ s) and rows respectively. 

Discrete Capabilities 

It is easy to imagine that a supplier's response on a discrete dimension is highly constrained 
by the values of the response on other dimensions, e.g., certain product characteristics come 
only in certain colors and package sizes. Consequently, it is not suitable to explicitly define 
a response but only to make available a supplier's constraints amongst the discrete variables. 
Consider 3 discrete dimensions where x\ € V\ — [a,b,c], *r 2 G T> 2 = [A,B,C,D], and x 3 G 
V3 = [a, 13,1,5], and assume the supplier has the following 3 constraints 

Ci (x 1; x 3 ) = {(a, a), (a, /3), (b, /?), (c, /?)}, 
C 3 (^i) = {6,c}. 

We first note that there are 4 feasible solutions (or product configurations the supplier can 
meet): \x x , x 2 , >f 3 ] = [6, 4,0], [c,A,/?L or [c,D,/3]. Feasible solutions to the con- 

straints define the response for the discrete variables. 

We indicate a supplier's or buyer's collective set of discrete constraints by C^(x) and 
C^ b \x) respectively. 

6.1.6 The Optimization Problem 

Having defined the necessary components, we now define the optimization task which deter- 
mines the continuous x* and discrete x* parameters of the trade. 

Since the trade must be acceptable to the supplier, we maximize the buyer's utility over 
a supplier's capabilities. Equivalently, we minimize the distance from the buyer's ideal values 
as 

[x*,**] = arg min {(x( s > - / ^(^)) t C- 1 (^ (s) )(x (s) ~ »{x {s) )) + Qd^(" (s) ) + Q r R w (™)\ 
where 

x^ = (i - MW^D^cixW) + (i - m^{^))- l m^^ 

subject to the constraints over continuous variables 

G x (x< a >)x< s) = gi(" (s) ), G 2 (*rW)x« < g 2 (*r (s) ) 

and the constraints over the discrete variables C^iy^), C^Xv^). In the above, we have 
defined the (c[ s) + cf } ) x n c and (c 2 s) + 4 &) ) x n c matrices Gi(x^) and G 2 (x^) by 

-[$$]. -c^> -[$$>]. 
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The (c^ + cf- 1 )- and (c[ s ^ + )-vectors gi(xW) and g2(^ (s) ) are defined by 

The optimization is accomplished by iterating two distinct phases. Phase one sets the 
continuous parameters optimally for a given setting of the discrete variables. We define the 
5 functions 

di(x,x) = (x-^D'c^WKx-M*)) +R(r;x) and d 2 (x) = Q d Z w (x {s) ), 
The first phase of the optimization is the continuous problem: 12 

x*(x) = arg min di(x (s) ,«) subject to Gi(x)x = gi(x) and G 2 (x)x < g 2 (x). (5) 

A detailed discussion on the solution of the phase 1 optimization problem is found in appendix 
10 D. The second phase determines the best value of the discrete variables by minimizing over a 
function of x alone 

x* = arg min di (x(x) , x) + d 2 (x) subject to C (6) (x) A C (s) (x) . (6) 

Further details on the phase 2 optimization are given in Appendix E. Once x* has been 
determined, we find x* as x* = x(x*). 

15 6.1.7 Aggregation 

Often a buyer may be willing to divide an order between multiple suppliers in order to aggre- 
gate the required demand or to obtain better deals. In this section, we detail how the present 
invention supports this aggregate optimization. 

Aggregation can only occur over the continuous variables where values may be subdivided. 

20 Each continuous variable xi must be parcelled out amongst a set of suppliers. Consequently, 
we extend our notation to xi — >■ Xi t k giving the contribution of the fcth supplier to continuous 
dimension i. The kth supplier may come from a (perhaps proper) subset of all suppliers. 
We indicate the set of potentially contributing suppliers as K, and the number of potentially 
contributing suppliers as |/C|. 13 

25 We restrict the discrete variables to be the same across all potentially aggregated suppliers, 

i.e., we do not generalize X* — > x^jt. This simplifying assumption is made for two reasons. 
First, the size of the discrete optimization problem is smaller and so optimization be performed 
faster. Second, it may be difficult to elicit from the buyer the allowed discrete alternatives for 
each supplier. Nevertheless, this generalization is straightforward should the need arise. 

30 This simplifying assumption requires that the union of discrete supplier constraints Cjc (x) = 

12 No optimization is required over range variables, since these are specified up front by both buyer and seller 
and merely add to the total distance. 

13 All work to this point is thus seen as the special case \fC\ = 1. 
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f\ keK C^\ye) yields a feasible solution when combined with the buyer's discrete constraints 
C^(>c). A necessary (but not sufficient) condition for satisfaction is then that each constraint 
satisfaction problem k having constraints C (6) (>f) A c[ s \x) has a feasible solution. 14 Hence- 
forth, we will assume that the set of suppliers, /C, satisfies this condition. If not, those suppliers 
5 violating the constraints C^(x) A C k \}f) are eliminated from consideration in JC. 

Discrete Search 

We must search over the subsets of JC for feasible solutions, which is a combinatorial problem. 
Fortunately, given a complete labelling of variables, determining the largest subset is easy. 
For any given labelling of all discrete variables, if each A V k G k C fC is satisfiable, 

10 then the union C (6) A c£ s) where = f\ k£K c[ s) is also satisfiable under the same labelling. 
The largest subset of variables is found by adding all k which have feasible solutions with 
the buyer. We needn't worry about smaller subsets because the continuous optimization will 
assign zero values to those if appropriate. Consequently, for any given labelling x we let JC{*c) 
represent the maximal subset of suppliers for which C^{»e) A Ck(x) is satisfiable. It is this 

15 set of suppliers which enter into the continuous optimization. The number of participating 
suppliers is denoted by |/C(*r)|. 
Continuous Optimization 

Optimization over the continuous variables is carried for each labelling x. Generally speaking, 
the buyer's utility will not be an explicit function of x^k but only of X{. We assume a linear 
20 relationship between these two quantities so that 15 



The n c \JC{x)\ vector x is denned as x* = [xi,--- ,x n J where xf = - ,£i ; |jc(ae)|]- The 

n c x n c \JC(x)\ matrix E ^ is assumed known and typically has the form 16 

'£[ 0* ••■ 0*" 
w= 0* & ••• 0* 

0* 0* ■•• 

25 where 0 is the K- vector of all zeros and & is the linear combination relating Xi to the so- 
under our assumptions for 3, X{ = £*Sq. In cases where the buyer wants to accumulate the 
results from suppliers (e.g., aggregating quantities) f = 1 is the |/C(j*)|- vector of all Is. In 
other cases the buyer may take £ = l/\K(je)\ so that the time of delivery becomes the average 

14 It is easily seen this requirement is not sufficient by having C (t) = {(>ci,l), (ki,2)}, C[ s) = {(xri, 1)}, and 
<#> = 2}}. 

15 We can also relate x and x by x = Sx+i/ for some constant vector v, which will not cause any complications. 
However, there seems to be little practical reason to do so. 

16 With no additional computational complexity, we can allow S to depend on x. 
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time of delivery across the suppliers. Constraints on x become constraints on x via 



Gi(*)x = gi(x) and G 2 (x)x < g 2 (x) (7) 

where G x (x) = G^E and Gi(x) = Gi(x)S. We might also expect the buyer to add 
additional linear constraints, such as requiring the latest shipment from any supplier to arrive 
earlier than a certain date, or requiring all deliveries to arrive the same day. There can also 
be constraints specific to particular suppliers, e.g., the buyer doesn't want any more than 100 
units from supplier 5. These can be handled simply as constraints on the individual x i;k and 
added as extra rows to Gi(x), G 2 (x), gi(x), and g 2 (x). 

With aggregation, the quadratic form to be minimized is (Hx-/i(>f)) C~ 1 (>t) (Sx -/*(.*)) 
subject to the constraints given in Eq. (7). This minimization can be carried out through a 
straightforward generalization of the method given in Appendix D. 



6.2 Implementation 

In this section we outline an implementation of the entire e-procurement invention. We begin 
with a high-level description of the architecture, then fill in the details by describing a complete 
object model. 

6.2.1 High-level Architecture of the Invention 

There are at least two modes in which the invention may be used. First, the invention may 
reside at the site of large buyers, and suppliers who wish to sell to the buyer may be required 
to submit their capabilities via a web interface to the buyer. The invention may also be used 
within a marketplace hosted by a third party. Buyers/sellers log onto the market, submit their 
preference/capabilities, and act on the results. The architecture is modular enough to support 
both modes of operation. 

In Figure 1 we present an architecture for the invention. We describe the architecture, 
starting from the optimization algorithm which finds matches between buyers and sellers and 
work our way outwards. 

A controller surrounds the optimization engine, feeding it buyer preferences and seller 
capabilities. If multiple optimization processes are running (perhaps on different machines), 
the controller can also do load balancing, forwarding the request to the least busy process. 
The controller decomposes preferences and capabilities into their constituent buyer- and seller- 
specific versions (see below), selects the most specific matching preference/ capability pairs, and 
sends them to the matching engine for optimization. The controller then collects responses 
from the matching engine and returns them to the buyer. Additionally, the controller logs all 
results into a database for recording purposes. 

Another layer, called the Connector in Figure 1, separates the graphical user interface 
(GUI) through which users communicate with the tool from the controller. This layer serves a 
number of functions. The connector transforms the description of preferences and capabilities 
from the GUI into a form suitable for the implementation of the matching engine. Part of 
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this transformation involves validation of appropriate input from the GUI layer so that no 
malformed input is ever sent to the controller. The Connector layer can also pull data from 
ERP or SCM systems and automatically infer preferences (using the total cost of ownership 
function) for the buyer. The enterprise abstraction layer insulates the Connector from the 
precise details of the manner in which the ERP and SCM data needs to be gathered. Total cost 
of ownership is evaluated in the simulation modules, which may either be running locally at the 
client's site or running centrally at the main server. These simulation modules pull operational 
data (the vector 0)from the enterprise abstraction layer. A preference optimization module 
(TCO) minimizes the total cost of ownership to determine the ideal trade and the flexibilities 
around the ideal trade. 

At the outmost level, a layer provides integration with the GUI and/or host system. A 
number of administrative systems are expected at this layer. Market administration services 
allow easy definition of trading spaces, the dimensions of negotiation, limits on continuous 
variables, allowed settings of the discrete variables, etc. User administration services allow an 
administrator to define buyers, passwords, spending limits, etc. Supplier services accomplish 
similar tasks on the supply side. Managers for preferences, capabilities, and match results 
ensure that these objects are properly stored in a database. This layer layer also dynamically 
generates the html necessary for presentation of the data via a web interface to buyers and 



For maximal portability, communications between the View and Connector are via XML 
documents. For maximal efficiency, communications between the Connector and matching 
controller are as serialized Java objects. 

6.2.2 An Object Model for the Invention 

The fundamental objects required for the invention are preferences from buyers, capabilities 
from sellers, and match results returned to all parties. The components of such objects have 
already been considered from a mathematical point of view, and we now describe one possible 
computer representation. 

In this section we describe a complete grammar for the object model. The following 
syntactic conventions are used: 

• (nt) denotes a non-terminal symbol nt 

• [obj] denotes an optional grammar segment obj 

• {obj} denotes 1, or many times the grammar segment obj 

• — >• denotes a production rule for non-terminal symbol. If there are multiple rules, say 
(a), (b), and (c), then these are denoted as 

{nt} -+ <«)«&) | <c). 
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In contrast, a production rule of the form 

(nt) -+{a),{b),(c) 

indicates that the non-terminal (nt) is composed of three grammar segments, (a), (b), 
and (c) 

5 • terminal keywords are in serif font 

Obvious non-terminal grammar elements like (string) and (integer) are not described. 
Supply Side 

To represent capabilities that apply to a specific buyer (perhaps for contractual reasons), we 
have defined a capability to be a list of (buyerSpecificCapability). With one exception, a 

10 buyer-specific capability applies only to one buyer - that buyer associated in the id field of 
the (buyerSpecificCapability). The exception occurs if the id field is * or wildcard. This 
indicates that the capability applies to all buyers. Using buyer-specific capabilities, sup- 
pliers can represent specific capabilities to certain buyers and generic capabilities applying 
to all other buyers. By not including a wildcard (buyerSpecificCapability) and only listing 

15 (buyerSpecificCapability) s applicable to specific buyers, sellers can also represent the fact 
that they will trade only with a subset of all buyers. In cases where both the wildcard 
(buyerSpecificCapability) and a (buyerSpecificCapability) applicable to a specific buyer ap- 
ply, the most specific (buyerSpecificCapability) is selected. 

A schematic of a (sellerSpecificPreference) is given in Figure 2. 

20 We begin at the top level of a capability: 

capability -> {(buyerSpecificCapability)} 

where 

(buyerSpecificCapability) -4 id: (id), 

discrete: {(discreteVarDescription)}, 

25 continuous: { (continuous VarDescription)}, 

range: {(range VarDescription)}, 
[discreteConstraint: (discreteConstraint)], 
i nsta n ce: { (discreteCapability Instance) } 

<jq [aggregation Participation: 0|1]. 

(id) identifies a buyer or group of buyers. Individual buyers are represented by some unique 
identifier (say an integer) and the group of all buyers is identified by the wildcard character 
'*'. So we have 

(id) — ► (integer) | * . 
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aggregationParticipation is a Boolean flag giving the supplier's willingness to participate in 
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aggregate orders to the identified buyer. 

Each of the variable constituent components is described by 



(discreteVarDescription) -4 name: (integer), 

allowedValues: {(integer)} 
(continuous Var Description) — > name: (integer), 

min: (double), 

max: (double) 

5 (range VarDescription) — >• name: (integer).. 

In its simplest form, a (discreteConstraint) is a list of more primitive constraints 

(discreteConstraint) —¥ {(primitiveDiscreteConstraint)} 

where each primitive constraint is composed as follows: 

10 (primitiveDiscreteConstraint) — >■ name: (string) 

variables: {(discreteVarName)}, 

includes: 0 1 1, 

values: (integerMatrix) 

15 (discreteVarName) is the name of the discrete variable involved in the constraint 

(discreteVarName) — > (integer). 

The includes field is a bit. If the bit is 1, then the combinations listed in the values field are 
the allowed values the variables may take on. If the bit is 0, then the combinations listed in 
values are the excluded combinations, i.e., everything in the powerset of the variables is allowed 
20 except those combinations listed in values. The order of the variable names is significant, since 
they will be assumed to be in the same order in values. If there are a variables involved in the 
constraint, and c constraints, then (integerMatrix) is an a x c matrix of integers: 

(integerMatrix) -» (integer Vector) , ■ ■ • , (integer Vector) 
25 (integer Vector) — t (integer), • • • , (integer) 

The (discreteCapabilitylnstance) component is described by 

(discreteCapabilitylnstance) — } mask: (discrete Var Mask), 

[ra ngeCa pa bi I ity : { (range Varlnstance) }] , 
continuousCapability: (continuousCapability) 
30 continuousConstraints: (continuousConstraints) 
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A (range Varlnstance) defines a range variable and has the form 

{range Varlnstance) — >• name: (integer), 
min: (double), 
max: (double). 

The { discrete VarMask) relates to the discussion of 6.2.2. As in Table 3 we have 

( discrete Var Mask) -» {(extendedVar Value)} 

where an (extendedVarValue) is either an integer from the domain of the discrete variable or 
the wildcard character '*': 

(extendedVarValue) —5- (integer) | * . 

(continuousConstraints) describes the hard linear constraints for the continuous variables. 
Since these constraints may be either inequality or equality, we have 

(continuousConstraints) — > [equality: (linearConstraints)] , 
[inequality: (linearConstraints)] 

Both the equality and inequality constraints are expressed through a matrix which is c x n c 
where c is the number of constraints, and a vector which is c x 1. Consequently we have 

(linearConstraints) -> matrix: (doubleMatrix) , 
vector: (doubleVector) 

A (doubleMatrix) is defined by 

(doubleMatrix) -> (doubleVector), • • • , (doubleVector) 
and a (doubleVector) is just what the name suggests - a vector of doubles: 
(doubleVector) — » (double) , • • • , (double) . 
The only remaining undescribed element above is (continuousCapabihty) whose description 

is 

(continuousCapability) — > breakpoints: (doubleListMatrix), 

valAtBreakPoints: (doubleListMatrix) 
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(doubleListMatrix) describes a n c x n c matrix whose elements are lists of (double): 



(doubleListMatrix) — > ( doubleList Vector },•• • , (do ubleList Vector) 
(doubleList Vector) — >• (doubleList), ■ • ■ , (doubleList) 
(doubleList) -^{(double)} 

It is assumed that the rows and columns of the matrix are in some canonical order so that 
we know which continuous variable is referenced. A natural order is the one defined in 
{ (continuous VarDescription) } 

Preferences 

Just as capabilities may be buyer-specific so too may preferences be seller-specific. The same 
rules determining which seller-specific preference to apply are followed. A schematic of a 
(sellerSpecificPreference) is given in Figure 3. 
We define a preference as follows 

(preference) — »• { (sellerSpecificPreference) }[, (aggregatedPreference)] 

i.e., a preference is a list of (sellerSpecificPreference) with an optional aggregated preference. 
We first describe (sellerSpecificPreference) and then consider (aggregatedPreference). 
The (sellerSpecificPreference) is composed as follows 

(sellerSpecificPreference) — >id: (id), 

d isc rete : { ( discrete VarD escript ion) } , 
contin uous: { ( continuous VarDescript ion) }, 
range: {(range VarDescription)}, 
dimensionWeights: (dimension Weights) , 
discreteTradeoff: (tradeoffTables) 
[discreteConstraint: (discreteConstraint)], 
i nstance: { (discretePreferencelnstance)} 

Of these elements, only (dimensionWeights), (tradeoffTables), and (discretePreferencelnstance) 
have yet to be defined. (dimensionWeights) gives the weights of all dimensions that indicate 
their importance. For convenience we break up the weights according to the three types of 
variables. Thus we have 

(dimensionWeights) — > range: (double Vector) , 
discrete: (double Vector) , 
continuous: (double Vector) 

A (double Vector) has been described previously. Each of the corresponding vectors is as long 
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as the number of range, discrete, or continuous dimensions. (tradeoffTables) is an ndiscrete x 
^discrete matrix of (tradeoflTable): 



{tradeoffTables) (tradeoffTableMatrix) 
(tradeoffTableMatrix) ->• (tradeoffTable Vector) ,•• • , (tradeoffTable Vector) 
(tradeoffTable Vector) ->■ (tradeoffTable), ■• • , (tradeoffTable) 

A (tradeoffTable) is simply a matrix of double values. 

Finally, we turn to the last undefined component of a (preference) . A (discretePreferencelnstance) 
3 follows: 

(discretePreferencelnstance) — >• mask: (mask), 

[ra ngeldea I : { (range Varlnstance)] , 
continuousldeal: (double Vector) , 
tradeoff Matrix: (doubleMatrix) , 
[continuousConstraints: (continuousConstraints)] 

The rangeldeal and continuousldeal fields give the desired range and continuous trade pa- 
rameters. The tradeoffMatrix field gives the positive definite matrix expressing the tradeoffs 
amongst the continuous variables. (continuousConstraints) have been described previously in 
the sell-side specification. 

To complete the specification of preferences, we conclude with the definition of (aggregatedPreference) 
Refer to the discussion of section 6.1.7 for details. 

(aggregatedPreference) participants: { (aggSpecification) } , 

contributionType: (contributionTypeVector) , 
additionalConstraints: (continuous Constraint s ) , 
d iscrete : { (discreteVarDescription) } , 
continuous: {(continuous VarDescription)}, 
range: {{range VarDescription)}, 
dimensionWeights: (dimensionWeights) , 
discreteTradeoff: (tradeoffTables) 
[discreteConstraint: (discreteConstraint)], 
instance: { (discretePreferencelnstance)} 

In the above definition, the previously defined elements maintain their meaning. The addi- 
tionalConstraints field allows the buyer to express constraints around the aggregation, such 
as "all orders must arrive on the same day," etc. participants lists the suppliers who can 
participate in the aggregation and their characteristics. Note that if the wildcard supplier 
participates, the order can potentially be aggregated across all suppliers. {aggSpecification) 
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discrete mask 


output 


specificity 


[* * *J 
[xi * *] 

* >tz[ 


{ continuous 1, rangel} 
{continuous2, range2} 
{continuous3, range3} 


0 
1 
2 



Table 3: Example of discrete masks for specifying continuous and range variables which are 
dependent on discrete variables. n x and x 3 signify specific values for the first and third discrete 
variables. The specificity of each mask is indicated in the third column. 

describes information specific to a supplier participating in the aggregation. It is defined by 

(aggSpecification) — > id: (id). 

id identifies the participating supplier and constraints specific to that supplier defined in an 
5 accompanying (sellerSpecificPreference) will be used in the optimization. Additional informa- 
tion may be added as required. The contributionType field is used to define the f vectors used 
in aggregation. The (contributionTypeVector) consists of n c elements indicating the type of 
aggregation for each continuous dimension: 

(contributionTypeVector) -> (contributionType), • • • , (contributionType). 

Possible contribution types include 

(contributionType) -> sum, average, zero. 

sum sets £ = 1, average sets f = l/|£(x)|, and zero sets f = 0. 
15 Masking 

We have allowed constraints, ideal values, tradeoffs, and continuous capabilities to be depen- 
dent on discrete variables. In this section we describe an efficient manner in which to encode 
this dependence. 

The data structure must represent continuous and range variables for all valid discrete 
20 inputs. An efficient way to do this is to use hierarchical definitions. At the top of the hierarchy 
are the definitions of the continuous and range variables for the discrete values = [*,•• ■ , *]■ 
These values apply to all v unless more specialized masks are defined. A more specialized 
mask of the continuous and range variables is specified by defining values for some of the 
components >q. The more components that are defined, the more specialized the definition. 
25 The most specific mask is always used. An example definition for three discrete variables is 
given in Table 3. The response to the lookup \>c\ ^ 2 £3] is {continuous3, range3} if and only if 
X! = xiAx 3 = >r 3 , {continuous2, range2} if and only if xi = %iAx 3 ^ *r 3 , and {continuousl, 
rangel} otherwise. 
Match Results 

30 Returned to the buyer is a list of matches with different suppliers, which can be ranked and 
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viewed in many different ways in the GUI. A (matchResultList) is a list of matchResult: 

(matchResultList) — >■ {(matchResult)}. 

A match result may either be a (singleSupplierMatchResult) or an (aggregatedMatchResult) : 

(matchResult) -> (singleSupplierMatchResult) | (aggregatedMatchResult). 

A (singleSupplierMatchResult) represents the best trade with a single supplier and is composed 
of the following elements: 

(singleSupplierMatchResult) supplierld: (integer), 
utility: (double), 
feasible: 0|1, 
[costFactors: {(double)], 
continuous: {(double)}, 
discrete: {(discreteVarDescription)}, 
range: (range Varlnstance) . 

The supplierld indicates the supplier sourcing this trade and the utility field indicates the utility 
of the trade (which can be used to rank the trades), feasible is a bit indicating whether or not 
a feasible trade with this supplier was found. The continuous, discrete, and range fields list the 
respective trade parameters determined by the matching algorithm. The optional cost factors 
field lists the constituent costs contributing to the total cost of ownership C 0 evaluated at the 
trade point returned in the (singleSupplierMatchResult). 

An (aggregatedMatchResult) returns the optimal trade when the buyer has requested 
aggregation. It is composed of the following elements: 

(aggregatedMatchResult) — )• utility: (double), 
feasible: 0|1, 
[costFactors: {(double)], 

supplierTradeParameters: {(supplierTradeParameters)}. 

As before, the utility field gives the utility of the aggregate trade, and the feasibility flag indicates 
whether or not a feasible aggregate trade was found (there may not be if the constraints 
were too stringent). costFactors can also be returned in C Q is sufficiently general to handle 
aggregated trades. Finally, (supplierTradeParameters) lists the trade parameters for each 
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supplier involved in the aggregation. It is defined as follows: 

(supplierTradeParameters) -> supplierld (integer), 

continuous: {(double)}, 
discrete: { (discrete VarDescription)}, 
5 range: (range Varlnstance). 

6.3 Summary 

We have described an efficient computational procedure in which to encode buyer's trading 
preferences and hard constraints, supplier's delivery capabilities and constraints, and optimize 

10 to find those matches between one buyer and one or many sellers that maximize the buyer's 
utility. By optimizing against both qualitative and quantitative factors, and exploiting the 
trading flexibilities possessed by both parties, the system determines better trades. The tool is 
particularly useful as companies move their direct material purchasing online. By optimizing 
across flexibilities, win-win trades are discovered for both trading parties. 

15 The representation of trading preferences is designed to be expressive yet easily elicitable 

from a buyer, and computationally tractable. The representation of supplier capabilities was 
chosen to parallel the manner in which suppliers already think of their delivery capabilities and 
seamlessly includes volume discounts and incremental costs. These supplier capabilities may 
be part of an online catalog. The representation of the negotiation space is rich, supporting 

20 three types of variables. 

We have outlined a manner in which preferences may be inferred automatically through 
models of the purchasing company. Such models incorporate many cost factors, taking the 
total cost of ownership into account. The system provides trades which minimize the total 
cost and represent significant new savings. 

25 The invention can operate both at a buyer's site, where suppliers input their capabilities 

through an HTML interface to the world wide web or as an embedded part of an electronic 
market hosted by a particular web site. The invention may operate at regularly scheduled 
intervals or sporadically in lieu of current request for quotations (RFQ). The buyer may 
broadcast a RFQ event to suppliers, indicating a time within which suppliers must respond. 

30 At the close of the event, the buyer can use the present invention to assist in the analysis of 
the supplier responses. 

Complex algorithms have been specified which should permit most matching optimization 
to occur in near real-time. The rapidity of optimization, combined with graphical what-if 
tools, allows for analysis and exploration of trades, which should significantly improve the 

35 quality of purchasing decisions. 
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6.4 An e-Commerce Infrastructure for Value Chains 

In this section we describe in detail how the proposed infrastructure delivers on the promises 
made in the Summary of the Invention. We begin by describing major innovations in the 
present invention and how they are all used synergistically. 

6.4.1 Major Innovations 

The most broad invention combines at least four advances: 

1. multidimensional automated markets (hereafter simply markets) which capture many 
aspects of value. 

2. algorithms and interfaces which implicitly allow for consumers to express the preferences 
over the multiple dimensions of value 

3. linked markets allowing for complex assembly of products 

4. specification of the constraints (both logical and numerical) inherent between markets 
to allow for coordinated buys and sells between markets 

In addition, inventions described in a patent application titled, "An Adaptive and Reliable 
System and Method for Operations Management", Application Number 09/345,441 filed July 
1, 1999 (the contents of which are herein incorporated by reference)can also be used in con- 
junction with the present invention. These other inventions are: 

5. the use of models and optimization algorithms to optimally determine the best bids to 
submit to the automated market 

6. the use of subset relations is-a, has-a etc. to automatically construct the constraints 
between markets 

6.4.2 Integrated Picture 

The internet and e-commerce are changing the way consumers and businesses trade with each 
other. One interesting trend is the development of centralized economic hubs (e.g. eSteel, 
ChemDex) through which all e-commerce in a particular domain flows. This trend is expected 
to continue and become more prevalent. The invention described here is the tool that drives 
these economic hubs. 

Before describing the components and innovations in detail we describe the overall archi- 
tecture of the integrated system. For concreteness we focus on the invention as applied to the 
personal computer (PC) supply chain. We note that the ideas can be applied equally well to 
supply chains in other industries. 

The infrastructure is composed of three major components all running on different com- 
puters. A central server (or more likely servers) coupled to the economic hub serves as the 
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source of consumer demand. Other sets of servers act as the trading markets themselves re- 
laying information between buyers and sellers. Finally, running at each market participant's 
site is client software coupled to the relevant markets. 

Consumer Demand: Consumer demand is pulled through the supply chain through a set of 
coupled reverse auction 17 markets. An interface iteratively elicits multi-dimensional consumer 
preferences which is then forwarded to a top level market where suppliers offer to fulfill the 
consumer request. The functioning of the consumer interface is described in more detail in 
section 6.4.3. 

Markets: Markets within the supply chain are represented by servers which link trading 
partners. The market broadcasts demand, collects responses, and forwards results back to the 
source of the demand. See section 6.4.4 for more information on the functioning of markets, 
particularly the coupling between markets. 

Client Companies: Each trading participant (whether buying or selling) is represented 
by client software connected to the appropriate markets. The client software both initiates 
requests and responds to requests from other market participants. The client software running 
at the company also maintains a memory of all past transactions and can eliminate those that 
are out of date. A user interface allows companies to define the operation of their interface to 
the markets. See 6.4.5 for further details. 

Supplier Push As has been mentioned the infrastructure is organized around demand 
pulled through the supply chain by the consumer. Interestingly, the same pull infrastruc- 
ture can also be used by suppliers to push demand. In addition to requesting and responding 
to requests, the market interface running locally at company sites could also be used to ad- 
vertise capability (capability not in response to any particular demand) to the market which 
will then forward the advertised capability to any participant connected to the market. 

6.4.3 Climbing in Preference Space: the consumer interface 

The consumer drives the integrated infrastructure by generating demand. In our example, 
a graphical user interface (GUI) at the PC economic hub provides a centralized interface 
through which all computers are purchased. The GUI provides a number of advantages to 
both consumers and computer vendors. A computer is described across multiple dimensions of 
value. Some of these dimensions include: price, memory size, processor speed, hard drive side, 
removable storage, monitor size and resolution, financing, warranty, delivery date, etc. For 
each dimension the user can identify a preference. For example the consumer may wish 128Mb 
of memory, a 600MHz Pentium III processor with a 15Gb hard drive for $1500 delivered within 
3 weeks. Additionally, it may be important for the consumer to express constraints that must 
be satisfied. To keep things as simple as possible we express an optional inequality constraint 
for each dimension, e.g. price must be less than $1600 and the hard drive must be at least 
12Gb. 

17 A reverse auction is like an auction except that a request is broadcast to suppliers who bid on fulfilling the 
request. Unlike traditional auctions the proposed market clears in one shot, i.e. each supplier has only one 
chance to bid. 
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Optional Dimensions: This example brings up an interesting feature of the dimensions of 
value. Imagine for a moment that one of the dimensions was portability with possible values 
being high portability, medium portability, and low portability. In each of these cases we might 
consider lightweight laptops, heavy laptops, and desktops. Any one of these choices might 
5 then invoke additional situational dimensions. For example, if we select high portability then 
battery life (in hours) becomes an important dimension which is only applicable to laptops. 
There are a number of ways to treat these optional dimensions. The simplest possibility is to 
have separate markets for laptops and desktops with appropriate dimensions for each. A better 
solution is to have a GUI which only turns on the optional dimensions when appropriate. If 

10 low portability is selected by the user the battery life dimension remains lightly greyed out 
and inaccessible. If medium or high portability is required then the battery life dimension is 
activated. In the technical details section we show how this can easily be accomplished. 

Dimension Weighting: Optionally, the user may also specify an importance to each di- 
mension. If hard drive size and price are paramount these two dimensions have high value 

15 while all other dimensions have low value. The weighting is used to identify similarity between 
computers; two computers which are quite different only in one dimension which is not highly 
weighted by the consumer are considered quite similar. See the technical details section. 

Threshold Constraints: A final way in which consumers may specify their needs is through 
constraints expressed for each dimension. The constraints are simple to state. The consumer 

20 specifies a threshold and whether or not he wants to be above or below that threshold. For 
example the consumer may specify that the monitor size has to be at least 17 inches and he 
will not pay any more than $1500. Knowledgeable consumers may define all these settings 
directly or the software may infer settings for less experienced consumers based upon simple 
questions (like the major uses of the computer). 

25 Branding: Consumers may also have brand preferences for a computer or any of its com- 

ponents. The proposed GUI allows for users to express these preferences. There are two 
alternative ways in which branding might work and depending on the application either may 
be appropriate. In the first alternative brand name is a hard constraint and no computers are 
shown to the consumer unless the brand name matches the consumers desire. In the second 

30 alternative brand name is a soft constraint and the consumer merely prefers one brand over 
another. 

In either case the GUI allows the user to specify for each dimension whether or not brand 
matters and if so which brand/brands is/are preferred. 18 If brand is a hard constraint this 
constraint is propagated to the top level market and all suppliers of the top level market so 
35 that no responses not respecting the constraint are ever returned to the consumer. If brand is 
a soft constraint it is incorporated into the distance function d{ as described in the technical 
details section. Otherwise identical components which match the desired brand are closer than 
those components which do not match the desired brand. 

Iterative Hillclimbing and Recombination: Once the consumer has specified a desired com- 

ls We might also allow consumers to rank their brand preferences but this extension seems unnecessary in 
light of the iterative nature of product exploration. 
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puter defined as a point (and perhaps weights and constraints) in the high dimensional value 
space, a query is sent to a top-level automated market. Sellers of computers {e.g. Micron, 
Gateway, Dell) attached to the market respond to the consumer request by returning offers (one 
or many) to the market for candidate computers again represented in the high-dimensional 
space of value. Responses to the consumer are filtered to satisfy the specified constraints. 
The GUI allows the consumer to sort responses based on any of the dimensions (see the 
function in the technical details section) or according to a global measure of distance between 
two computers (see the d function in the technical details section). 

The consumer may modify any of the responses (a mutation) and use this as a resubmis- 
sion to the market. This method of exploring the space of possible computers according to 
the consumer's desires and the supplier's availability can broadly be summarized as climbing 
in preference space. The climbing can begin from the most recently visited computer config- 
uration or from any configuration that has been saved in a history list. 19 It is also perfectly 
feasible to allow the consumer to recombine desirable computer configurations. For example 
the user may ask for configurations which are "in between" two previously examined configura- 
tions and combine the advantages of both configurations. We also might allow the consumer to 
submit more than a single configuration to the top level market and receive responses around 
all submissions. 

Benefits of the Approach: This procedure offers a significant benefit over other approaches 
to multidimensional markets. In systems like OptiMark 20 the consumer (in this case a financial 
trader) must a priori specify preferences over all dimensions (price and quantity). The explicit 
specification of preference over even two dimensions is a difficult task and it is unreasonable to 
expect consumers to be able to do this. The present innovation allows consumers to iteratively 
explore preference space thereby implicitly defining preference in a friendly manner. 

This method of interactive consumer exploration is suitable for any complex product and 
we expect the application to be broad. The present invention has many other applications 21 . 

Technical Details We have mentioned that computer configurations can be sorted on all 
dimensions according to how closely they match the consumers preferences. In this section 
we describe how this sorting might work and how it interacts with the weights and brand 
preferences attached to each dimension. 

The value space defines a notion of distance or similarity between products (in this case 
computers) defined as follows: If there are a total of D dimensions 22 of value then a computer 
c is described as a D-vector of the settings for each dimension. Let i label the components 
of c so that for example c x = memory is C2 = processor speed etc. For example, the memory 
dimension might assume any of the values ci = 64Mb, c\ = 128Mb, or c\ = 256Mb and 
similarly for all other dimensions. Let a» be a binary variable (i.e. a; takes the value 0 or 1) 

19 History lists extend both within and between sessions. 

20 U.S. Patent 5,845,266, "Crossing network utilizing satisfaction density profile with price discovery features" 
dated 12/1/98. 

21 e.g. credit cards where the dimensions might be interest rate, yearly fee, grace period, benefits (like frequent 
flyer miles, etc.) 

22 The D dimensions include both required and optional dimensions. 
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indicating which dimensions are actually used. If cj is the dimension corresponding to battery 
life in hours then aj — 1 indicates this dimension is appropriate while aj = 0 indicates the 
dimension is not appropriate. The vector whose components are the CJj is represented by a. 

A consumer preference also includes optional weights and preferred brand. The weight 
(or importance) of dimension i is indicated by All w { are non- negative and normalized so 
that J2?=i a i w i = 1 - The vector of all weights is represented by w. The brand preferences 
for dimension i are expressed as a (possibly empty) set Bi which includes desired brands. 
B is the set of all S;. A complete consumer description is thus represented by the vector 
C = {c, a, w, B} with components Ci = {ci,a{,Wi,Bi}. 

The distance d(c, c') between computers c and c' is then defined as 

Ei=i °i 

where measures the distance along the single dimension i. To allow for standard 

comparison between computer configurations we normalize all di so that they lie between 0 
and 1. 

We recall from previous discussion that for soft brand constraints the matching of a brand 
to its desired value is incorporated into the distance function. Thus we generally write the 
distance function di as the sum of two terms, one for brand and one for similarity. If hi is an 
indicator variable whose value is 1 if the consumer has identified that brand matters for the 
ith dimension {i.e. Bi is non-empty) and 0 otherwise, then we write 

di(Ci,C$ = OibngHBuBi) + (! - oabi)4'\ci,^). (9) 

The oti parameter determines the relative importance between brand name and similarity. 
Lacking better information we might take all ai = 1/2. 23 

For simplicity we assume that the brand distance function, df^ (a, c£), is 0 if there is overlap 
in the preferred brands for Cj and c\ and 1 otherwise, i.e. 



df{B i% B[)^\ 

1 0 otherwise 



(10) 



The similarity distance measure, d^(ci,cQ, is only slightly more complicated. The allowed 
values for most dimensions are discrete over a finite range. If the ith dimension is bounded 
and ranges from c™ m to c?* ax then we might define 

^'fe^J^L (id 

though other choices are certainly possible and may be more appropriate depending on the 
nature of the dimension. In many situations it is appropriate that the distance function is not 

23 Better values for a z might be identified from real consumer interaction with the GUI. 
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symmetric, i.e. ^ a) (ci,cJ) ^ ^i'Vi*^)- Exam P les of asymmetry are easy to find. A consumer 
may be more than willing to accept an extra few Mb of memory but will be very dissatisfied 
receiving a few Mb less. As a concrete example we may offer open ended ranges on most 
choices. The range is defined by a single point x\ and an indication of whether the range is 
above or below that point. In the case of the range [x[, oo) we can take the distance function 
measuring the distance from the range to be d(x i ,x' i ) = (x\ - x i )6(x' i - a?»). In the opposite 
case where the range is [0,a£] the distance function can be d(x i ,x' i ) = {xi - xQOfa - xQ. 
Accordingly, the present invention also includes possibilities for d^(ci,c£). 

Distance Ranking: The distance metric of Eq. (8) allows for a global ranking (after 
configurations not satisfying the constraints have been filtered out) of vendor responses based 
upon their similarity to the original consumer request (smaller distances are more similar). 
If C req describes the consumer specified computer and a labels responses from vendors in 
response to the request then the automated market sorts the responses based upon d(C req , C a ) 
and returns a specified number back to the consumer interface. The interface also allows 
sorting (in order of decreasing distance) on each dimension according to the dj. The consumer 
then iterates the process starting from any of the returned offers, modifying them slightly 
and re-submitting to the market. In this way the consumer can move about in the space of 
preferences climbing towards ever more desirable computer configurations. When the consumer 
finally does identify a configuration that he desires to purchase it is a simple matter of hitting 
a buy button. 

6.4.4 Integrating the Supply Chain: cascading markets 

It is natural to consider the extension of high-dimensional automated markets to other trans- 
actions made in the supply chain. In fact, this may be essential to fulfill the kind of real-time 
consumer preference exploration described in section 6.4.3. These other markets will also 
typically be high dimensional but the dimensions will differ from the top-level market. For 
example there may be a hard drive market where buyers of hard drives (whether computer 
vendors, consumers, or value added retailers) meet sellers (e.g. suppliers like Seagate, Western 
Digital). The dimensions of the hard drive market might include: price, size (in Gb), type 
(IDE vs SCSI), quantity, delivery date, etc. 

Operation of Cascading Markets: The basic idea is straightforward. Before a computer 
vendor responds to consumer requests in the top-level market the vendor may send its own 
requests to submarkets to purchase products and services it needs to meet the consumer 
demand. Suppliers to the computer vendors may in turn go to other submarkets before 
responding to computer vendor requests. A single consumer request may then result in real- 
time cascading inquiries to a tree of nested markets all the way down to the beginning of 
the supply chain. We call the market which originated a request the parent market and the 
submarket a child market. It is a requirement that any possible path through the cascading 
markets is acyclic. 

When the consumer moves on to examine a new configuration all relevant market par- 
ticipants (i.e. those that responded to the initial customer request or any concommittant 



subsequent request) are informed and their offers are released. If, on the other hand, the 
consumer elects to purchase the configuration all participants are informed and are held to 
their offers. When combined with emerging XML standards for e-commerce (e.g. RosettaNet 
in the IT industry) the necessary "paperwork" for confirmed purchases can be generated by 
the computer and sent between all relevant parties (including the consumer after a credit card 
check) . 

Note that it remains unlikely that every consumer purchase will trigger a cascade of trans- 
actions throughout the tree of markets because it is not profitable in most supply chains to 
order and transport in single unit batches. Thus, only when inventories need to be replenished 
will companies in the supply chain go to the automated market in response to a request. 

Auxiliary Information: It is useful to pass additional information between markets beyond 
merely the parameters describing the trade. There are many reasons extra information may 
be desired and we describe two applications of auxiliary information. Most simply, in cases 
where it is important to identify who the buyer or seller actually are an additional identifier 
may be passed labeling each trader. 

A more serious application addresses a potential problem of nested markets. A problem 
with cascading markets is that demand can spuriously be amplified as it moves between mar- 
kets and away from its originating source. Consider a consumer demand for 100 computers. 
If two of the computer suppliers are both low on hard drives they might both submit a re- 
quest to the hard drive market for an additional 100 drives. If the hard drive company were 
to estimate demand based solely upon the hard drive market then in this case the company 
would falsely assume that the demand was twice as high as it really is. This problem is easily 
fixed by supplying auxiliary information along with trade parameters. A request that comes 
into a toplevel market is assigned an identifier indicating the market and a unique identifier. 24 
Every request to submarkets that is in response to a request in the parent market is again 
assigned both market and trade identifiers. These identifiers are appended to the identifiers 
from the parent market when the request is broadcast to all suppliers listening to the market. 
Returning to the hard drive example, if the top-level market has market identifier ASSEMBLE 
and the hard drive market has identifier HARDDRIVE and if the original consumer request had 
trade identifier 4561133 and the market trade identifiers were 56891 and 56892 respectively then 
the identifiers which are passed to all hard drive suppliers in response to both hard drive de- 
mands are ASSEMBLE:4561133,HARDDRIVE:56891 and ASSEMBLE:4561133,HARDDRIVE:56892 
respectively. In this way all hard drive suppliers can see that the request originated with the 
same demand. In fact such mechanisms might also allow for improved forecasting of demand 
since all demand is accompanied by its path through the entire supply chain which may contain 
additional useful information. 

Expiration Dates: One of the advantages of the present e-commerce infrastructure is the 
real-time setting of trade parameters. Companies can use the system to optimize their op- 
erations continuously. To achieve the maximal benefit from continuous optimization it is 
important that the offers made by any market participant have a bounded lifetime. What 

24 All trades in that market will have distinct identifiers. 
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is optimal now may not be in a day (or even an hour). Thus an additional application of 
auxiliary information is the inclusion of expiration dates after which the trader parameters 
are invalid. Unless perishable trades are executed before the expiration date they are removed 
from the system. 

5 Other B2B Market Issues The cascading submarkets are all business to business markets 
and as such can operate under different rules than the business to consumer market. In this 
section we address what some of these differences might be. 

Clearing Mechanisms: We have proposed that the B2B markets clear in a fashion similar 
to the B2C market. One important difference is that these markets probably must clear in 

10 one shot in order to deliver real-time feedback to the end consumer. There are certainly 
other ways in which the B2B markets clear. One possibility is to use optimization tools to 
define satisfaction surfaces just as in OptiMark. In this way the parameters of the trade are 
determined by software to optimize some joint criteria. Ideally, we would like protection for 
other B2B market clearing mechanisms. 

15 Aggregation: Another potentially important function of automated markets is the aggre- 

gation of orders to meet demand. For large orders it may be the case that no single supplier 
can meet the requested demand. Market software can automatically collect responses to fulfill 
orders. This software, in addition to coordinating between purchases in different markets (see 
section 6.4.5) can also be used to coordinate multiple orders within the same market. 

20 Many-to-many markets: Thus far the market mechanisms we have described all result in 

one-to-one trades (except if we aggregate orders in which case a trade may be many-to-one). 
For future reference we also note that the markets themselves may be more complex and used 
for many-to-many trades. 

6.4.5 Coordination Between Markets: inter-market language 

25 The assembly of a product or service from constituent products or services often requires 
coordinated purchases from a number of supply markets. 

Logical Constraints: As illustration, there is no point in ordering a part if transportation 
can not also be arranged at the appropriate time. Both events must occur or else neither of 
the purchases should be made. This is an example of a Boolean logical constraint (in this case 

30 the logical and function) that exists between multiple markets. Of course in more complicated 
situations there may be more complex logical constraints between markets. 

For added expressivity we propose using linear logic as the language in which to express 
logical constraints. Linear logic offers compelling advantages over Boolean logic since it sub- 
sumes Boolean logic and is resource sensitive exactly as real supply chain transactions are. 

35 Efficient implementations of linear logic exist so that checking of linear logical constraints is 
straightforward. The present invention includes support for both Boolean and linear logical 
constraints. 

Numerical Constraints: There is another important class of relationships that may exist 
between markets. In most complex processes certain events have to occur sequentially; the 
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process can not progress until some condition is fulfilled. Our shipping example also provides 
an example of this. There is no point arranging the purchase of transportation if the only 
available pick up date is before the item to be shipped can be completed. This is an example of 
a numerical constraint that exists between the parameters of transactions on different markets 
5 and optionally some description of the state of the company. Other examples might include 
cost constraints which require that the total expenditures on all purchases to be less than 
some threshold. 

A Constraint Language: Any automated market solution to supply chain coordination 
problems will require these constraint mechanisms and generically these constraints will be 

10 either logical (as in the and constraint) or numerical (as in the cost constraint). We propose an 
Inter-Market Language (IML) in which to express these constraints so that the nested market 
solution always respects the user specified constraints. 

IML constraints are specified as follows. Let be a logical variable which is true if a 
transaction is made on the ith market (and false otherwise) and let mj(x;) be a logical variable 

15 which is true if the transaction occurs with multidimensional parameters Xj. 25 Furthermore 
let s be a vector of numerical values (for example describing current inventory levels, lead 
times, etc) and~s be a vector of logical values describing the state of the company. The two 
classes of IML constraints which both must simultaneously be satisfied are then expressed as 



where B(-) is an arbitrary logical function, there are M relevant submarkets, and the functions 
N=(-) =0 and iV<(-) < 0 express the equality and inequality constraints respectively. 

Operational Use of Constraints: Operationally the coordination amongst different market 

25 transactions is achieved as follows. Software representing the interests of a vendor (and there- 
fore probably running at the vendors site) expressed in IML listens to servers representing 
various markets. When a market request that the vendor wishes to respond to is heard the 
software examines the companies internal state {s, s} and determines what items, if any, need 
to be purchased on child submarkets. The software defines the purchase by defining a point 

30 or the desired parameters for each trade in each market. These requests are sent to servers 
representing each child market. The client software running at each company then waits for 
responses from each child market. If there are M total markets queried and n* responses from 
the ith market then there are a total of ni^i n i possible combinations of deals. The software 
then examines each deal filtering out those combinations which fail to satisfy the constraints 

35 of Eqs. (12) and (13). Thus we refer to this aspect of the client software as filtering. Those 
combinations which pass the constraints can then be ranked (either by a human user or by 
additional company specific criteria (see reliability below and in section 6.4.7) added to the 
client software running at the vendor site) and sent back as responses to the original request. 
We have operationally described how constraints and filters may be used to coordinate 

2S m t and mj(x;) are related by mi — V X( mj(x) where by V x; we mean the logical or of rai(xi) over all 
possible values of the parameters x,. 
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B(mi, ■ • • ,m,M,~s) 
N = {xi,--- ,x M ,s) =0, iV<(x l5 --- ,x M ,s) <0 



(12) 
(13) 
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purchases from automated markets. The converse problem is the coordination of sell orders 
to automated markets. This is less likely to be problematic since usually there is only one 
sell order being responded to at any given time and thus there are no sell-side coordination 
problems. Nevertheless, if there are sell-side coordinating activities that must occur the same 
5 mechanism can be used. A separate set of logical and numerical constraints can be defined in 
1ML for sell orders. Any responses to requests can then be checked against these constraints 
to ensure no responses are sent which violate the selling practices expressed through the 
IML constraints. Sell-side filters can be used to enforce preexisting contracts like specially 
arranged pricing for preferred customers. This may be particularly important when responses 

10 to requests are submitted automatically by software coupled to the companies ERP software. 
Typically such software will appropriately construct a response but when no human is present 
to check responses sell-side filters provide a final sanity check. 

Other Applications of Constraints: An important consideration for any company in a sup- 
ply chain is reliability. The company must be assured that its suppliers meet their obligations. 

15 In current supply chains this consideration is important enough that companies often limit 
themselves to a handful of trading partners with whom they maintain long term relationships. 
In order to fully reap the benefits of automated markets where many vendors may fulfill a 
need we need to resolve the issue of reliability. Filters are one way to address this issues. If 
for some historical reason a particular company doesn't like some suppliers the company may 

20 add these to be filtered. Thus the internal state {s, s} may include references to unwanted 
suppliers. Alternatively, reliability may be used to rank bids that pass all filters. 

6.4.6 Practical Concerns 

In this section we show how real world concerns of scalability, real-time performance, reliability, 
and security can be assured within the proposed e-commerce framework. 

25 Scalability The utility of an automated market is directly dependent on the liquidity (vol- 
ume of trades) of the market. The greater the number of participants the greater the value to 
any individual participant. It is thus essential that the system scale gracefully to large num- 
bers of market participants. Because the proposed e-commerce infrastructure is distributed 
scaling to large sizes is not a problem. 

30 The market is driven by consumer demand pulled through the supply chain. Interfaces to 

consumer interaction can run locally on the consumer's computer as Java language applets. 
The central server (or servers) for consumer interaction is then only responsible for relaying 
consumer information (expressed in compact XML documents) into the top level market and 
returning the list of consumer responses to the appropriate user computer. The workload 

35 is shared across multiple machines with more machines involved the greater the number of 
consumers. 

Similarly, the automated markets themselves serve mainly as relays of information and 
have little computational demands placed upon them. The markets send out requests, col- 
lect responses, and inform winning bidders of actual purchases. The bulk of the necessary 
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computation is in the software which determines the parameters of the bids and responses. 
Since this software runs locally on each company's computer the system naturally exploits the 
parallelism inherent with increased vendor participation. 

Real-Time Response To be most useful to consumers it is important that consumer feed- 
back is in real or near real time. This may be problematic on the rare occasions in which a 
query triggers a cascade down through submarkets. Before the top level market can clear (i.e. 
for responses to come back to the consumer) the tier one market must clear and before the 
tier one market can clear the tier two market must clear, etc. Fortunately, it is likely that all 
B2B markets will be monitored by software so that turnarounds on all markets will be fast. 
Market time-outs beyond which responses will not be accepted ensure timely responses from 
all submarkets. 

Nevertheless, in situations where the response time is slow we can offer the consumer an 
immediate but approximate response and warn the consumer that the response may not be 
accurate. This can be accomplished by caching previous responses and estimating the response 
by interpolation of similar previous responses. Quick algorithms for the interpolation are easy 
to devise as long as the cached data can be accessed quickly as in any modern database. 

6.4.7 Reliability 

For any company to rely on the proposed infrastructure to run its supply chain it must 
be assured that the system is always up and running. For the system to be functioning 
for a company its own local server has to be up as well and all markets must be running. 
Each company is responsible for ensuring that its local servers are operating. The markets 
themselves will be served by a network of computers sharing the responsibility for managing 
transactions. The software which manages transactions can ensure that if a connection is 
broken during mid-transaction {e.g. if a computer goes down) it will be correctly reconstituted 
and resent. Such software is currently available from software companies. 

6.4.8 Security 

Whenever real dollars are being exchanged electronically there will always be important issues 
around security. We address each of these issues. 

Commitment Assurance: It is important the every participant is confident that contra par- 
ties honor their commitments whether in terms of payment or delivery. We have described one 
mechanism that allows companies to automatically filter out those companies with whom they 
have had bad past experiences. It might also be expected that word of any party defaulting 
on its commitments would rapidly spread to all other participants making any default very 
costly for the offender. 26 If these mechanisms are insufficient the infrastructure may require 
each participant to pay an up front fee from which offenders compensate their contra partners. 

26 We might post defaults publicly so that other parties rapidly become aware of this knowledge and add the 
offender to their filters. 
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Authentication: The system must ensure that all parties are who they say they are. This 
is an old problem faced by internet participants and has been solved satisfactorily by systems 
like PGP etc. 

Privacy: It is also important that companies private data remains private. If a competitor 
could intercept a companies responses to the market it would have a distinct and unfair 
advantage. All communications will be strongly encrypted to assure that this can not occur. 

Information Visibility: Similarly we may want different companies to have differing access 
to information. 27 For example data about consumer demand may be obtainable by a partici- 
pant but only at a fee. This information may be served based on the authentication of the 
participant. 

6.4.9 Optimization of the Supply Chain 

The final ingredient in the infrastructure has been alluded to but not explicitly discussed. To 
achieve maximal benefit from the flexibility the proposed invention introduces at all levels of 
the supply chain, real-time optimization systems can be used. This innovation is described in 
a patent application, titled, "An Adaptive and Reliable System and Method for Operations 
Management", Application Number 09/345,441 filed July 1, 1999 (the contents of which are 
herein incorporated by reference). 

Coupling to ERP Software: The optimization systems monitor (in real-time) the state of 
the company and extrapolate from models of the company in its supply chain the optimal 
responses to requests that come in on the automated markets that the company supplies to. 
To achieve the best benefit, the servers or filters running at each company are coupled to ERP 
software to determine the best real-time response. The optimal response will be a function 
of the current instantaneous state of the company (as recorded in the ERP software) coupled 
with modeling and planning capabilities. A model of the likely future short term evolution 
of the state of the company may be used to best exploit the flexibility inherent in meeting 
the demand from the automated market. Such models might also take into account future 
expected demand and future responses from suppliers. 

Another important optimization that might yield significant payoffs is in the tradeoff be- 
tween making a response that is very close to meeting a requester's demand versus the savings 
obtained by perhaps being further from the requesters true desires. 

6.4.10 Self- Organizing Supply Chains 

In the procedure outlined above it is the responsibility of each company to define its busi- 
ness rules as expressed though its operating constraints and ranking criteria. This is not a 
unreasonable request and is a formalization of the company's current operating procedures. 
However, the present system does suggest a novel extension for the supply chain of the future. 
We can view the activities along a supply chain resulting in a finished good or service as an 
exercise in theorem proving in linear logic. Each company in the supply chain expresses its 

27 Tliis will certainly be the case if we wish to charge for access to some of this information. 
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capabilities as predicates and functions (to include constraint) in linear logic. The finished 
good or service is a logical expression of a set of predicates. A path or method of assembly 
of the final good can then be represented as a proof of the logical expression representing the 
finished good or service. If the good can not be built then no such proof will be found. The 
precise manner in which the product is assembled is read off the proof net describing all the 
necessary steps. The beauty of this approach is its simplicity and the availability of theorem 
provers for linear logic. Though theorem proving for linear logic can be NP complete proofs 
of supply chains may be much simpler. To summarize: the mechanism of self-organization in 
the supply chain is proving the final good or service. 

What are some of the axioms describing operations in the supply chain? There is shipping, 
manufacturing from components, etc. The state of the company includes inventory levels, lead 
times, etc. 

6.4.11 Learning in the Infrastructure 

There are many ways in which the data naturally recorded by the proposed infrastructure can 
be usefully used. 

Prediction of customer demand: All the markets (both B2B and B2C) are rich sources 
of data for mining. Minimally, within any market we can use the data for both forecasting 
and inference of true customer desires (i.e. what they really want and now just what they are 
buying). The quantitative and nature of this market data should make mining it simpler than 
traditional data sources. 

Inference of consumer preferences: The present invention also improves upon the efficiency 
of the end consumers search by aiding the search with learning tools. Based on the consumers 
interaction with the top-level market software can extrapolate to guess preferences and speed 
the search. 

While the above invention has been described with reference to certain preferred embod- 
iments, the scope of the present invention is not limited to these embodiments. One skilled 
in the art may find variations of these preferred embodiments which, nevertheless, fall within 
the spirit of the present invention, whose scope is defined by the claims set forth below. 
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9 Abstract 



A system is described which allows buyers to define their preferences and sellers to define 
their capabilities, then determines which trading points maximize the utility of the buyer. The 
system suggests trades by exploiting the flexibilities and tradeoffs encoded by both parties, thus 
providing win-win trades. A second level of optimization ranks the trades with all suppliers, 
allowing the buyer to rapidly determine the best alternatives. The system allows for rich 
negotiation spaces and supports continuous, discrete, and range or interval decision factors. 
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A Distance Function for Range Variables 



The distance between buyer range variable n = (r^n) and seller range variable rj = (r_j,rj) 
is related to the overlap defined as 



overlap (ri, r 



where M(x;ri) is a Gaussian in x centered at m = + rf)/2 with standard deviation o-j = 
<y(n - Li) and a is a tunable parameter. The range of integration is bounded by the seller's 
capabilities. Evaluating the integral, we find 



over lap ir^rj) - 



erf 



- <jjfii — ojHj + erf 



- crjfii - <Jjy.j - erf 



a^2(af + a]) 



By adjusting a we can control the hardness of the boundary at the ends of the range (smaller 
a corresponds to sharper dropoff outside the range). The maximum overlap will occur when 
fj,j = ^ and (j j = 0. Substituting those values in the above equation, we get 



maxOverlap = 



erf I 



V2j 



The distance d(ri,rj) is related to the overlap so that it varies from 0 (for identical range 
variables) to oo (for maximally distant range variables). Consequently, we define 



R(n,rj) = - In 



[overlap(rj,7-j)l 
maxOverlap J 



B From Rankings to Numerical Values 

The numerical values in the tables Z(x t ,JCj) can be formed from the buyer's rankings in a 
20 number of ways. Since we are minimizing distance, we want preferred options to contribute 
little Z and progressively less favorable options to have progressively larger Z values. To 
approach this, we first consider the problem of denning a Z' value normalized to lie between 
0 and 1 and then take Z = - ln[l - Z'\ 

One possibility for Z' is a linear scaling. If di is the integer indicating the number of values 
25 in T>i, r\ a \>ij) indicates the buyer-specified rank (smaller rank means more preferable) for 
the option ^ = a e T>i conditional on jej, then we might set Z'fa — a,Xj) = (r\ a \vj) - 
l)/{d{ — 1). Alternatively, we might use a non-linear scaling and set Zfa = cc,Vj) — (1 — 
_ /3 d *- 1 ) for some 0 E [0,1] (exponential) or set Z(yi = a,vj) = ((rj a) (uj) - 
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type 


Z' 


linear 
exponential 
algebraic 


In 

-In 


n Ri ~ 1 
[Ri-r\ a \vj)\ 





Table 4: Possible numerical ranking functions. 



l)/(Ri — l)) n for some n > 0 (algebraic). The corresponding scaling functions for Z are given 
in Table 4. 

For those allowed values of Xi which are left unranked by the buyer, we can assign them 
a common numerical value greater than the lowest ranked value. For example if n(xj) is the 
5 largest ranked value we might assign all remaining Xi values to have Z(xi, xj) = Vi(xj)/(di — 1) 
in the linear case or Z{x i ,x i ) = or Z{x h Xj) = (l-l/(r;(^) + l))/(l- 

1/di) in the non-linear cases. 



C Determination of Scaling Factors 

In this section we supply the details necessary to determine the scaling factors Q r and Q d - In 
10 the equations defining the average distances we have need of the integrals f v du exp[— (u — 
u)*C -1 (u - u)] and f v du (u - u)*C -1 (u - u)exp[-(u - u)*C -1 (u - u)]. In general, these 
integrations are not analytically tractable; however we will assume that the region of integra- 
tion defined by x and x extends at least 2 standard deviations out in all directions so that the 
volume of integration is effectively infinite. In this case we can evaluate the integrals in closed 

15 form: 

f v du exp[-(u - uJ'C-V - u)] « i/^rr 

and 

f v du(u- uYC- l (u - u)exp[-(u - u^CTV - u)] » ^^C^. 

To this point we have not been specific about what is meant by the average over the range 
20 variables. In fact, we will assume that, however we define this average, it is independent of 
what the buyer requested and thus can be removed from inside the sum over the discrete 
variables. In this case we will concentrate on setting the discrete and continuous components 
equal. We have {d c )/n c = {d d )ln d where 

jd d ) = Q d Ev dd exp[-Q d d d ] f v du exp[-rf c ] ^ (dj _ 1 E v exp[-Q rf rf rf ] f v du d c exp[-rf e ] 
n d ~ n d Y!, v e:x P[-Qdd d ] J v du exp[-d c ] n c n c E v exp[-Q d d d ] f v du exp[-d c ] 
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and cancelling denominators Q d is determined by solving 

®±Y d d exp[-Q d d d ] [ du exp[-d c ] = — Vexp[-Q d d d ] f du <f c exp[-d c ]. 

This can be simplified further given our assumptions about the Gaussian integral: 

Qd d d exp[— Qddd] _ 1 exp[— Q d dd] 
~d^ i/detC-Kv) ~ 2^ ^/detC-Hv)' 

5 In most cases there will also be at least one finite root, which is what we need. If there 

are a multiplicity of roots, we select that root having the smallest value to ensure that the 
distribution defined by the utility function is as flat as possible. If there are no real roots, 
then we select the Qd which minimizes the error. 

D Phase 1 Optimization 

10 In this section we detail the solution to the phase 1 optimization problem defined in Eq. (5). 

D.l General Formulation 

We begin by writing down the objective function: 

di(x,x) = (x< a > -ai(x))*C" 1 (xW -/i(x)) + J R(r;x) 
where (see Eq. (15)) 

x« = (i-mW(*w {&>}))- V*^^ 

(15) 

as long as x^{bfl) < < ^(b^l + 1). This must be optimized over x subject to 
the constraints Gi(>f)x( s ) = gi(^) and G 2 (^)x^ s ) = g2(x) for a given it. Since x is 
fixed throughout the optimization R{v;*c) is just a constant which can be evaluated and 
then ignored as far as determination of x is concerned. We simplify our notation so that 
di(u) = (xW-^C-^xW -p)+R(t). Similarly, 

x« = (l-MW({^}))- 1 c({6g},{^>}) + (i-MWaftSJ}))-^*)^^})^) 
= 8({6Si ) } ) {6S}) + S({6g,{6J3})xW 

where we have defined the n c -vector s = (I - M( s )) -1 c and the n c x n' c matrix S = (I - 
25 M^^) _1 M( 6 ). 28 n' c < n c is the number of continuous dimensions upon which the supplier 
response depends. Putting these pieces together we have 

di (xW) = X W*Q({&g> { & S» x(6) + q4x(6) + const (16) 
28 In the next section we consider S further, taking advantage of any redundant variables to simplify the QP 
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where Q is an n' c x n' c matrix given by 

Q({&&{&$}) = ^({^i,^})^^^,^}) 

q is an n'-vector given by 

<1({6J3, = 2S*({ftW, {6$ })C- 1 (s{bg, {&<»}) - 

5 and the constant term is 

const = ( 8 ({6SS, {bg}) - -Ai)^" 1 (s({b$, {bg}) - + J2(r) 

It is easily shown that the matrix Q is positive semidefinite. di(u^) must be optimized 
subject to the constraints 

GiS({6g, {6<2})xW = gl - G l8 ({6« {bg}), (17) 

10 G 2 S({6«, {6jg})x« < g 2 - G 2 s({bg, {bg}), (18) 

x<*>(») <s({6S;2,{6S|*}) + S({6SJ,{6S})x (6) <x«(i) Vi, (19) 

xW(i) <xW <xW(i) Vi (20) 

where xW(*)* = [a* (6$) ■ • • and xto(»)* = [a^&JJ + 1) • ■ • + 1)] are 

15 the vectors of breakpoints for the (k E [1, n c ]) functions. Similarly, for the buyer variables, 
= [^(ftW) ... x n (hg>)]and xW(0* = [^(bg + 1) -•■ s n (&g> + 1)] are the 
vectors of breakpoints for the (A; G [l,ra c ]) functions. 

di is optimized in two steps. First we assume a value for the integer variables and 
{bf^} fixing the constraints in Eqs. (17) and (18). We then minimize the quadratic d\ sub- 
20 ject to these constraints using a standard quadratic programming (QP) algorithm 29 to find 

4p\({^M©)- This result is then used to find x *(^{ 6 lfc}'{ 6 S» = + 

S({6$}, {b?}})y!St{{b\'l},{b { ^}). By searching exhaustively over {bg} and {bg} we deter- 
mine 

^( 3 <) = argmin * (x* {&$}, {bg}), *) 

25 Before turning to the manner in which this optimization is accomplished in section D.3, we 
discuss a few improvements to the input into the QP algorithm. 

D.2 The Quadratic Program 

In this section we discuss improvements in the specification of the quadratic programming 

task to make it better conditioned and more efficient. 

30 It is possible that the matrix Q used as input to the QP is singular or near singular, 

29 Numerous third parties supply quadratic programming packages, some based on interior-point methods 
which are appropriate, in the present context where matrix Q may be singular. 
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depending on the rank of Si. Consequently, it is important to use an interior point QP 
algorithm for which this is not a problem. 

D.2.1 Variable Elimination 

There may be rows or columns of S which are entirely zero. It is important to take advantage 
5 of this structure to eliminate unnecessary computation. 

If row i is zero, this indicates that supplier variable is independent of any buyer variable 
and is simply equal to Sj. If there are any constraints involving only supplier variable i, then 
these can be checked before calling the QP solver to eliminate wasteful computation. 30 If the 
constraint is satisfied, we can eliminate it and then call the QP function. If column j of S 
10 is all zero, this means that there is no dependence of u( s ) on uf^ and so there is no need to 
optimize over that buyer variable. We will, however, still have to check that the constraints on 
uW are satisfied. This observation allows us to simplify the quadratic programming problem 
by eliminating this unnecessary variable. 

We formulate the problem to take advantage of both these observations. Assume there are 
15 m < n c non-zero rows of S and m' < n' c non-zero columns. Relabel the rows and columns so 




where Si is of size m x m' and, with a slight abuse of notation, the 0's indicate matrices of 
the appropriate size. In like fashion we write 



20 




where x^ is an m vector, x^ is an n c - m vector, x^ is an m' vector, and x^ is an n' c - m' 
vector, x^ indicates those components of x( s ) which vary with buyer variables while x^ 
indicates those variables which are constant. x[ b) are those buyer variables which affect at 
least one supplier variable and x^ are those buyer variables which do not affect any supplier 
25 variables. If we also partition sass*= js* s|] then 

x^ s) = si + Sixf and x^ s) = s 2 . 

Finally, we will also have need of the partitioned covariance matrix 

c -i _ c r,i c ri 

[^2,1 ^2,2] 

where (C^l )* = C^. w since C" 1 is symmetric. 

30 Additionally, the QP doesn't like to have a row of zeros in the constraint matrix. 
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Including the dependencies on ^t, the quadratic form we must optimize is 

[ (S1 - + (*f >ra (S2 - »)<] [ WaC y Wi WaC L? Wa j [ S2 _ ^ 2 j 

where we have made the obvious definitions 

W=f Wl °] and „=K 
[ 0 W 2 J [fr 

Expanding the quadratic, the constant term is 

The linear term is 

q = (( Sl - /x 1 )*WiCr i lWiS« 1 + 2(s 2 - ^'WsCjJWxSOx^ 
and the quadratic term is 

Q = (xftysiWiC^WiSiX^ 

timiz 

GijiSixf 5 = gi - Gi ; isi - Gl; 2 S 2 , 
G^iSiX^ < g 2 - G 2; isi - G 2;2 s 2 , 
x ( i S) W ^si + Six^ <xj s) (i) Vt, 

where we have simplified notation by dropping the dependence on breakpoints and defined 
Gi = [G i; i G i;2 ] 

i.e., G i; i is the columns of Gi involving the active buyer variables while G i;2 is formed from 
the inactive buyer variables. Similarly, the constraints which can be checked prior to even 
solving the QP are 

x 2 s) w <s 2 <4 s) w Vi - 



D.2.2 QP Normalization 

There are two types of normalizations we perform for increased numerical stability. If the scales 
of the continuous dimensions vary greatly (i.e., one dimension has an ideal value of 10 -2 and 
another dimension has ideal value 10 6 ), then the matrix C -1 becomes badly conditioned, which 
makes S badly conditioned and thus makes the QP unstable. Similarly, if the breakpoints have 
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radically, different lengths, the QP problem can also become unstable. We consider how to 
remedy each of these problems. 
Dimension normalization 

The quadratic form to be minimized is (x-/i)*C -1 (x-p). We modify this problem by defining 
5 the dimensionless variables fa through x\ = fam in terms of which 

(x - tfC- 1 (x - ix) = (<f> - lYC- 1 ^ - 1) 

where C -1 = DC^D and D = diag(ji). When making this variable transformation we must 
assure that the constraints on x are adjusted so that Gix = gi and G 2 x < g 2 become the 
constraints Q x <$> = gi and Qify < g 2 on cj> where 0i,2 = Gi, 2 D. If there are globalbounds of 
10 the form x < x < x, then these are translated to 0 < <p < <j> where = D _1 u and <f> = D _1 u. 

In addition to these changes, we also need to modify the capabilities appropriately. Recall 
Eq. (2) which reads (after dropping the dependence on the discrete variables) 

^E^f') +/&>(*?>)■ < 21 > 

k 

In terms of the dimensionless variables we have 

,^_ v /g(4*V*) , /iMV*) (22) 

15 * V ^ /i< 

Thus, when we use the dimensionless variables, the function f^k needs to be modified so that 
the y dimension is divided by m and the x dimension is multiplied by fj, k . 

Once the optimization has determined an optimal 0, we return to the original variables 
through x = D^. 

20 Breakpoint range normalization 

For larger problems it is sometimes important to normalize the input variables for the QP 
solver. This is carried out in order to improve the QP search performance in situations 
where there is large variation in the range of the segments 4> k b \bk) < 4$ < 4>t\ b k + 1) and 
the magnitude of their values. This can be accomplished by introducing the normalization 
25 variables 4>[ b) = {<ff - (h)) / (<j>[ b ) (h + 1) - 4>[ b) (h)). Here 4>[ b) will have minimum value 
of zero and maximum of one. We can now express as 

</>^ = A^ b) +p (23) 

where A is the diagonal matrix A = diag(0 ( 1 bl+1) - #\ - 4>^\ ■ ■ ■ , - <& m) ) 

and v t = [4>f-\4^\--- ^\ 
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Expanding Eq.(16) using the definition of <f> in Eq.(23), we obtain 
<*i(w; {&*,*}) =^ (&)t ADQ({fe i)fc })DA^ 6) + q*({6i >fc })DA^ + p'DQ({& i , fe })DA^ (b) + 
ip t DQ({6 ijfc })Dp + q*({6^})Dp. 

Hence the quadratic programming problem to be solved in terms of 0 is 

QP(*< 6 >; {bi, k }) = \± {b)t Q({bi, k m ib) + q t ({^»u (b ), (24) 

where Q = ADQ({6 i)fc })DA and q* = (q < ({6 i>fc }) + p t DQ({& iiJfe }))DA. Similarly, the equality 
and inequality constraints of Eqs. (17) and (18) expressed in terms of 4> are 

6i({&i )fc })£ (6) = &i and G 2 ({&^})0 (6) < g2, (25) 

where 

G^fc}) = G 1 S({& i)fc })DA, 

G 2 ({fe ijfc }) = G 2 S({6 i)fc })DA, 
gi = gi - G lS ({b ijk }) - GiS({6i, fc })Dp, 
g 2 = g2 - G 2 s({bi, k }) - G 1 S({6 i , fc })Dp. 

D.3 Optimization Over The Breakpoint Variables 

We have detailed how to determine x^(x, {6 -J}, g iven a specification of all {&$} and 

{bfl}, but how do we determine the best feasible choice of {bJJ} and Before we answer 

that question, we have to examine the infeasibilities that might arise due to improper choices 
of{^}and{fcg}. 

D.3.1 Elimination of {b itk } Infeasibilities 

Most of the settings of {&$} or { 6 i**} wiU not y ield a feasible solution because Eq. (20) or 
Eq. (19) will be violated. The reason is that, for any k, it is highly likely to get two segments 
4 6) (M < 4 6) ^ x k ] ( b i,k + 1) and 4 b) ( & i,fc) < x f < x fi h 3,k + 1) that do not overlap (i.e., 
the optimizer will not be able to find a xf^ value that satisfies both bounds, and hence no 
feasible solution will be found). The same is also true for supplier response segment 

To eliminate those infeasibilities, we construct a super set of breakpoints {zo b k } such that 

»?> = 0(4" = i). *?'<»£! = *).- = 

i=l 

Hence each piecewise linear function f ijk will be described by a set of breakpoints zu k = 
{»*(!), ■ • • ,a:fc(«fc)} and function values f it k{x k (l)), ■ • • , fi, k {x k {K k )). These function values 



at the new breakpoints are linearly interpolated from the old breakpoint-function pairs. After 
this reduction we have K k - 1 different b k e [l,Kk ~ 1] values, which eliminates the need for 
constraint Eq. (20). 

The continuous supplier variables and constraint Eq. (19) are more problematic. Sup- 
pressing any dependence on *e, we know 

k<i k 

If {bfl} and {bf}, ■■■ , bf]^} are known, then this expression allows us to bound the values 
xf^ may take on, and thus bound the interval in which the feasible 6^ lie. If we define the 



jsmax(a;,y) if s > 0 jsmin(a;,y) if s > 0 

smaxfs, x, y) = < and smm(s, x, y) = < 

[smm(x,y) if s < 0 [ s max(a;, y) if s < 0 

then the minimum and maximum values that xf^ may take on are: 

mm(^) = • {&S2» + E smin (^( 6 2)'4 s) (^)»4 s) (^ + 1)) 



max(xW) = Ci ({6g,- • • ,6^}, {&<g}) + E^Kji^S)*^^^).^^ + *)) 

+ Esmax(^(^),^(6£^(^ + 1)) (27) 

A; 

respectively. The above equations imply that we can fully translate a complete specification of 
buyer variable breakpoint intervals to a corresponding set of feasible supplier variable break- 
point intervals. We will search over the breakpoints using a branch-and-bound algorithm, 
and will need to be able to infer the restriction on the values of feasible supplier breakpoint 
intervals given a partial specification of buyer breakpoint interval variables. This inference, 
which narrows the domains of the supplier breakpoint variables, is called constraint propaga- 
tion. Before describing how the branch-and-bound algorithm search over breakpoint interval 
variables is accomplished, we show how the constraint propagation may be obtained. 
Narrowing of Supplier Breakpoint Domains Through Propagation 
Let {if" 1 be a vector of integers specifying a set of breakpoint intervals for the fcth buyer 
variable. \fi k \ = \ if the interval is uniquely determined, and \fi k \ — K k b) - 1 if no restrictions 
on the interval have been specified. We index the Zth component of 0$ as pfj, I E [1, nf^ - 1]. 
A set of breakpoints for all buyer dimensions is written as {fif*}. In analogous fashion, we 
define P k *\ I G [l,*^ _ 1 ]' and i^k} for the supplier breakpoint interval variables. 
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Prom a specification {/3jj. 6) } we infer the feasible domains {P k s) } using max(a;[ s) ) and 
min(ar^) where: 



and 




We begin with i = 1 (where no supplier variables contribute to x^) and determine the interval 
h = [min(x ( 1 s) ),max(4 s) )]- We then form the vector fi ( ' ] with components {P < £j\[xi\pi;i),Xi ) (Pi;l+ 
1)] n Ji ^ 0}. Equations (28) and (29) can then be iterated to find the interval I 2 = 
[min(4 s) )> max(4 s) )] using the previously calculated min(a; ( 1 s) ) and max(x[ s) ). With J 2 known, 
the components of P^ are determined as above. Iterating this procedure, we can find the 
potentially narrowed domains of all {P^}- 

We now return to the optimization problem: we need to search over all Yl k 4 &) x Yl k 
combinations of breakpoints to determine the combination that minimizes di (x opt {b k ^}, }),*)■ 
The simplest manner in which this accomplished is via exhaustive enumeration over all {b^} 
and {b[ s) }. However, once all b[^ have been specified, we need only search exhaustively over 
the corresponding inferred as above. There are, however, many ways to speed up the 

search over the buyer breakpoint intervals. We describe two methods. 

D.3.2 Local Search 

This method is extremely fast but not guaranteed to return the global optimum across all 
buyer breakpoint intervals. Here, initial random and feasible 31 {b[ s) }, and {b k b) } sets are first 
created. After the continuous optimization is carried out using the initial {b k ^}, and {b[^}, a 
new {b^}, and {b[ b ^} set is generated from the old one by examining whether each generated 
4 6) ({&£°},{4 6) }) solution is on the boundaries of the interval x[ b \b k ) < xf < x {b \b k + 1) 

or 4 S) (M < 4* ] < 4 S) ( 6 * + !)• If 4 6) is on the lower bound of the interval > &i 6) = h f ~ !; 
similarly, if it is on the upper bound, b[ b) = b k b) + 1. If xf* is not on either end, then 
b[ b) = bf 1 . Of course this done while maintaining b k b) in the interval 1 < b[ b) < - 1. 
Supplier breakpoint variables are updated in a similar manner. The procedure is repeated 
until no decrease in distance is attained. 

31 Feasible in the sense that they do not violate the constraints in equations (17), (18), (19), and (20). 
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D.3.3 Branch-and-Bound Search 

Branch-and-bound is an efficient procedure used to exhaustively search a space and eliminate 
configurations that can be proven to be suboptimal. The theory behind branch-and-bound 
is well understood [5]. The idea is to structure the space as a tree of partially-completed 
possibilities and to explore the tree so that good complete solutions are constructed early on. 
Branches of the tree that can be shown to yield solutions which cannot be better than the 
currently best solution need not be explored. Heuristics indicate the order in which to explore 
the tree and provide a lower bound on the fitness of partially completed solutions for pruning 
of suboptimal branches.. 

In the present context, we carry out branch-and-bound search over the space of possible 
buyer breakpoint combinations. The tree is structured beginning from a node at which none 
of the have been assigned. Prom this node, the tree branches into the nf } - 1 possible 
assignations for bf ] . From any of these nodes, the tree branches into the nf possible values 
for bf\ The tree branches at each node to assign the remaining variables. At the leaves of 
the tree, all buyer breakpoint variables have been assigned. 

From any node, the search proceeds by choosing the most promising node to consider, 
as indicated by a variable ordering heuristic and considering the possible values for that 
variable in the order indicated by the value ordering heuristic. A complete solution and 
its accompanying distance is obtained when the search ends at a leaf of the tree. To avoid 
unnecessary exploration, a lower bound heuristic indicates when a partially completed solution 
(node) need not be expanded because all completions will have higher distance than the best 
located so far. We describe these three heuristics for search over buyer breakpoints. 

Lower Bound 

Given a set of buyer breakpoint interval vectors {fif^}, we can propagate to {^ s) } and the 
intervals [min(a;[ s) ), max(x[ s) )]. A lower bound on the value d x can be obtained by ignoring 
the supplier capabilities and minimizing 



subject to 



d 1 (x) = (xW-/x) t C- 1 (x( s) -/i) 

Gix<"> = gi 
G 2 x (s) < g 2 
min(a£°) < x\ s) < max(x\ s) ), 



(30) 



(31) 
(32) 
(33) 



The lower bound thus obtained is denoted by lb({^^})- 
Variable and Value Ordering 

The lower bound lb({fi^}) can be used to define an ordering over variables and the values 
any variable may take on. We define the functions L k j which operate on the set {fi®} by 
locking the fcth buyer breakpoint variable to 0 k;l . At any node we can determine the lower 
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bounds lb(L k>l {p( bS >}). We can select to explore the variables in order of increasing lb* where 



i.e., the average over all possible completions for b[ b \ 32 Within any given variable the (say 
k*) values can be explored in order of increasing lb(L k * . If we are willing to do 

more computation, we can explore the branches in an interleaved order according to all the 
lb(L fc ,K{^ (6) }))- 
Pruning and Search 

During the branch-and-bound search, branches can be pruned if the node has lb(L fci ;({j8 (6) })) 
greater than the distance of the lowest examined leaf node. 

Search begins at the root node where each fif* = [1, • • • - !]• We applying the 
ordering heuristics above to move to a node where one breakpoint variable has been fixed. 
From that node, we move down the tree, fixing a value at a time until we reach a leaf. We 
then progressively backtrack up the tree, eliminating branches whose lower bounds are greater 
than the best leaf examined. 

Once the optimal set of buyer breakpoints has been determined, we search exhaustively 
over the supplier breakpoints feasible with this choice according to Equations (26) and (27) . 

Improvements 

We have not yet exploited the fact that, from any node having {ffl}, the constraints in 
Equation (19) may narrow the domains of the buyer breakpoint variables deeper in the tree. 
The computational cost of narrowing the breakpoint domains may be warranted if the pruning 
thus obtained is sufficient. 

E Phase 2 Optimization 

In this section we detail an efficient algorithm to optimize Eq. (5): 

arg min d x (x(x), *) + d 2 (n) subject to C<® (x) A C w (x). 
Our approach is a marriage of constraint programming and local search. 
E.l Local Search over ^-neighborhoods 

We begin with a definition. The d-neighbors of a discrete configuration x are all those con- 
figurations, V, which differ from x in d attributes. Attribute % differs in two configurations 
if Hi ^ Thus the distance varies from 0 if the configurations are identical, up to rid if all 
attributes are different. The ^-neighborhood of *c is the union of all d = 1 neighbors, d = 2 
32 We might also take the minimum over all possible completions lb A ({y3 (6) }) = min (e[1))/S cb) n lb(Ljfe,i({^ (i,) })). 
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(a) 



Table 5: Configurations examined in the D = 3 neighborhood search around (a) the initial 
configuration 11111 and (b) the initial configuration 0+0+0 (duplicate search configurations 
are indicated in bold) . 

neighbors, up to d = D neighbors. The D-neighborhood around x is indicated by J\f D (x)- We 
note that, for any configuration, its n d neighborhood is the entire discrete space. 

Using the D-neighborhood, we apply a standard hillclimbing algorithm which respects all 
the discrete constraints. D is specified before the algorithm begins and serves as an efficient 
5 manner to trade off running time with the quality of the final solution. The hillclimbing 
algorithm takes as input a starting configuration x t and finds the configuration x m e A/"d(x 4 ) 
which minimizes di(x(x), x) + d 2 (x) subject to C^(x) A The algorithm then iterates 
again starting from x t+ i. The hillchmbing algorithm terminates when it finds a configuration 
x T which has a distance less than or equal to all of its £>-neighbors. If D = n d , this is 
10 the global minimum. Given a function which exhaustively searches a ^-neighborhood for 
the lowest distance feasible configuration, a ^-neighborhood hillclimbing algorithm is easily 
implemented. 

There is one noteworthy point about the implementation of an efficient ^-neighborhood 
hillclimber. Consider a D = 3 search over n d = 5 discrete variables. Without loss of gen- 

15 erality, we label the initial configuration as 11111. In the D = 3 neighborhood search, the 
configurations examined are given in Table 5(a) as indicated by the *s. Let us assume that the 
best neighbor in the D neighborhood of 11111 differs in components 2 and 4. We indicate this 
best neighbor as 1+1+1, which serves as the starting configuration for a new D = 3 search. 
The configurations examined in the second neighborhood search are given in Table 5(b). Prom 

20 Table 5(b) we see that the search branches 2, 7, and 10 are redundant, having already been 
performed in the D neighborhood search around lllll. 33 Thus these branches may be pruned 
in the second D neighborhood search. A search branch covers a previous search and may be 
pruned if the vector of positions of the wildcard characters covers the vector of positions at 
33 The duplication of search effort is most acute when D is large and the fittest neighbor is found very close 
by (small d). 
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which the previous best solution was found. The binary relation covers is defined as 

yt covers ^ ^ G x V x- G x'. 

The covers relation is transitive but not reflexive. Thus, to increase the efficiency, we can keep 
a list of the previous branches searched and ignored search branches that have already been 
covered. 

E.l.l Improvements 

Throughout the present discussion we have assumed that D remains constant through the 
hillclimbing process. There are good reasons why this should not be the case. With any 
reasonable landscape correlation structure, we know that lower distance neighbors are more 
likely to be found early in the hillclimbing process. Thus a reasonable heuristic is to start 
with a larger D early on and gradually lower it with time. 

E.2 Exhaustive Search of a D-neighborhood 

To complete the discussion of the phase 2 optimization, we describe a computationally efficient 
way in which to find the lowest distance feasible configuration within a ^-neighborhood. We 
use a branch-and-bound method with heuristics tuned for the present problem. The factors 
tuned include the variable order in which we search, the order of the values any given variable 
may take, and a heuristic which says when certain branches of the search space may be safely 
ignored. 

We use heuristics which early on in the search generate low-distance configurations, and 
later prune parts of the search that we can prove must have higher distance than the best 
configuration generated so far. 

E.2.1 Heuristics And Lower Bounds 

We focus on obtaining good heuristics for the discrete contribution to the total distance, then 
later incorporate the continuous and range contributions since these are simpler. 

At any point in the search, some variables have been assigned and some have not. We 
divide up the sum over the Zij{xi,Xj) into four contributions 

Q d Z v ,{x) = Q d H WiWjZiji^x^ + Qd ^ mwjZiji^Xj) 

i locked j locked i locked j floating 

i floating j locked i floating j floating 

= Z M + + Z f ,i + Zf, f 

Since we are minimizing and all Zi tj (>Ci, xj), «;», and Q d are positive, we can prune any search 
node which has Zy > than the distance of the lowest distance configuration found thus far, 
since additional terms will only add to the sum. 
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We obtain a lower bound on the distance by lower-bounding each of the above four terms. 
Just taking the lowest value in the tables can yield differing values for any given variable, e.g., 
the values for a single variable ^: 

>c* = arg min Z^-fa,^), x* = arg min Z jt i(xj, x;), and x* = arg min Z u {ni) 
may all be different. To correct this we write 

QdZ w = Z^ + Z ti t+ Y Z§{x*) 

i floating 

where and Ztj are calculated as before and 

Z$(x i ) = Q d (u$Z i , i {x i )+w i Y v>j{ZiA*i>"j) + Z jA*j>Xi)))- 

j locked 

The best choice for the floating variable xj is 

x* = argmin Zf?(^)- ( 34 ) 

In this way all the possible choices for Xj are at least consistent, and the true lower limit is 
still guaranteed to be higher than this (since the Z u contribution can only be higher). 

This approach also naturally suggests a variable ordering and a value ordering within any 
variable. We can try the variables in increasing order of their corresponding Z${x-) so that 
variable t, which has the smallest is tried first. For any given variable x* we can also 

order the possible value assignations for X, according to the array of values Z$(xi) so that 
the value Xj = x| is tried first and the value which maximizes z[ l J(tq) is tried last. 

With little extra work we can impose a consistency requirement raising the lower bound 
for Z f)f . We decompose the sum as follows 

Zf f = Q d Y w iWj Z iJ (x i ,K j ) = Q d Y, + 

i,j floating:.?^ %3 floating:j<i 

Instead of taking the minimum component- wise as previously suggested, i.e instead of 

#f,f = Qd Y ™£ i w i w 3 Z ij X j)] , 
i,j floating:;^ 

we take the minimum so that 

Zi f = Qd Y rmnlwiW^Zi^^ + Zj^Kj^i))}. (35) 

i,j floating: j<i 

which enforces pairwise consistency between variables % and j. 
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The lower bound, lb, we have obtained is thus 



lb = Z l } + Zl f + Y^Z®(xt) 

where x* is determined from Eq. (34) and Z£ f is determined from Eq. (35). 

Integration Including Continuous and Range Contributions 

5 In this section we summarize the complete heuristic specification, including the effects of the 
continuous and range variables. As above, we use the definitions 

i,j locked 
j locked 

10 * J floating: j<i 

so that a lower bound to the complete discrete contribution to the distance is Z\^\ + Z% { + 
Y^i -^i*f(^*) where x* = arg min H . Z$(j<i). 

We now extend this result to include the continuous and range variables. The continuous 
and range variables depend on the discrete variables through the mask. Since we expect there 
15 will be relatively few masks we exhaustively enumerate all the discrete possibilities. Conse- 
quently, we define d\ (.^floating Nocked) as th e distance contribution made by the continuous 
and range variables (the d\ contribution to Eq. (5)). We define the continuous + range 
contribution for any given floating variable i as 

d l\xi) = n mi 11 , d l (^floating Nocked) 
floating \i 

20 where the min is taken over all floating variables other than i and where floating variable i is 
locked to value 

Since each continuous/range/discrete distance contributes additively to the total distance, 
we can rank the possible values for any given floating variable i according to Z^(>q) + d^{vj). 
The best (according to this heuristic) setting for floating variable i is 

25 *? = arg min {z$(>a) + 4° («<)}■ 

The heuristic then suggests a variable ordering in order of increasing Z^(x*) + d±\v*) so 
that we first select the floating variable with the lowest sum and last select the variable with 
the highest sum. 

Finally, we turn to the lower bound on the distance. Noting that the continuous and range 
30 distance contribution does not have the form of a sum over the dimensions (as the discrete 



77 



contribution does), we write the lower bound, lb, as 34 

lb = Z x , x + J2(z®te) + + ^f*f 

where i* = arg min^ (^(x?) + d^K*)) (i.e., i* labels the floating dimension with the least 
zfl + d^ 1 contribution to the distance). 

E.2.2 Improvements 

We point out that an extensive literature exists on parallel branch and bound algorithms. 
Many of the techniques developed are applicable to the present problem. 



l 5 x , x < is the Kronecker delta function which is 1 if and only if x = x' . 
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Architecture 

with Testing Interfaces and Testing Systems Shown 




Figure 1: An architecture for the invention. 
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Figure 3: Schematic of a supplier-specific preference with examples indicating potential input. 
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A Method and System for Discovery of Trades Between Parties 



1 Related Applications 

This application claims priority to provisional application no. 60/168,754 filed on December 
6, 1999, titled, "An El-Commerce Infrastructure for Value Chains", the contents of which 
are herein incorporated by reference. This application also claims priority to provisional 
application no. 60/194,880, titled, "Method and System to Mediate Commerce", filed on 
April 6, 2000, the contents of which are herein incorporated by reference. 

2 Field of the Invention 

The invention relates to a method and system for discovery of trades between parties. In 
particular, the invention is a system which allows buyers to define their preferences and sellers 
to define their capabilities, then determines which trading points maximize the utility of the 
buyer. The system suggests trades by exploiting the flexibilities and tradeoffs encoded by both 
parties, thus providing win-win trades. A second level of optimization ranks the trades with 
all suppliers, allowing the buyer to rapidly determine the best alternatives. The system allows 
for rich negotiation spaces and supports continuous, discrete, and range or interval decision 
factors. 

3 Background of the Invention 

The present invention relates to methods of automatic exploration and exploitation of the 
flexibilities possessed by negotiating parties to uncover improved win-win agreements. The 
invention describes computationally efficient mechanisms that are applicable whether there 
are one or many selling parties. The precise number and types of negotiating dimensions are 
irrelevant as long as they are numerical. Thus the present invention applies equally to the 
optimal determination of terms in the purchase of a commodity or an arbitrarily complex 
artifact. 

The past 5-10 years have seen remarkable growth in software tools to help firms with 
enterprise- wide planning (ERP software) and supply chain management (SCM software). 
While these tools do a wonderful job at integrating disparate data sources within and between 
firms, the opportunity exists for significant further cost reductions. 

1 



The same time period has also seen a tremendous rise in the widespread use of the 



internet by both consumers and businesses. Forecasters are predicting that within a few 
years e-commerce between businesses (B2B) and between consumers and businesses (B2C) 
will grow to in excess of a trillion dollars per year in annual revenues. 
5 Electronic markets have proliferated over the last few years with the advent of B2C 

(business-to-consumer) and B2B (business-to-business) electronic commerce. Such market 
places have yielded significant cost savings by lowering the transaction costs between buy- 
ers and sellers. Buyers have also profited through increased competition between suppliers. 
However, electronic markets have hurt suppliers, since the zero-sum negotiation over price has 

10 been at their expense. The present invention describes a tool whereby cost savings for both 
parties are derived from the discovery of win-win trades. Fundamentally, the system works 
by allowing trading parties to describe their desired trade across multiple dimensions and to 
express their flexibility around this ideal trade. Through an algorithmic exploration of their 
flexibilities, the present invention can discover trades that are near the ideal trades of both 

15 parties, enabling both to win. 

The adoption of B2B and B2C electronic commerce was facilitated by the migration of 
catalogues online. This familiar method of presentation ameliorated the significant cultural 
change to electronic trade. For the foreseeable future, electronic commerce will be dominated 
by online catalogs. At present, online catalogues are direct translations of their hardcopy 

20 counterparts where the attributes of a product are described and a price quoted. Inevitably 
however, online catalogs will become more expressive. Catalog entries will be able to represent 
price breaks for large quantity orders, lot sizes, etc. Thus it is important that any software (like 
the present invention) that uncovers mutually beneficial trading scenarios is able to operate 
with such catalogs. Consequently, in the present invention there is an asymmetry between 

25 buyer (usually a human) and seller (usually an online catalog). 

One of the reasons catalogs have come to dominate electronic commerce is that the types 
of goods that can be represented in catalogs are simple. Whether the product is pens or paper 
clips, different vendor's offers differ little from each other (a pen is a pen is a pen), and a quick 
scan of a catalog gives a buyer enough information to make an informed purchase. These types 

30 of goods are low margin and inexpensive. In contrast, the vast amount of purchasing between 
businesses involves materials which are directly connected with business operations - car parts, 
turbines, etc. Such direct goods are the future of electronic commerce. Unlike present-day 
engines, any truly useful procurement tool must be able to support direct materials with 
complex attributes and complex inter-relationships between its components. 

35 Electronic commerce offers unprecedented opportunities for more informed decision-making 

for both buyers and sellers. The past few decades have seen the widespread adoption of enter- 
prise resource planning (ERP) systems, to the point that now almost every major company 
has some form of ERP software. ERP functions as the digital nervous system of a company, 
transmitting and logging information between the company's many different business func- 

40 tions. ERP software keeps track of inventory, monitors the state of purchase orders, signals 
when a company should reorder direct and indirect materials, and a myriad of other functions. 
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Consequently, ERP databases are a rich source of information to optimize a company's oper- 
ations. Yet today this information is rarely used to make more informed buying and selling 
decisions. The present invention can utilize such information sources to optimize a company's 
interactions with suppliers and customers. 
5 One important manner in which this optimization can occur is through an analysis of all 

cost factors. Current buying and selling practices often focus on limited goals, e.g., minimize 
the total purchase price. Myopic purchasing strategies often result in higher total cost of 
ownership when all cost factors relevant to a product in its lifetime of use are included. These 
other cost factors can be significant. Why save the money in taking delivery two days late 

10 if the receiving docks will be full at that time and an additional shift needs to be hired to 
clear the docks? Why order the cheaper drill bit if it is much more expensive to replace when 
it breaks? The present invention improves trades by minimizing the total cost of ownership 
of a product, yielding significant savings to its users. Many total cost factors are difficult to 
quantify - e.g. what is the cost of dealing with a unionized versus a non-unionized supplier? 

15 Consequently, the present invention supports qualitative (best guesses and intuition) as well 
quantitative factors. 

All companies are situated in a supply or value chain. At each step in the chain, a 
company purchases from its suppliers, transforms these inputs, and sells the output to its 
customers. The termination of the supply chain is the sale of the final product to the end 

20 consumer. Since the only influx of external capital comes from the end consumer, companies 
have realized that they compete not only as individuals but also as entire supply chains. 
As result, software products have recently become available which attempt to streamline the 
operations of links within the entire supply chain. This software, variously called supply chain 
optimization (SCO) or advanced planning and optimization (APO), operates on the basis of 

25 forecasted demand at various points within the supply chain. Based on these predictions, plans 
are generated telling companies how much to produce and how to schedule their operations. 
SCO systems are a valuable source of intra-company information - data the present invention 
capitalizes on. Because SCO software relies on forecasted demand, it is only as helpful as 
the forecast is accurate, and, unfortunately, in many cases demand is very difficult to predict. 

30 How can the software know that laundry detergent will go on special at grocery stores in the 
Northeast in 7 weeks? As a result of the volatility in demand and the many other unpredictable 
perturbations that plague supply chains, companies keep significant buffers in the form of 
inventories. In addition to planning, businesses must also be able to adapt to unplanned 
effects. Such adaptation requires flexibility and a means to exploit that flexibility. The present 

35 invention exploits the flexibility of trading parties to streamline the operations of supply chains 
by smoothing the boundaries between trading parties. 

The present invention is therefore a system to allow trading parties to express trading 
desires and constraints across many and varied different factors. These trading preferences 
are informed by many different data sources to optimize for a company's internal operations 

40 and its connections to it's supply chain through an analysis including total cost factors. The 
flexibility expressed by all trading parties is exploited to locate win-win opportunities for all 
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parties if they exist. 

4 Summary of the Invention 

We describe the present invention in its application to facilitating trade between buyers and 
sellers, but note that the mechanisms described are much more general. We can easily imagine, 
5 for example, using the present invention to match individuals (with the desires and skills) to 
projects. 

The inspiration for the present invention comes from utility theory developed by economists 
since the 1960's. Since we are interested in multiple dimensions of negotiation, we draw from 
the multi-attribute utility theory literature. 1 Utility is an abstract concept which has been 
10 formalized in various ways. For the present purposes utility, u, is a number between 0 and 1 
representing a party's willingness to trade. Larger values indicate a greater willingness. 

4.1 The Negotiation Space 

In any negotiation the parties must come to agreement on the factors requiring negotiation. 
We call these factors dimensions or variables. As an example, when purchasing a car, the buyer 

15 may be concerned with price, time of delivery, and color. Each factor price, time, and color 
is a dimension. Most dimensions can be classed as one of three types: continuous, discrete, 
or range/interval. A continuous dimension is one like price for which the buyer's utility varies 
smoothly across that dimension. The buyer's utility at $23 001.00 is almost the same as the 
utility at $23 000. Color is a discrete dimension. Since the car may only be available in 

20 black, white, and silver, the domain of this dimension is the finite set of values {black, white, 
silver}. Moreover, the buyer's utility may be quite different for the three colors. The third 
class of dimensions is called interval dimensions. An interval dimension arises often in B2B 
negotiations. If a machined part is built to some tolerance (e.g., the inner diameter of a screw 
is between 24.5 and 25.5 mm), the range of variability in the dimension is specified as an 

25 interval. In the language of statistical quality control, a certain percentage of the machined 
parts will fall in this range. These three broad classes of variables capture almost all the types 
of attributes relevant to B2B negotiation. 

The present invention operates over any number of continuous, discrete, and range or 
interval variables. We call the negotiation space X and any point in the negotiation space 

30 (x,xr, r) E X. It is important to recognize that the single trading point (x,>f,r) may have 
multiple components, e.g., price = $23 000, time of delivery = 3 weeks, color = black. 

In the present invention, the space of negotiation is agreed upon by all parties involved 
prior to the commencement of any negotiation. We can, however, imagine more dynamic 
situations in which dimensions are introduced and discarded over time. 
1 For a good introduction to multi-attribute utility theory see [1]. 
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4.2 The Buyer's Utility Function 

A party defines it's utility function over this space so that every (x, x, r) is assigned a util- 
ity number indicating the party's willingness to trade. We indicate the utility function as 
u((x,*,r)). A great deal of work has been done on the appropriate form for utility func- 
5 tions. In the present invention, we take a simple form for the utility function for two reasons. 
First, we would like the form of the utility to be conducive to rapid computation. Sec- 
ond we would like the utility to be simple enough to be easily understood by and elicited 
from users of the invention. With no loss in generality, we write the utility function as 
u((x,*r,r)) = exp(-d((x,x,r))) where d(x) is interpreted as the distance of trading point 

10 (x,>f, r) from the most preferred trade. 

So that we can operate against seller catalogs, only the buyer needs to define a utility 
function. Across the continuous dimensions, the buyer's utility is defined by specifying the 
most preferred (or ideal) continuous dimensions and the manner in which utility drops off as 
we move away from this ideal. For the discrete dimensions, the utility is specified in tabular 

15 form since there are a finite number of alternatives. Again, the buyer must specify it's ideal 
discrete values and how utility decays away from those values. In section 6.1 we describe how 
this is accomplished. The range dimensions contribute to utility similarly; the buyer specifies 
an ideal range and the utility decays for ranges other than the ideal according to their distance 
from the ideal. 

20 The utility function can also express tradeoffs between variables, e.g., I may take delivery 

in 5 weeks if the price drops to $20 000, or I may accept the white car if I can take delivery in 
2 weeks. The tradeoffs may be between pairs of continuous dimensions (as in the first case), 
between pairs of discrete variables, or between continuous and discrete variables (as in the 
second case). 

25 4.2.1 Normalization and Weighting 

When utility is defined over different types of variables, it is important to normalize the 
contributions of each variable so that the buyer can weight the importance of the various 
contributions to utility. This is a difficult problem. How should a buyer's color preferences 
be normalized so that they can be traded off against time of delivery? The present invention 

30 solves this problem by requiring that the average distance of any negotiation variable from 
its ideal value is the same for all dimensions. Since the buyer is more interested in those 
regions of the negotiation space where the utility is high, the average is weighted by utility. 
This procedure defines a manner in which to define a baseline where all dimensions contribute 
equally. Given this baseline, the buyer can then weight the various contributions and obtain 

35 useful results. 

4.2.2 Utility Elicitation 

Since utility is fundamental to the present invention, its elicitation from the buyer is important. 
Utility may be defined using any of a number of sources: 
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1. 



graphical user interfaces associated with the invention 



2. 



standard benchmark criteria applicable to the domain 



3. formal methodologies like the analytical hierarchical process [2], or discrete choice anal- 
ysis [3] 



We expand briefly upon method 4. As discussed in the background section, it is important 
to buyers to minimize their total cost of ownership. If we have a function representing these 
costs as a function of the negotiation variables, and perhaps other factors, this function can 
be used to infer a utility function which will act to minimize the total costs. Later we describe 
10 how this can be accomplished. 

4.3 A Supplier's Capabilities 

As noted previously, suppliers are treated differently from buyers so that the tool can operate 
against catalog information with no human intervention required on the part of the seller. In 
fact, we do not require sellers to define a utility at all. 

15 A supplier cannot offer deals at all points within the negotiation space X, e.g., he certainly 

can't offer the black car tomorrow for free! A capability then represents the ability of a supplier 
to deliver and defines a subspace of X. It can include such things as price discounts on large 
volume orders, variation in delivery time as a function of price, etc. Since these relationships 
are already specified by businesses in terms of simple rules like "the price per unit is $10.00 

20 if 1 to 999 units are ordered and $9.50 per unit if 1000 or more units are ordered", suppliers' 
capabilities are represented in the present invention by piecewise linear functions. 

4.4 Negotiation Constraints 

Both parties may have constraints which must be satisfied in order for them to trade. For 
example, the buyer may not buy the car unless he gets it within 6 weeks, or he may not 

25 purchase the car if it is available only in white. These are examples of continuous and dis- 
crete constraints, respectively. A continuous constraint sets a requirement on the continuous 
variables. In the present invention, continuous constraints must be either linear or quadratic. 
Discrete constraints involve discrete variables. A discrete constraint can be expressed as a list 
of the allowed (or disallowed) combinations of the discrete variables for which the trade will 

30 be acceptable. For example, if the buyer would accept either the black or the silver car, the 
constraint would list both these colors as viable. It is important to note that both continuous 
and discrete constraints may involve one or more variables. We can also express constraints 
involving both types of variables by allowing the continuous constraints to differ depending 
on the discrete variables. 
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inferred through models 
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4.5 Utility Optimization 

With the major components of the invention in place, we describe how the overall system 
works. As a procurement tool for the buyer, there are two levels of optimization. First, for 
any given supplier we maximize the buyer's utility, subject to the supplier's capabilities to 
5 find that trade which makes the buyer as happy as possible. Since we are optimizing within 
a supplier's capabilities, the supplier has expressed a willingness to complete the trade at 
whatever point is determined to be optimal. The tool then optimizes across suppliers to rank 
them according to utility at the optimal point. A graphical user interface allows a buyer to 
investigate the trades suggested by the tool by sorting according to any dimension or by the 
10 overall utility. 

Utility, while a useful concept in assessing an overall score, may be of limited use to a buyer 
due to its abstract meaning. Consequently, we can also apply the total cost of ownership 
function to the results to rank order the suggested trades according to their various cost 
components. Recall that for any trade x e X, the total cost of ownership function returns 
15 the various cost contributions. This additional information aids the buyer in his purchasing 
decision. The utility number for each trade is still useful because the total cost of purchase 
function includes only those cost factors which can be quantified, whereas the utility also 
includes "softer" qualitative factors. 

4.5.1 Aggregation 

20 In addition to optimizing against one supplier at a time, the present invention can also be used 
to optimize against an arbitrary aggregation of suppliers. This is important if, for example, 
no single seller can supply the large volume requested by a buyer. In this mode of operation, 
the buyer specifies sets of suppliers participating in the aggregation and the dimensions over 
which aggregation can occur, and the tool finds the optimal combination in which to distribute 

25 the volume dimension over the allowed suppliers. 

4.6 An e-Commerce Infrastructure for Value Chains 

This patent application also describes an integrated solution for B2C and B2B e-commerce 
that would be built on top of ERP and SCM software and which would provide a number of 
compelling benefits to companies. Amongst the benefits are: 

30 • multidimensional markets which allow consumers to implicitly define their preferences 

over many criteria. This allows both consumers to express what it is they really value, 
allows companies to position themselves clearly in the space of value, and allows for 
efficient matches between trading partners 

• optional anonymity of market participants and their trading desires when that is appropriate 2 

35 • explicit pricing of the flexibility possessed by the consumer and all businesses in the 

supply chain which allows for more robust operation of the entire supply chain. This 
2 Anonymity is not, however, a pre-requisite of the proposed invention. 
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concept is very different from other types of markets (e.g. auctions, reverse auctions, 
exchanges) where transactions are specified exactly. The flexibility introduced by any 
party, whether consumer or supplier, is propagated and exploited through the entire 
supply chain. 

5 • capture and quantification of true consumer demand leading to improved forecasting 

and product development by suppliers 

• automated markets that integrate supply chain networks through coordination across 
and within company boundaries. Coupling of the automated markets with local (i.e. 
at the company level) optimization tools fed by real-time company data allows for 
10 optimization and cost savings across the entire supply chain- 

It should be recognized that supply chains may be very different in the near future- 
Current supply chains are based on physical objects made valuable through a sequence of 
transformations resulting in a product purchased by an end consumer. With the move to an 
information economy the supply chain of the near future may not involve physical goods at all. 
15 In particular the entire supply chain may consist of value adding operations converting raw 
data to consumer-desired information. Such supply chains will have the same coordination 
problems current ones do. Our proposed solution applies equally well to these future supply 
chains and by supply chain we mean this more general notion. 

5 Brief Description of the Figures 

20 Figure 1 shows an [an] architecture for the invention. 

Figure 2 shows a schematic of a buyer-specific capability with examples indicating potential 
input. 

Figure 3 shows a schematic of a supplier-specific preference with examples indicating potential 
input. 

25 6 Detailed Description 
6.1 Theory 

In this section we outline the mathematical foundations of the optimization process in sufficient 
detail to allow for computer implementation. 

6.1.1 The Negotiation Space 

30 In Table 1 we define the parameters which collectively define the space of negotiation X. For 
each of the n c continuous variables, we specify an allowed range over which that continuous 
dimension may vary as Xi £ Xi = [x^Xi], where x is the n c - vector of lower continuous bounds 
and x is the n c -vector of upper continuous bounds. Each discrete variable assumes a value 
from within its domain k\ E T>i. Without loss of generality, we label the domain of discrete 
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n c 
n<i 
n r 


number of continuous dimensions 
number of discrete dimensions 
number of range dimensions 


X 
X 


n c -vector of values for continuous dimensions 
rid- vector of values for continuous dimensions 


%i 


value of ith continuous variable 
value of ith discrete variable 



Table 1: Definition of the negotiation search space. 



variable i by T>i = [1, • • • , d;] where dj > 0 is an integer giving the number of possible values 
discrete variable Xi may assume. 

With these definitions, we define the space of negotiation by the tensor product X = 
X± <2> ■ ■ ■ <£> X nc ® T>i <g> ■ • ■ ® T> nd . Range variables are treated separately and not negotiated 
5 over. 

6.1.2 The Utility Function 

The utility function is a mapping from X into the interval [0,1]. As indicated earlier we 
assume the utility to have the form u(x,x) = exp[— d(x,x)] where d(x,x) is interpreted as a 
distance. In what follows we will assume that in its simplest form the distance function has 
10 the form 

d(x,*r,r) = (x -/*(*:)) ^(^(x-ziO*)) + Z(x) + R(r; x). 

Each contribution to the distance function is positive. We consider each contribution to the 
distance in turn, beginning with the range variable contribution R{r;x). 

First, we note that the range distance depends on the setting of the discrete variables. 

15 This allows the buyer to express different preferences for the range variables depending on 
discrete factors. The total range distance is summed up over all possible range variables 
so that R(r;x) = Yh=i R i( r i^ H )- Tne vector r indicates the preferred values for all range 
variables. If range variable i is specified as the interval n = (24, n) (where > r,) then r is 
an n r -vector of such tuples. The distance contribution, i?;, from one range variable will depend 

20 on the application. If the range variables are meant to represent the tolerances on machined 
parts where issues of statistical quality control are important, then the distance between two 
ranges might be related to the overlap between Gaussian distributions. If n is interpreted 
as a Gaussian having mean (r^ + fj)/2 and standard deviation proportional to n - then 
an appropriate range distance is given in Appendix A. Other choices for the range distance 

25 function are certainly possible. 

The continuous distance is quadratic and determined by the positive semidefinite n c x n c 
matrix C _1 . We have allowed this matrix to vary with the setting of the discrete variables and 
indicated this explicitly through C _1 (x). 3 The n c - vector fi may also depend on x and indicates 
the point at which the utility is maximal - fi is thus identified with the ideal value for the 

30 continuous variables. The precise quadratic form is convenient, but, using recent developments 
3 Where no confusion will arise we eliminate the dependence of C _1 on yc for notational simplicity. 
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in interior point methods, other convex functions are also computationally tractable [4]. 

The discrete distance is determined through the function Z(*c) which maps the discrete 
space V\ ® ■ ■ ■ ® V nd onto the positive real line [0, oo]. In keeping with the assumption that 
distance is a function of only pairs of components Xi, xj, we assume the discrete distance has 



Each contribution Zij is a table consisting of didj entries, where Z^{xi, x 3 ) can be interpreted 
as the distance if discrete dimension i has value Xj conditioned on discrete dimension j having 
value Xj. The diagonal terms Zi offer an unconditional distance. The most preferred value 
for the ith discrete dimension is that for which Zi(xi) = 0. 

Rather than require the user to enter the distances explicitly, there are numerous ways in 
which the distances can be generated automatically based upon a buyer's ranking of preferred 
values. Further details can be found in Appendix B. 
Weighting of Dimensions 

In many cases it is important for simple modifications of the distance function to re-weight 
the contributions to the total distance. If w c is an n c -vector of weights for the continuous 
dimensions, we can accomplish this by letting C" 1 = W C C _1 W C where W c is the diagonal 
matrix W c = diag(w c ). 5 In a similar way we modify the discrete distance to Z w -ij(xi, Xj) = 
w d .iW d .jZi t j(xi,Xj) where w d is the rc d -vector of weights for the discrete variables and w d;i 
is its ith. component. The range contribution is also modified so that R w; i(ri) = w r -iRi(ri) 
where w r is the n r -vector of weights for the range variables and w r -i is its ith component. For 
convenience the weights are normalized so that (l*w c ) 2 + (l*w<£) 2 + l*w r = I. With little 
additional complexity the dimension weights can be made dependent on the setting of the 
discrete variables but we will assume throughout that the weights are constant. 
With these modifications, the total distance function becomes 



where Z w (x) = J2i=i ™d;i{ w d;i z i(*i)+T,%i&i) ™d;jZi,j{^ and i^(r) = ££i «V ; ii2tfo) 
Assigning weighting factors is useful only if the relevant contributions have been previously 
normalized so that they are all roughly the same magnitude. This serves as the baseline for 
which all weights are equal. The question immediately arises as to what criteria to use to 
weight the distance contributions. 

We shall determine scaling factors, Q d > 0 and Q r > 0, so that the average distances per 
dimension of the discrete, range, and continuous contributions are equal, where by average we 
mean a utility-weighted average over the entire space of possible trades. This weighting places 

more emphasis on the better trades 

If d c , d d , and d r are the continuous, discrete, and range contributions to the total distance, 
then after multiplication by the scaling factors d = d c + Q d d d + Q r d r . The scaling factors are 

4 Later we shall generalize this distance to include weighting of dimensions. 

5 M = [m lt j] — diag(v) is the diagonal matrix formed by setting m;,; = Vi and rm t j — 0 for j ^ i. 



the form 4 




i=l j=l(&) 



d(x,x) = (x-/x(x)) t C- 1 (> < r)(x-/i(>«:)) + Z w (x) +i? t0 (r). 



(1) 
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determined through the utility weighted average distances denned by 

_ £ x £„ J v du Qrdr exp[-Q r d, - Q d d d - d a ] _ £„ ex P [-Q d d d ] £, d, exp[-Q,d„] JV du exp[-d e ] 
<<W = £„ / v du «p[-Q,d, - C?^ - d.] W ' E„ «P[-«d««,l] £„ exp[-Q.d.] J v du exp[-d c ] 

S« £r JV du Qd^ exp[-O r d,. - Q d d d - d c ] _ „ S„ dd expI-Qddrf] £ r e*p[-Qrdr] Sv exp[-d e ] 
<di> = E r / v du exp[-Q.d, - Q d d d - d c ] Qd £« exp[-Q d d d ] E r ex P [-Q r d r ] / v du exp[-d c ] 

_ £ r JV du d c exp[-Q r d r - Q d d d - d c ] = exp[-Q d dj]i: r exp[-Q r d r ] JV du d e exp[-d c ] 
(cJc> - E v E, JV du expt-Qr-d,. - Qdd^ - d c ] ■£„ exp[-Q d d d ] £ r expt-C^d,.] f v du exp[-d c ] 

A few comments on the above equations are in order. First, indicates the repeated sum 
' • ■ YLx over a11 possible discrete trades. $2 r indicates a sum over all the range variables 
and the integral over volume V indicates integration over the continuous trading volume of 
interest. Finally, we have not included a scaling factor Q c on the continuous distance, since 
this can be made equal to 1 if we reinterpret Q r as Q r /Q c and Q d as Q d /Q c - e Each of the 
averages is an explicit function Q d and Q r . 

The requirement on equal average contributions determines the two unknowns Q r and 
Q d through the equations: (d r )/n r = (d c )/n c and (d d )/n d = {d c )/n c . These two nonlinear 
equations are coupled in terms of Q r and Q d and must be solved simultaneously for Q r and 
Qd. Further details are found in Appendix C. 

6.1.3 Constraint Specification 

Buyers and sellers may express constraints over both continuous and discrete variables. 
Continuous Constraints 

For simplicity (and because additional expressiveness is rarely required) we assume that the 
buyer's constraints over the continuous variables are linear. 7 This allows a buyer to express a 
constraint, e.g., the time of delivery must be within 10 days or I will not trade, i.e., t < 10. 
We allow for both inequality and equality constraints which can be expressed as G^ 6) x = g^ 
and G^ 6) x < g 2 b) . If there are cf ] equality constraints then has cf } rows. Similarly, G^ 6) 
has 4 6) r o ws if there are cf> inequality constraints. We allow the constraints to depend on the 
setting of the discrete variables, and to be explicit we often write G^(x), g£ 6) (x), G 2 6) (x), 
and g<»(x)- 
Discrete Constraints 

We use a standard methodology to represent and process constraints over discrete variables [5]. 
Abstractly, a constraint over a (perhaps proper) 8 subset of the discrete variables is represented 
as a list of all the allowed combinations the variables may assume. An example representation 
of a pair of discrete constraints is given in Table 2. There are two solutions to this set of 
constraints: xi = 1, x 2 = 2, x 3 = 3 and xi = 3, x 2 = 2, x 3 = 1. We indicate these solutions 
as {(x 1? 1), (x 2 , 2)}, (x 3 , 3)}} and {(xi, 3), (x 2 , 2)}, <x 3 , 1)}} respectively. Each solution where 
all the variables have been identified with specific values from their domains is called a labelling. 
Computationally efficient representations are used to ensure that only feasible combinations 
6 We singled the continuous variables out, since there will always be continuous variables in any trading 
scenario. 

7 With little additional complexity we can also handle quadratic constraints, see [4]. 
8 A proper subset of a set is the set itself. 
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X2 


X3 


1 


2 


3 


1 


3 


2 


2 


1 


3 


2 


3 


1 


3 


1 


2 


3 


2 


1 



ill 
2 



(b) 



(a) 



Table 2: An example set of constraints involving 3 variables where the domains of all variables 
are V\ = T> 2 = T> 3 = {1, 2, 3}. Constraint (a) requires that the values assumed by xi, x 2 , and 
x 3 are all different from each other, and constraint (b) requires that the value assumed by x 2 
is even. See text for the solution to both constraints. 

of variables are ever processed. Numerous third-party libraries offer constraint programming 
functionality. 9 

6.1.4 Utility and Total Cost of Ownership 

The buyer's utility function and associated constraints may be difficult for many users to 
5 define. In this section we show how models of the buyer's business can be used to define 
utility in a natural manner. 

We imagine a function which provides an estimate of the total cost of ownership for any 
given purchase. Cost contributions to this function might include piece part costs, freight costs, 
setup costs, quality assurance costs, repair costs, etc. It is important to include all quantifiable 

10 costs associated with the lifetime of use of the purchased product because it is this function 
we will be minimizing. Significant savings may be obtained by taking a longer-term view 
of the purchase. Revenues (negative costs) generated from the purchase are also included 
in the function so that the function represents some measure of profitability associated with 
the purchase. We write the total cost of ownership function as (7 o (x,**,r;0). We explicitly 

15 indicate the dependence on the negotiated trade parameters x, *r, and r, as well as other 
factors f$. The other factors might include forecasted demand, current inventory levels, etc. 
These factors will vary over time, and they can be extracted from the buyer's ERP and supply 
chain management systems (SCM) in real-time just before the purchase to ensure continuous 
real-time optimization. See section 6.2.1 for further details. 

20 Minimization of C 0 (x,>f, r;/3) defines an ideal trade dependent on current conditions: 

Xo P t(£), *o P t(|8), r op tOS). If desired, these can be used to define fj, = x opt (/3, r = r opt (/3) 
and the desired ideal discrete configuration x op tW) (having distance contribution Z = 0). 
Moreover, the tradeoffs between continuous dimensions around this minimum can be obtained 
through calculation of the Hessian matrix H = [hij] where the i,j matrix element is given by 

25 ^ ~ dXjd Xj lx=x opt (/J),,,=><„ pt (/J),r=r opt ( / 8) 

9 See for example www.mozart-oz.org or www.ilog.com. 
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We then identify C" 1 with H. In this way, little trading flexibility is obtained in directions 
where total cost of ownership rises rapidly, while significant flexibility is obtained in directions 
where total cost of ownership increases slowly. 

In summary, a total cost of ownership model defines both the most preferred trade param- 
eters and the flexibility possessed around the preferred trade. The model pulls dynamically 
from real-time data sources to provide the most up-to-date optimization based on total costs 
of ownership and other important qualitative factors the buyer may wish to describe in the 
utility function. The same function and its constituent costs may also be used to help analyze 
proposed trades from suppliers. 

6.1.5 Supplier Capabilities 

As discussed in the introduction, suppliers represent their capabilities through a specification 
of the subspace of X in which they will trade. A supplier's capabilities must specify the 
allowed continuous, discrete, and range variables. The allowed range variables are expressed 
as the pairs (r_j,rj), one for each range variable. For example, if a supplier produces 25mm 
inner diameter screws to within a tolerance of 0.5 mm, then the range variable is simply 
(24.5,25.5). These are compared with the buyer's ideal range and contribute to the distance 
function through the R w (^) function. 

Capabilities over continuous and discrete variables are more complex. 

Continuous Capabilities 

Continuous capabilities are viewed naturally as responses to a buyer's request. Thus we 
distinguish between a buyer's requested continuous vector and a seller's response x.^ . A 
vector- valued function, f (x^jX^U) returns the response based on the buyer's request and 
also, perhaps, other previously defined supplier responses. Component fa of f defines the «th 
continuous variable, i.e. xf* = /i(x^,xW). 

Currently, suppliers are used to quoting price discounts for large volume orders and these 
price discounts are expressed as piecewise linear functions. Consequently, we restrict fa to 
have the following form (where we distinguish between the functions depending on the buyer 
and seller variables): 

k 

An example of how this may be used to define a supplier response is the following: We 
assume three continuous dimensions - price, volume, and time of delivery and indicate these 
as [xi,X2,xz] = \p,v,t]. Then a response may be formed as 

= 

tW = / t (t;W,^) = /*,.(««) 

p (s) = / P (« w ,* w ) = f P Av (s) ) + UAt [s) )- 

The function returns the volume a supplier will fulfill as a function of what the buyer 
asked for. If the supplier can deliver any volume, this will be the identity function. If the 
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supplier delivers only in certain lot sizes, this function may have a staircase shape, etc. The 
f tv function indicates the time it will take a supplier to deliver a certain volume. So, for 
example, if larger shipments require longer transportation, then this dependence is given by 
this function. Finally, we turn to the price determination. In this example the price depends 
on the quantity «W being shipped and the f PjV might represent price discounts for large volume 
orders. There is also an incremental price contribution based on the time of delivery. If faster 
delivery is more expensive, this is reflected in f p>t . 

For a given setting of the discrete variables, each f$ (x^\^) and f^(x^\^)is 
a one-dimensional piecewise linear function. Consequently, the functions can be specified 
by listing the breakpoints. If f$ (x^ } , ) has «jj breakpoints, then we list these as 
{ x k ] ( b i S k = 1 )' a; i S) ( & i S 2 = 2 )>"' '^fc^ifc = K< fk)Y° and function values at these break- 
points 'are {f$(xplb% = i)U$(xPib# = 2)),... ,/<fM } (fcg = fM- S j milarl * 

ffl, is a piecewise linear function denned by the nfl breakpoints {xf\b^ ] k = l),xf ] (bfl = 
2), • ■ ■ , 4 b) (&g = «S2)} and function values at these breakpoints {/$ (zf (6$ = 1)) , (*f (6$ 
2)), • • • , f[ b l (x^ibfl = «$))}. The breakpoints are indexed by the integers 6$ and b\% 

An interval is specified by assigning a value ftjj G [l,^ _ *] and b tl 6 ~ ^ so 

that 11 

4 S) (©<4 S) <4 S) (C ) + 1 ) v * and 4 6) (^)<4 6) <4 6) (^ + 1 ) v - < 3 ) 

Within each interval the functions are linear, so we have 

where Ci(vW, {&$}, {&$}) = £ fc c$0* (s \ &$) + <$(6$). In the above equations, the inter- 
cepts and slopes are given explicitly by 



4" cff + D/g (4" < fc .,*» - 4" (M/fl (4" + d) 



(.)„(.), _ /g(4"("a + i))-/g(4 ,) (M) 
™*<W- 4»' ( 6,!;» + i)-4"(^) 



respectively. An analogous result holds for the cjj and 

To eliminate any cyclic dependence on xf^ we must impose an ordering on df* so that 
can only depend on x^ where j < i. Consequently, we can write 

^u Ci( ^)a6S^...,ei}.{^}) + E^('' w .^)4 s) +E^(€)4' ) - 



10 The breakpoints axe assumed to be in increasing order. 

"Note that many choices for or {&^} would be inconsistent with these constraints. A simple way t 

fix this is to define a union of all breakpoints (s) and and have all fi, k have the same breakpoints. 



Written as a matrix equation, the above becomes 



(I - MW(*W, {6g}))xW = c(xW, + M«({(|})xW 



where c(xW, {&«}, {6$}) = [d(^), {6^}, •■■ ^(xW, {&«}, {&$})] \ 

*«=[*<•> ••• zfe>] ',*<*>=[*?> ••• *£>]', and 



0 0 ••• 0 

MW(xW,{&W})= m&(*M,6g) 6$) •■• 0 



0 



0 



0 
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m 



and 



M«({i>< 6 >}) = 



•4MJ) »4>S) 





In most cases will depend only on a subset of the variables in x( 6 ) . If x^ depends on 
n' < n of the x( 6 ) variables, then M( 6 ) is an n x n' matrix. In the example given everything 
10 depending only upon the volume the buyer requested. 

Since M^^x^)) is lower triangular and can be inverted in time 0(ri), we can rapidly 
express x( s ) as 



x(*> = (I-M( s )(*( s ),{&^ 



15 will be used in section 6.1.6 which formulates the optimization problem. 

We also allow a supplier to express additional linear constraints so that, for example, 
he may represent that he does not deliver on Sunday. Thus the supplier may define the 
matrices G^(x), G^\x), and the vectors g]^(x), g^C^) such that Gj^x( s ) = g^ and 
G^x^x < g^. G^(x) and G^x) have c[ s ^ and rows respectively. 

20 Discrete Capabilities 

It is easy to imagine that a supplier's response on a discrete dimension is highly constrained 
by the values of the response on other dimensions, e.g., certain product characteristics come 
only in certain colors and package sizes. Consequently, it is not suitable to explicitly define 
a response but only to make available a supplier's constraints amongst the discrete variables. 
25 Consider 3 discrete dimensions where *c\ G T>\ = [a,b,c], >c<2 € T>2 = [A,B,C,D], and >cz G 



(4) 



as 



long as the 6^2 are chosen to also satisfy x^\b[ s l) < < x^(b^ + 1). These constraints 
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V3 — [a, f3, 7, 5], and assume the supplier has the following 3 constraints 



Ci m) = {(«, a), (a, 0), (b, /?), (c, /?)}, 
C 2 (^2,^3) = {(A 7), (A/?)}, 
5 C 3 (^i) = {&,c}. 

We first note that there are 4 feasible solutions (or product configurations the supplier can 
meet): [x u x 2 ,X3] = [b,A,0\, [b,D,0\, [c,A,p], or [c,D,0\. Feasible solutions to the con- 
straints define the response for the discrete variables. 

We indicate a supplier's or buyer's collective set of discrete constraints by C^{x) and 
10 C^ h \)€) respectively. 

6.1.6 The Optimization Problem 

Having defined the necessary components, we now define the optimization task which deter- 
mines the continuous x* and discrete *c* parameters of the trade. 

Since the trade must be acceptable to the supplier, we maximize the buyer's utility over 
15 a supplier's capabilities. Equivalently, we minimize the distance from the buyer's ideal values 
as 

[x*,x*] = arg min {(x<"> - /i(x^))*C- V s) ) (* {s) - /*(* W )) + QdZ w (x (s) ) + QrRwW} 
where 

x« = (I - MW(xW)) _1 c(^) + (I - mW(xM)) _1 mWx< 6 > 
20 subject to the constraints over continuous variables 

Gi(^)x w = gi(^ (s) ), G 2 (^ a ))xW < g 2 (* W ) 

and the constraints over the discrete variables C^(v^), C^(y^). In the above, we have 
defined the (cj s) + cf ] ) x n c and (d 2 s) + dp) x n c matrices Gi(x^) and G 2 (^ (s) ) by 

^, =[};<:;:;], - o,^ 

25 The (4 s) + 4 6) )- and (4 s) + 4 &) )-vectors Si(^ (s) ) and g20* (s) ) are defined by 

The optimization is accomplished by iterating two distinct phases. Phase one sets the 
continuous parameters optimally for a given setting of the discrete variables. We define the 
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functions 



d!(x,x) - (x - p(x))*C- 1 (x))(x - /i(x)) + J?(r;x) and d 2 (*r) = Q d Z M (^), 
The first phase of the optimization is the continuous problem: 12 

x*(x) = arg min di(x (s) ,x) subject to Gi(x)x = gi(x) and G 2 (x)x < g 2 (x). (5) 

A detailed discussion on the solution of the phase 1 optimization problem is found in appendix 
D. The second phase determines the best value of the discrete variables by minimizing over a 
function of *c alone 

x* = arg min di(x(x),x) + d 2 {x) subject to C (6) (x) A C^ s) (x). ( 6 ) 

Further details on the phase 2 optimization are given in Appendix E. Once x* has been 
determined, we find x* as x* = x(x*). 

6.1.7 Aggregation 

Often a buyer may be willing to divide an order between multiple suppliers in order to aggre- 
gate the required demand or to obtain better deals. In this section, we detail how the present 
invention supports this aggregate optimization. 

Aggregation can only occur over the continuous variables where values may be subdivided. 
Each continuous variable x t must be parcelled out amongst a set of suppliers. Consequently, 
we extend our notation to Xi -)> 5;, k giving the contribution of the fcth supplier to continuous 
dimension i. The fcth supplier may come from a (perhaps proper) subset of all suppliers. 
We indicate the set of potentially contributing suppliers as K, and the number of potentially 
contributing suppliers as |/C|. 13 

We restrict the discrete variables to be the same across all potentially aggregated suppliers, 
i.e., we do not generalize Xi x^. This simplifying assumption is made for two reasons. 
First, the size of the discrete optimization problem is smaller and so optimization be performed 
faster. Second, it may be difficult to elicit from the buyer the allowed discrete alternatives for 
each supplier. Nevertheless, this generalization is straightforward should the need arise. 

This simplifying assumption requires that the union of discrete supplier constraints C/c(x) = 
Afce/C vields a feasible solution when combined with the buyer's discrete constraints 

C^ b \x). A necessary (but not sufficient) condition for satisfaction is then that each constraint 
satisfaction problem h having constraints C W (x) A (x) has a feasible solution. 14 Hence- 
forth, we will assume that the set of suppliers, /C, satisfies this condition. If not, those suppliers 
violating the constraints (x) A (x) are eliminated from consideration in K. 

12 No optimization is required over range variables, since these are specified up front by both buyer and seller 
and merely add to the total distance. 

13 All work to this point is thus seen as the special case \JC\ = 1. 

14 It is easily seen this requirement is not sufficient by having C (6) = {(*n, 1), (xi,2)}, C[ s) = {(hi, 1)}, and 
C« = {(*i,2>}. 
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Discrete Search 

We must search over the subsets of K for feasible solutions, which is a combinatorial problem. 
Fortunately, given a complete labelling of variables, determining the largest subset is easy. 
For any given labelling of all discrete variables, if each A C< s) V k £ k C K is satisfiable, 
then the union C< 6 > A C ( K s) where c£° = /\ keK is also satisfiable under the same labelling. 
The largest subset of variables is found by adding all k which have feasible solutions with 
the buyer. We needn't worry about smaller subsets because the continuous optimization will 
assign zero values to those if appropriate. Consequently, for any given labelling x we let £(x) 
represent the maximal subset of suppliers for which cW(x) A C K (x) is satisfiable. It is this 
set of suppliers which enter into the continuous optimization. The number of participating 
suppliers is denoted by |/C(x)|. 
Continuous Optimization 

Optimization over the continuous variables is carried for each labelling x. Generally speaking, 
the buyer's utility will not be an explicit function of but only of a*. We assume a linear 
relationship between these two quantities so that 15 



The n c |AC(x)| vector 5 is defined as x* = [xi,--- ,x n J where xf = [x i;U --- , Xi-,\K(*)\]- 
n c x rc c |/C(x)| matrix S is assumed known and typically has the form 16 





0* • 


■ 0*' 




el • 




0* 


0* • 





where 0 is the if -vector of all zeros and & is the linear combination relating Xi to the x i]k - 
Under our assumptions for S, x t = In cases where the buyer wants to accumulate the 

results from suppliers (e.g., aggregating quantities) £ = 1 is the |/C(x)|-vector of all Is. In 
other cases the buyer may take £ = 1/|/C(x)| so that the time of delivery becomes the average 
time of delivery across the suppliers. Constraints on x become constraints on x via 

Gi(x)x = gi(x) and G 2 (x)x < g 2 (x) (7) 

where Gi(x) = Gi(x)S and Gi(x) = Gi(x)S. We might also expect the buyer to add 
additional linear constraints, such as requiring the latest shipment from any supplier to arrive 
earlier than a certain date, or requiring all deliveries to arrive the same day. There can also 
be constraints specific to particular suppliers, e.g., the buyer doesn't want any more than 100 
units from supplier 5. These can be handled simply as constraints on the individual Xi- k and 
added as extra rows to Gi(x), G 2 (x), gi(x), and g 2 (x). 

ls We can also relate x and x by x = Sx+i/ for some constant vector v, which will not cause any complications. 
However, there seems to be little practical reason to do so. 

16 With no additional computational complexity, we can allow S to depend on k. 
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With aggregation, the quadratic form to be minimized is (3x -fi(x)) C^ 1 (>f)(Sx-/i(^)) 
subject to the constraints given in Eq. (7). This minimization can be carried out through a 
straightforward generalization of the method given in Appendix D. 

6.2 Implementation 

In this section we outline an implementation of the entire e-procurement invention. We begin 
with a high-level description of the architecture, then fill in the details by describing a complete 
object model. 

6.2.1 High-level Architecture of the Invention 

There are at least two modes in which the invention may be used. First, the invention may 
reside at the site of large buyers, and suppliers who wish to sell to the buyer may be required 
to submit their capabilities via a web interface to the buyer. The invention may also be used 
within a marketplace hosted by a third party. Buyers/sellers log onto the market, submit their 
preference/capabilities, and act on the results. The architecture is modular enough to support 
both modes of operation. 

In Figure 1 we present an architecture for the invention. We describe the architecture, 
starting from the optimization algorithm which finds matches between buyers and sellers and 
work our way outwards. 

A controller surrounds the optimization engine, feeding it buyer preferences and seller 
capabilities. If multiple optimization processes are running (perhaps on different machines), 
the controller can also do load balancing, forwarding the request to the least busy process. 
The controller decomposes preferences and capabilities into their constituent buyer- and seller- 
specific versions (see below), selects the most specific matching preference/ capability pairs, and 
sends them to the matching engine for optimization. The controller then collects responses 
from the matching engine and returns them to the buyer. Additionally, the controller logs all 
results into a database for recording purposes. 

Another layer, called the Connector in Figure 1, separates the graphical user interface 
(GUI) through which users communicate with the tool from the controller. This layer serves a 
number of functions. The connector transforms the description of preferences and capabilities 
from the GUI into a form suitable for the implementation of the matching engine. Part of 
this transformation involves validation of appropriate input from the GUI layer so that no 
malformed input is ever sent to the controller. The Connector layer can also pull data from 
ERP or SCM systems and automatically infer preferences (using the total cost of ownership 
function) for the buyer. The enterprise abstraction layer insulates the Connector from the 
precise details of the manner in which the ERP and SCM data needs to be gathered. Total cost 
of ownership is evaluated in the simulation modules, which may either be running locally at the 
client's site or running centrally at the main server. These simulation modules pull operational 
data (the vector j0)from the enterprise abstraction layer. A preference optimization module 
(TCO) minimizes the total cost of ownership to determine the ideal trade and the flexibilities 
around the ideal trade. 
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At the outmost level, a layer provides integration with the GUI and/or host system. A 
number of administrative systems are expected at this layer. Market administration services 
allow easy definition of trading spaces, the dimensions of negotiation, limits on continuous 
variables, allowed settings of the discrete variables, etc. User administration services allow an 
administrator to define buyers, passwords, spending limits, etc. Supplier services accomplish 
similar tasks on the supply side. Managers for preferences, capabilities, and match results 
ensure that these objects are properly stored in a database. This layer layer also dynamically 
generates the html necessary for presentation of the data via a web interface to buyers and 
sellers. 

For maximal portability, communications between the View and Connector are via XML 
documents. For maximal efficiency, communications between the Connector and matching 
controller are as serialized Java objects. 

6.2.2 An Object Model for the Invention 

The fundamental objects required for the invention are preferences from buyers, capabilities 
from sellers, and match results returned to all parties. The components of such objects have 
already been considered from a mathematical point of view, and we now describe one possible 
computer representation. 

In this section we describe a complete grammar for the object model. The following 
syntactic conventions are used: 

• (nt) denotes a non-terminal symbol nt 

• [obj] denotes an optional grammar segment obj 

• {obj} denotes 1, or many times the grammar segment obj 

• -> denotes a production rule for non-terminal symbol. If there are multiple rules, say 
(a), (b), and (c), then these are denoted as 

(nt) -Ma>|(&)|<c). 

In contrast, a production rule of the form 

(nt) ^(a),(b),{c) 

indicates that the non-terminal (nt) is composed of three grammar segments, (a), (b), 
and (c) 

• terminal keywords are in serif font 

Obvious non-terminal grammar elements like (string) and (integer) are not described. 
Supply Side 

To represent capabilities that apply to a specific buyer (perhaps for contractual reasons), we 
have defined a capability to be a list of (buyerSpecificCapability). With one exception, a 
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buyer-specific capability applies only to one buyer - that buyer associated in the id field of 
the (buyerSpecificCapability). The exception occurs if the id field is * or wildcard. This 
indicates that the capability applies to all buyers. Using buyer-specific capabilities, sup- 
pliers can represent specific capabilities to certain buyers and generic capabilities applying 
to all other buyers. By not including a wildcard (buyerSpecificCapability) and only listing 
(buyerSpecificCapability) s applicable to specific buyers, sellers can also represent the fact 
that they will trade only with a subset of all buyers. In cases where both the wildcard 
(buyerSpecificCapability) and a (buyerSpecificCapability) applicable to a specific buyer ap- 
ply, the most specific (buyerSpecificCapability) is selected. 

A schematic of a (sellerSpecificPreference) is given in Figure 2. 

We begin at the top level of a capability: 

capability — >■ {(buyerSpecificCapability)} 

where 

(buyerSpecificCapability) ->• id: (id), 

discrete: {(discreteVarDescription)}, 
continuous: {(continuous VarDescription)}, 
range: {(range VarDescription)}, 
[d isc reteCo n st ra i nt : (discr eteConstraint )] , 
i nstance: { (discreteCapabilitylnstance) } 
[aggregation Participation: Oj 1] . 

(id) identifies a buyer or group of buyers. Individual buyers are represented by some unique 
identifier (say an integer) and the group of all buyers is identified by the wildcard character 
So we have 

(id) -> (integer) | * . 

aggregationParticipation is a Boolean flag giving the supplier's willingness to participate in 
aggregate orders to the identified buyer. 

Each of the variable constituent components is described by 

(discreteVarDescription) -4 name: (integer), 

a I lowed Values: {(integer)} 
(continuous VarDescription) name: (integer), 
min: (double), 
max: (double) 
(range VarDescription) -> name: (integer).. 



21 



In its simplest form, a (discreteConstraint) is a list of more primitive constraints 



(discreteConstraint) — > {(primitiveDiscreteConstraint)} 



where each primitive constraint is composed as follows: 



(primitiveDiscreteConstraint) — >• name: (string) 
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variables: {(discreteVarName)}, 

includes: 0 1 1, 

values: (integerMatrix) 



(discreteVarName) is the name of the discrete variable involved in the constraint 
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(discreteVarName) — V (integer). 



The includes field is a bit. If the bit is 1, then the combinations listed in the values field are 
the allowed values the variables may take on. If the bit is 0, then the combinations listed in 
values are the excluded combinations, i.e., everything in the powerset of the variables is allowed 
except those combinations listed in values. The order of the variable names is significant, since 
15 they will be assumed to be in the same order in values. If there are a variables involved in the 
constraint, and c constraints, then (integerMatrix) is an a x c matrix of integers: 

(integerMatrix) — ¥ (integer Vector), ■ • • , (integer Vector) 
(integer Vector) -> (integer) , • ■ ■ , (integer) 

20 The (discreteCapabilitylnstance) component is described by 



(discreteCapabilitylnstance) — > mask: (discrete Var Mask) , 
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[rangeCapability: { (range Varlnstance)}], 
continuousCapability: (continuousCapability) 
conti n uou sConst ra i nts: (continuousConstraints) 



A (range Varlnstance) defines a range variable and has the form 
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(range Var Instance) — *• name: (integer). 

min: (double), 
max: (double). 



The (discreteVarMask) relates to the discussion of 6.2.2. As in Table 3 we have 



(discreteVarMask) -4 {(extendedVar Value)} 
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where an (extendedVar Value) is either an integer from the domain of the discrete variable or 
the wildcard character '*': 

(extendedVar Value) — > (integer) | * . 

(continuousConstraints) describes the hard linear constraints for the continuous variables. 
5 Since these constraints may be either inequality or equality, we have 

(continuousConstraints) -» [equality: (linearConstraints)] , 
[i neq ua I ity : (linearConstraints)] 

Both the equality and inequality constraints are expressed through a matrix which is c x n c 
10 where c is the number of constraints, and a vector which is c x 1. Consequently we have 

(linearConstraints) — v matrix: (doubleMatrix), 
vector: (double Vector) 

A (doubleMatrix) is defined by 
15 (doubleMatrix) -> (double Vector) , ■ • • , (double Vector) 

and a (double Vector) is just what the name suggests - a vector of doubles: 
(double Vector) — >• (double), • • • , (double). 
The only remaining undescribed element above is (continuousCapability) whose description 

20 is 

(continuousCapability) -> breakpoints: (doubleListMatrix) , 

valAtBreakPoints: (doubleListMatrix) 

(doubleListMatrix) describes a n c x n c matrix whose elements are lists of (double) : 

25 (doubleListMatrix) -» (doubleList Vector) , • • • , (doubleList Vector) 

(doubleList Vector) — >■ (doubleList), • • • , (doubleList) 
(doubleList) ->{(double)} 

It is assumed that the rows and columns of the matrix are in some canonical order so that 
30 we know which continuous variable is referenced. A natural order is the one defined in 
{ (continuous VarDescription) } 

Preferences 

Just as capabilities may be buyer-specific so too may preferences be seller-specific. The same 
rules determining which seller-specific preference to apply are followed. A schematic of a 
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(sellerSpecificPreferen.ee) is given in Figure 3. 
We define a preference as follows 



(preference) -> {(sellerSpecificPreference)}[, (aggregatedPreference)] 

i.e., a preference is a list of (sellerSpecificPreference) with an optional aggregated preference. 
We first describe (sellerSpecificPreference) and then consider (aggregatedPreference). 
The (sellerSpecificPreference) is composed as follows 

(sellerSpecificPreference) — Hd: (id), 

d iscrete: { (discrete VarDescription) } , 
conti n uous: { (continuous VarDescription) } , 
range: {(range VarDescription)}, 
dimensionWeights: (dimension Weights), 
discreteTradeoff: (tradeoffTables) 
[discreteConstraint: (discreteConstraint)], 
instance: {(discretePreferencelnstance)} 

Of these elements, only (dimensionWeights), (tradeoffTables), and (discretePreferencelnstance) 
have yet to be defined. (dimensionWeights) gives the weights of all dimensions that indicate 
their importance. For convenience we break up the weights according to the three types of 
variables. Thus we have 

(dimensionWeights) range: (double Vector), 
discrete: (double Vector), 
continuous: (doubleVector) 

A (doubleVector) has been described previously. Each of the corresponding vectors is as long 
as the number of range, discrete, or continuous dimensions. (tradeoffTables) is an n discrete x 
^discrete matrix of (tradeoffTable): 

(tradeoffTables) -» (tradeoffTableMatrix) 
(tradeoffTableMatrix) -> (tradeoffTable Vector), ■ • ■ , (tradeoffTable Vector) 
(tradeoffTable Vector) -> (tradeoffTable), • ■ ■ , (tradeoffTable) 

A (tradeoffTable) is simply a matrix of double values. 

Finally, we turn to the last undefined component of a (preference). A (discretePreferencelnstance) 
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is composed as follows: 



(discretePreferencelnstance) -> mask: (mask), 
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[ra ngeldea I : { (rangeVarInstan.ee)] , 
continuousldeal: (double Vector) , 
tradeoffMatrix: (doubleMatrix) , 
[continuousConstraints: (continuousConstraints)] 



The ran gel deal and continuousldeal fields give the desired range and continuous trade pa- 
rameters. The tradeoffMatrix field gives the positive definite matrix expressing the tradeoffs 
10 amongst the continuous variables. (continuousConstraints) have been described previously in 
the sell-side specification. 

To complete the specification of preferences, we conclude with the definition of (aggregatedPreference) 
Refer to the discussion of section 6.1.7 for details. 



25 In the above definition, the previously defined elements maintain their meaning. The addi- 
tionalConstraints field allows the buyer to express constraints around the aggregation, such 
as "all orders must arrive on the same day," etc. participants lists the suppliers who can 
participate in the aggregation and their characteristics. Note that if the wildcard supplier 
participates, the order can potentially be aggregated across all suppliers. (aggSpecification) 

30 describes information specific to a supplier participating in the aggregation. It is defined by 



id identifies the participating supplier and constraints specific to that supplier defined in an 
accompanying (sellerSpecificPreference) will be used in the optimization. Additional informa- 
35 tion may be added as required. The contributionType field is used to define the f vectors used 
in aggregation. The (contributionType Vector) consists of n c elements indicating the type of 



(aggregatedPreference) -» participants: {(aggSpecification)}, 
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contributionType: (contributionType Vector), 
additionalConstraints: (continuousConstraints), 
discrete: { (discreteVarDescription) } , 
continuous: { (continuous VarDescription)}, 
ra nge: { (range VarDescription) } , 
dimension Weights: (dimension Weights), 
discreteTradeoff: (tradeoflTables) 
[discreteConstraint: (discreteConstraint)] , 
i nsta n ce : { (discretePreferencelnst ance ) } 
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(aggSpecification) id: (id). 
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discrete mask 


output 


specificity 


[. ♦ ,j 


{continuousl, rangel} 


0 


U * ,] 


{continuous2, range2} 


1 


[>Cl * x 3 J 


{continuous3, range3} 


2 



Table 3: Example of discrete masks for specifying continuous and range variables which are 
dependent on discrete variables. >c\ and x 3 signify specific values for the first and third discrete 
variables. The specificity of each mask is indicated in the third column. 

aggregation for each continuous dimension: 

(contributionTypeVector) — > (contributionType), ■ • ■ , (contributionType) . 

Possible contribution types include 

5 (contributionType) — >■ sum, average, zero. 

sum sets £ = 1, average sets £ = l/\K-(x)\, and zero sets £ = 0. 
Masking 

We have allowed constraints, ideal values, tradeoffs, and continuous capabilities to be depen- 
10 dent on discrete variables. In this section we describe an efficient manner in which to encode 
this dependence. 

The data structure must represent continuous and range variables for all valid discrete 
inputs. An efficient way to do this is to use hierarchical definitions. At the top of the hierarchy 
are the definitions of the continuous and range variables for the discrete values *£* = [*,■■• , *]. 

15 These values apply to all v unless more specialized masks are defined. A more specialized 
mask of the continuous and range variables is specified by defining values for some of the 
components Xj. The more components that are defined, the more specialized the definition. 
The most specific mask is always used. An example definition for three discrete variables is 
given in Table 3. The response to the lookup [xi x 2 x 3 ] is {continuous3, range3} if and only if 

20 *c\ = x\ AX3 — x 3 , {continuous2, range2} if and only if x\ = yt\ Ax 3 ^ X3, and {continuousl, 
rangel} otherwise. 

Match Results 

Returned to the buyer is a list of matches with different suppliers, which can be ranked and 
viewed in many different ways in the GUI. A (matchResultList) is a list of matchResult: 

25 (matchResultList) -> {(matchResult)}. 

A match result may either be a (singleSupplierMatchResult) or an (aggregatedMatchResult): 

(matchResult) — > (singleSupplierMatchResult) | (aggregatedMatchResult) . 

30 A (singleSupplierMatchResult) represents the best trade with a single supplier and is composed 
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of the following elements: 



(singleSupplierMatchResult) — > supplierld: (integer), 
utility: (double), 
feasible: 0|1, 
[costFactors: {(double)], 
continuous: {(double)}, 
discrete: { (discrete VarDescription)}, 
range: (range Varlnstance) . 

The supplierld indicates the supplier sourcing this trade and the utility field indicates the utility 
of the trade (which can be used to rank the trades), feasible is a bit indicating whether or not 
a feasible trade with this supplier was found. The continuous, discrete, and range fields list the 
respective trade parameters determined by the matching algorithm. The optional cost factors 
field lists the constituent costs contributing to the total cost of ownership C 0 evaluated at the 
trade point returned in the (singleSupplierMatchResult). 

An (aggregatedMatchResult) returns the optimal trade when the buyer has requested 
aggregation. It is composed of the following elements: 

(aggregatedMatchResult) — > utility: (double), 
feasible: 0|1, 
[cost Fa ctors : { ( double) ] , 

supplierTradeParameters: {(supplierTradeParameters)}. 

As before, the utility field gives the utility of the aggregate trade, and the feasibility flag indicates 
whether or not a feasible aggregate trade was found (there may not be if the constraints 
were too stringent). costFactors can also be returned in C 0 is sufficiently general to handle 
aggregated trades. Finally, (supplierTradeParameters) lists the trade parameters for each 
supplier involved in the aggregation. It is defined as follows: 

(supplierTradeParameters) — ¥ supplierld {integer), 

continuous: {(double)}, 

discrete: { (discrete VarDescr iption) }, 

range: (range Varlnstance) . 



6.3 Summary 

We have described an efficient computational procedure in which to encode buyer's trading 
35 preferences and hard constraints, supplier's delivery capabilities and constraints, and optimize 
to find those matches between one buyer and one or many sellers that maximize the buyer's 
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utility. By optimizing against both qualitative and quantitative factors, and exploiting the 
trading flexibilities possessed by both parties, the system determines better trades. The tool is 
particularly useful as companies move their direct material purchasing online. By optimizing 
across flexibilities, win- win trades are discovered for both trading parties. 
5 The representation of trading preferences is designed to be expressive yet easily elicitable 

from a buyer, and computationally tractable. The representation of supplier capabilities was 
chosen to parallel the manner in which suppliers already think of their delivery capabilities and 
seamlessly includes volume discounts and incremental costs. These supplier capabilities may 
be part of an online catalog. The representation of the negotiation space is rich, supporting 

10 three types of variables. 

We have outlined a manner in which preferences may be inferred automatically through 
models of the purchasing company. Such models incorporate many cost factors, taking the 
total cost of ownership into account. The system provides trades which minimize the total 
cost and represent significant new savings. 

15 The invention can operate both at a buyer's site, where suppliers input their capabilities 

through an HTML interface to the world wide web or as an embedded part of an electronic 
market hosted by a particular web site. The invention may operate at regularly scheduled 
intervals or sporadically in lieu of current request for quotations (RFQ). The buyer may 
broadcast a RFQ event to suppliers, indicating a time within which suppliers must respond. 

20 At the close of the event, the buyer can use the present invention to assist in the analysis of 
the supplier responses. 

Complex algorithms have been specified which should permit most matching optimization 
to occur in near real-time. The rapidity of optimization, combined with graphical what-if 
tools, allows for analysis and exploration of trades, which should significantly improve the 

25 quality of purchasing decisions. 

6.4 An e-Commerce Infrastructure for Value Chains 

In this section we describe in detail how the proposed infrastructure delivers on the promises 
made in the Summary of the Invention. We begin by describing major innovations in the 
present invention and how they are all used synergistically. 

30 6.4.1 Major Innovations 

The most broad invention combines at least four advances: 

1. multidimensional automated markets (hereafter simply markets) which capture many 
aspects of value. 

2. algorithms and interfaces which implicitly allow for consumers to express the preferences 
35 over the multiple dimensions of value 

3. linked markets allowing for complex assembly of products 
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4. specification of the constraints (both logical and numerical) inherent between markets 
to allow for coordinated buys and sells between markets 



In addition, inventions described in a patent application titled, "An Adaptive and Reliable 
System and Method for Operations Management", Application Number 09/345,441 filed July 
5 1, 1999 (the contents of which are herein incorporated by reference)can also be used in 
conjunction with the present invention. These other inventions are: 

5. the use of models and optimization algorithms to optimally determine the best bids to 
submit to the automated market 

6. the use of subset relations is-a, has-a etc. to automatically construct the constraints 
10 between markets 

6.4.2 Integrated Picture 

The internet and e-commerce are changing the way consumers and businesses trade with each 
other. One interesting trend is the development of centralized economic hubs {e.g. eSteel, 
ChemDex) through which all e-commerce in a particular domain flows. This trend is expected 
15 to continue and become more prevalent. The invention described here is the tool that drives 
these economic hubs. 

Before describing the components and innovations in detail we describe the overall architecture 
of the integrated system. For concreteness we focus on the invention as applied to the personal 
computer (PC) supply chain. We note that the ideas can be applied equally well to supply 
20 chains in other industries. 

The infrastructure is composed of three major components all running on different computers. 
A central server (or more likely servers) coupled to the economic hub serves as the source 
of consumer demand. Other sets of servers act as the trading markets themselves relaying 
information between buyers and sellers. Finally, running at each market participant's site is 
25 client software coupled to the relevant markets. 

Consumer Demand: Consumer demand is pulled through the supply chain through a set of 
coupled reverse auction 17 markets. An interface iteratively elicits multi-dimensional consumer 
preferences which is then forwarded to a top level market where suppliers offer to fulfill the 
consumer request. The functioning of the consumer interface is described in more detail in 
30 section 6.4.3. 

Markets: Markets within the supply chain are represented by servers which link trading 

partners. The market broadcasts demand, collects responses, and forwards results back to the 

source of the demand. See section 6.4.4 for more information on the functioning of markets, 

particularly the coupling between markets. 

35 Client Companies: Each trading participant (whether buying or selling) is represented 

by client software connected to the appropriate markets. The client software both initiates 

17 A reverse auction is like an auction except that a request is broadcast to suppliers who bid on fulfilling the 
request. Unlike traditional auctions the proposed market clears in one shot, i.e. each supplier has only one 
chance to bid. 
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requests and responds to requests from other market participants. The client software running 
at the company also maintains a memory of all past transactions and can eliminate those that 
are out of date. A user interface allows companies to define the operation of their interface to 
the markets. See 6.4.5 for further details. 

5 Supplier Push As has been mentioned the infrastructure is organized around demand 
pulled through the supply chain by the consumer. Interestingly, the same pull infrastructure 
can also be used by suppliers to push demand. In addition to requesting and responding to 
requests, the market interface running locally at company sites could also be used to advertise 
capability (capability not in response to any particular demand) to the market which will then 
10 forward the advertised capability to any participant connected to the market. 

6.4.3 Climbing in Preference Space: the consumer interface 

The consumer drives the integrated infrastructure by generating demand. In our example, 
a graphical user interface (GUI) at the PC economic hub provides a centralized interface 
through which all computers are purchased. The GUI provides a number of advantages to 

15 both consumers and computer vendors. A computer is described across multiple dimensions of 
value. Some of these dimensions include: price, memory size, processor speed, hard drive side, 
removable storage, monitor size and resolution, financing, warranty, delivery date, etc. For 
each dimension the user can identify a preference. For example the consumer may wish 128Mb 
of memory, a 600MHz Pentium III processor with a 15Gb hard drive for $1500 delivered within 

20 3 weeks. Additionally, it may be important for the consumer to express constraints that must 
be satisfied. To keep things as simple as possible we express an optional inequality constraint 
for each dimension, e.g. price must be less than $1600 and the hard drive must be at least 
12Gb. 

Optional Dimensions: This example brings up an interesting feature of the dimensions of 
25 value. Imagine for a moment that one of the dimensions was portability with possible values 
being high portability, medium portability, and low portability. In each of these cases we might 
consider lightweight laptops, heavy laptops, and desktops. Any one of these choices might 
then invoke additional situational dimensions. For example, if we select high portability then 
battery life (in hours) becomes an important dimension which is only applicable to laptops. 
30 There are a number of ways to treat these optional dimensions. The simplest possibility is to 
have separate markets for laptops and desktops with appropriate dimensions for each. A better 
solution is to have a GUI which only turns on the optional dimensions when appropriate. If 
low portability is selected by the user the battery life dimension remains lightly greyed out 
and inaccessible. If medium or high portability is required then the battery life dimension is 
35 activated. In the technical details section we show how this can easily be accomplished. 

Dimension Weighting: Optionally, the user may also specify an importance to each 
dimension. If hard drive size and price are paramount these two dimensions have high value 
while all other dimensions have low value. The weighting is used to identify similarity between 
computers; two computers which are quite different only in one dimension which is not highly 
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weighted by the consumer are considered quite similar. See the technical details section. 

Threshold Constraints: A final way in which consumers may specify their needs is through 
constraints expressed for each dimension. The constraints are simple to state. The consumer 
specifies a threshold and whether or not he wants to be above or below that threshold. For 
5 example the consumer may specify that the monitor size has to be at least 17 inches and he 
will not pay any more than $1500. Knowledgeable consumers may define all these settings 
directly or the software may infer settings for less experienced consumers based upon simple 
questions (like the major uses of the computer). 

Branding: Consumers may also have brand preferences for a computer or any of its 

10 components. The proposed GUI allows for users to express these preferences. There are 
two alternative ways in which branding might work and depending on the application either 
may be appropriate. In the first alternative brand name is a hard constraint and no computers 
are shown to the consumer unless the brand name matches the consumers desire. In the second 
alternative brand name is a soft constraint and the consumer merely prefers one brand over 

15 another. 

In either case the GUI allows the user to specify for each dimension whether or not brand 
matters and if so which brand/brands is/are preferred. 1 8 If brand is a hard constraint this 
constraint is propagated to the top level market and all suppliers of the top level market so 
that no responses not respecting the constraint are ever returned to the consumer. If brand is 

20 a soft constraint it is incorporated into the distance function dj as described in the technical 
details section. Otherwise identical components which match the desired brand are closer than 
those components which do not match the desired brand. 

Iterative Hillclimbing and Recombination: Once the consumer has specified a desired computer 
defined as a point (and perhaps weights and constraints) in the high dimensional value space, 

25 a query is sent to a top-level automated market. Sellers of computers (e.g. Micron, Gateway, 
Dell) attached to the market respond to the consumer request by returning offers (one or 
many) to the market for candidate computers again represented in the high-dimensional space 
of value. Responses to the consumer are filtered to satisfy the specified constraints. The GUI 
allows the consumer to sort responses based on any of the dimensions (see the dj function 

30 in the technical details section) or according to a global measure of distance between two 
computers (see the d function in the technical details section). 

The consumer may modify any of the responses (a mutation) and use this as a resubmission 
to the market. This method of exploring the space of possible computers according to the 
consumer's desires and the supplier's availability can broadly be summarized as climbing in 

35 preference space. The climbing can begin from the most recently visited computer configuration 

or from any configuration that has been saved in a history list. 19 It is also perfectly feasible 

to allow the consumer to recombine desirable computer configurations. For example the user 

may ask for configurations which are "in between" two previously examined configurations 

and combine the advantages of both configurations. We also might allow the consumer to 

18 We might also allow consumers to rank their brand preferences but this extension seems unnecessary in 
light of the iterative nature of product exploration. 

1 History lists extend both within and between sessions. 
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submit more than a single configuration to the top level market and receive responses around 
all submissions. 

Benefits of the Approach: This procedure offers a significant benefit over other approaches 
to multidimensional markets. In systems like OptiMark 2 0 the consumer (in this case a financial 
5 trader) must a priori specify preferences over all dimensions (price and quantity). The explicit 
specification of preference over even two dimensions is a difficult task and it is unreasonable to 
expect consumers to be able to do this. The present innovation allows consumers to iteratively 
explore preference space thereby implicitly defining preference in a friendly manner. 

This method of interactive consumer exploration is suitable for any complex product and 
10 we expect the application to be broad. The present invention has many other applications 2 1 

Technical Details We have mentioned that computer configurations can be sorted on all 
dimensions according to how closely they match the consumers preferences. In this section 
we describe how this sorting might work and how it interacts with the weights and brand 
preferences attached to each dimension. 

15 The value space defines a notion of distance or similarity between products (in this case 

computers) defined as follows: If there are a total of D dimensions 2 2 of value then a computer c 
is described as a D- vector of the settings for each dimension. Let i label the components of c so 
that for example c\ = memory is c 2 = processor speed etc. For example, the memory dimension 
might assume any of the values ci = 64Mb, eg = 128Mb, or ci = 256Mb and similarly for all 

20 other dimensions. Let aj be a binary variable (i.e. aj takes the value 0 or 1) indicating which 
dimensions are actually used. If c, is the dimension corresponding to battery life in hours 
then a,- = 1 indicates this dimension is appropriate while a,- = 0 indicates the dimension is not 
appropriate. The vector whose components are the a, is represented by a. 

A consumer preference also includes optional weights and preferred brand. The weight 

25 (or importance) of dimension i is indicated by Wj. All Wj are non-negative and normalized 
so that Y^?_i ajWj — 1. The vector of all weights is represented by w. The brand preferences 
for dimension % are expressed as a (possibly empty) set Bj which includes desired brands. 
B is the set of all Bj. A complete consumer description is thus represented by the vector 
C = {c, a, w, B} with components Cj = \ci,ai,Wj,Bi}. 

30 The distance d(c, c') between computers c and c' is then defined as 

d(C;C , ) = EL^( C ,c;) (8) 

where dj(Cj,C') measures the distance along the single dimension i. To allow for standard 
comparison between computer configurations we normalize all dj so that they lie between 0 
and 1. 

20 U.S. Patent 5,845,266, "Crossing network utilizing satisfaction de nsity profile with price discovery features" 

dated 12/1/98. 

e.g. credit cards where the dimensions might be interest rate, yearly fee, grace period, benefits (like frequent 
flyer miles, etc.) . 

22 The D dimensions include both required and optional dimensions. 
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We recall from previous discussion that for soft brand constraints the matching of a brand 
to its desired value is incorporated into the distance function. Thus we generally write the 
distance function dj as the sum of two terms, one for brand and one for similarity. If bj is an 
indicator variable whose value is 1 if the consumer has identified that brand matters for the 
5 ith dimension (i.e. Bj is non-empty) and 0 otherwise, then we write 

di{CiA) = ogbidfXBuB'i) + (1 - aAK^(c»,cQ. (9) 

The cti parameter determines the relative importance between brand name and similarity- 
Lacking better information we might take all ctj = 1/2. 2 3 

For simplicity we assume that the brand distance function, d^ (cj, c'), is 0 if there is overlap 
10 in the preferred brands for a and c' and 1 otherwise, i.e. 

d\\Bi,B^ = l . (10) 

0 otherwise 



The similarity distance measure, d)^ is only slightly more complicated. The allowed 
values for most dimensions are discrete over a finite range. If the ith. dimension is bounded 
and ranges from cf in to cf 1 ^ then we might define 

though other choices are certainly possible and may be more appropriate depending on the 
nature of the dimension. In many situations it is appropriate that the distance function is not 
symmetric, i.e. d\ s) {a,c' 4 ) ^ d^jc'^a). Examples of asymmetry are easy to find. A consumer 
may be more than willing to accept an extra few Mb of memory but will be very dissatisfied 
receiving a few Mb less. As a concrete example we may offer open ended ranges on most 
choices. The range is defined by a single point x{ and an indication of whether the range is 
above or below that point. In the case of the range \x[, oo) we can take the distance function 
measuring the distance from the range to be d(xj, x';) = {x^ - Xi)0{x' { - Xj). In the opposite 
case where the range is \0,x'A the distance function can be d{xj,x'^ = (xj - x'^ejxj - x' f ). 
Accordingly, the present invention also includes possibilities for d^ {cj,d 4 ). 

Distance Ranking: The distance metric of Eq. (8) allows for a global ranking (after 
configurations not satisfying the constraints have been filtered out) of vendor responses based 
upon their similarity to the original consumer request (smaller distances are more similar). 
If C req describes the consumer specified computer and a labels responses from vendors in 
response to the request then the automated market sorts the responses based upon d{C ma , 
and returns a specified number back to the consumer interface. The interface also allows 
sorting (in order of decreasing distance) on each dimension according to the dj. The consumer 
then iterates the process starting from any of the returned offers, modifying them slightly 

23 Better values for eg might be identified from real consumer interaction with the GUI. 
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and re-submitting to the market. In this way the consumer can move about in the space of 
preferences climbing t owards ever more desirable computer configurations. When the consumer 
finally does identify a c onfiguration that he desires to purchase it is a simple matter of hitting 
a buy button. 

6 - 4 - 4 Integrating the Supply Chain: cascading markets 

It is natural to consider the extension of high-dimensional automated markets to other transactions 
made in the supply chain. In fact, this may be essential to fulfill the kind of real-time 
consumer preference exploration described in section 6.4.3. These other markets will also 
typically be high dimensional but the dimensions will differ from the top-level market. For 
example there may be a hard drive market where buyers of hard drives (whether computer 
vendors, consumers, or v alue added retailers) meet sellers (e.g. suppliers like Seagate, Western 
Digital). The dimensio ns of the hard drive market might include: price, size (in Gb), type 
(IDE vs SCSI), quantity, delivery date, etc. 

Operation of Cascadi ng Markets: The basic idea is straightforward. Before a computer 
vendor responds to con sumer requests in the top-level market the vendor may send its own 
requests to submarkets to purchase products and services it needs to meet the consumer 
demand. Suppliers t o the computer vendors may in turn go to other submarkets before 
responding to compute r vendor requests. A single consumer request may then result in 
real-time cascading inquiries to a tree of nested markets all the way down to the beginning of 
the supply chain. We c all the market which originated a request the parent market and the 
submarket a child market. It is a requirement that any possible path through the cascading 
markets is acyclic. 

When the consumer moves on to examine a new configuration all relevant market participants 
{i.e. those that respo nded to the initial customer request or any concommittant subsequent 
request) are informed and their offers are released. If, on the other hand, the consumer elects 
to purchase the configura tion all participants are informed and are held to their offers. When 
combined with emerging XML standards for e-cornmerce (e.g. RosettaNet in the IT industry) 
the necessary "paperwork " for confirmed purchases can be generated by the computer and 
sent betwe en all relevant parties (including the consumer after a credit card check). 

Note that it remain s unlikely that every consumer purchase will trigger a cascade of 
transactions throughout t he tree of markets because it is not profitable in most supply chains to 
order and transport in si ngle unit batches. Thus, only when inventories need to be replenished 
will companies in the supply chain go to the automated market in response to a request. 

Auxiliary Information: It is useful to pass additional information between markets beyond 
merely the parameters desc ribing the trade. There are many reasons extra information may 
be desired and we describ e two applications of auxiliary information. Most simply, in cases 
where it is important to identify who the buyer or seller actually are an additional identifier 
may be passed labeling each trader. 

A more serious applic ation addresses a potential problem of nested markets. A problem 
with cascading markets is that demand can spuriously be amplified as it moves between 
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markets and away f rom its originating source. Consider a consumer demand for 100 computers. 
If two of the com puter suppliers are both low on hard drives they might both submit a 
request to the hard drive market for an additional 100 drives. If the hard drive company 
were to estimate dem and based solely upon the hard drive market then in this case the 
compan y would falsely assume that the demand was twice as high as it really is. This 
problem is easily fix ed by supplying auxiliary information along with trade parameters. A 
request that comes i nto a toplevel market is assigned an identifier indicating the market and a 
unique identifier. 24 Every request to submarkets that is in response to a request in the parent 
market i s again assigned both market and trade identifiers. These identifiers are appended 
to the identifiers from the parent market when the request is broadcast to all suppliers 
listening to the marke t. Returning to the hard drive example, if the top-level market has 
market identifier ASSEMBLE and the hard drive market has identifier HARDDRIVE and if 
the original consumer request had trade identifier 4561133 and the market trade identifiers 
were 56891 and 56892 respectively then the identifiers which are passed to all hard drive 
suppliers in response to both hard drive demands are ASSEMBLE:4561133,HARDDRIVE:56891 
and ASSEMBLE:4561133,HARDDRIVE:56892 respectively. In this way all hard drive suppliers 
can see that the request originated with the same demand. In fact such mechanisms might 
also allow for improve d forecasting of demand since all demand is accompanied by its path 
through the entire supply chain which may contain additional useful information. 

Expiration Dates: One of the advantages of the present e-commerce infrastructure is 
the real-time setting of trade parameters. Companies can use the system to optimize their 
operations continuousl y. To achieve the maximal benefit from continuous optimization it is 
important that the o ffers made by any market participant have a bounded lifetime. What 
is optimal now may not be in a day (or even an hour). Thus an additional application of 
auxiliary information is the inclusion of expiration dates after which the trader parameters 
are invalid. Unless pe rishable trades are executed before the expiration date they are removed 
from the system. 

Other B2B Market Issues The cascading submarkets are all business to business markets 
and as such can operate under different rules than the business to consumer market. In this 
section we address what some of these differences might be. 

Clearing Mechanis ms: We have proposed that the B2B markets clear in a fashion similar 
to the B2C market. One important difference is that these markets probably must clear in 
one shot in order to d eliver real-time feedback to the end consumer. There are certainly 
other ways in which the B 2B markets clear. One possibility is to use optimization tools to 
define satisfaction surfac es just as in OptiMark. In this way the parameters of the trade are 
determined by softwar e to optimize some .joint criteria. Ideally, we would like protection for 
other B2B market clearing mechanisms. 

Aggregation: Anoth er potentially important function of automated markets is the aggregation 
of orders to meet demand. Forlargeorders it may be the case that no single supplier can meet 
24 All trades in that market will have distinct identifiers. 
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the requested demand. Market software can automatically collect responses to fulfill orders. 
This software, in addition to coordinating between purchases in different markets (see section 
6.4.5) can also be used to coordinate multiple orders within the same market. 

Many-to-many markets: Thus far the market mechanisms we have described all result in 
5 one-to-one trades (except if we aggregate orders in which case a trade may be many-to-one). 
For future reference we also note that the markets themselves may be more complex and used 
for many-to-many trades. 

6.4.5 Coordination Between Markets: inter-market language 

The assembly of a product or service from constituent products or services often requires 

10 coordinated purchases from a number of supply markets. 

Logical Constraints: As illustration, there is no point in ordering a part if transportation 
can not also be arranged at the appropriate time. Both events must occur or else neither of 
the purchases should be made. This is an example of a Boolean logical constraint (in this case 
the logical and function) that exists between multiple markets. Of course in more complicated 

15 situations there may be more complex logical constraints between markets. 

For added expressivity we propose using linear logic as the language in which to express 
logical constraints. Linear logic offers compelling advantages over Boolean logic since it 
subsumes Boolean logic and is resource sensitive exactly as real supply chain transactions 
are. Efficient implementations of linear logic exist so that checking of linear logical constraints 

20 is straightforward. The present invention includes support for both Boolean and linear logical 
constraints. 

Numerical Constraints: There is another important class of relationships that may exist 
between markets. In most complex processes certain events have to occur sequentially; the 
process can not progress until some condition is fulfilled. Our shipping example also provides 

25 an example of this. There is no point arranging the purchase of transportation if the only 
available pick up date is before the item to be shipped can be completed. This is an example of 
a numerical constraint that exists between the parameters of transactions on different markets 
and optionally some description of the state of the company. Other examples might include 
cost constraints which require that the total expenditures on all purchases to be less than 

30 some threshold. 

A Constraint Language: Any automated market solution to supply chain coordination 
problems will require these constraint mechanisms and generically these constraints will be 
either logical (as in the and constraint) or numerical (as in the cost constraint). We propose an 
Inter-Market Language (IML) in which to express these constraints so that the nested market 
35 solution always respects the user specified constraints. 

IML constraints are specified as follows. Let mj be a logical variable which is true if a 

transaction is made on the ith. market (and false otherwise) and let mi{x.i) be a logical variable 

which is true if the transaction occurs with multidimensional parameters x^. 2 5 Furthermore let 

25 m.j and m,(x t ) are related by rm = V x ,to,(x) where by V x , we mean the logical or of m 1 (xj) over all 
possible values of the parameters Xj . 
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s be a vector of numerical values (for example describing current inventory levels, lead times, 
etc) and ~s be a vector of logical values describing the state of the company. The two classes 
of IML constraints which both must simultaneously be satisfied are then expressed as 

%i,-,m w r«) (12) 



iV=(xi,-.- ,x M ,s) = 0, JV<( Xl) --- ,XAf,s) < 0 (13) 

where B(-) is an arbitrary logical function, there are M relevant submarkets, and the functions 
iV=(-) = 0 and JV<(-) < 0 express the equality and inequality constraints respectively. 

Operational Use of Constraints: Operationally the coordination amongst different market 
transactions is achieved as follows. Software representing the interests of a vendor (and 
therefore probably running at the vendors site) expressed in IML listens to servers representing 
various markets. When a market request that the vendor wishes to respond to is heard the 
software examines the companies internal state {s, §} and determines what items, if any, need 
to be purchased on child submarkets. The software defines the purchase by defining a point 
or the desired parameters for each trade in each market. These requests are sent to servers 
representing each child market. The client software running at each company then waits for 
responses from each child market. If there are M total markets queried and rij responses from 
the ith. market then there are a total of fjf-i n>i possible combinations of deals. The software 
then examines each deal filtering out those combinations which fail to satisfy the constraints 
of Eqs. (12) and (13). Thus we refer to this aspect of the client software as filtering. Those 
combinations which pass the constraints can then be ranked (either by a human user or by 
additional company specific criteria (see reliability below and in section 6.4.7) added to the 
client software running at the vendor site) and sent back as responses to the original request. 

We have operationally described how constraints and filters may be used to coordinate 
purchases from automated markets. The converse problem is the coordination of sell orders 
to automated markets. This is less likely to be problematic since usually there is only one 
sell order being responded to at any given time and thus there are no sell-side coordination 
problems. Nevertheless, if there are sell-side coordinating activities that must occur the same 
mechanism can be used. A separate set of logical and numerical constraints can be defined in 
IML for sell orders. Any responses to requests can then be checked against these constraints 
to ensure no responses are sent which violate the selling practices expressed through the 
IML constraints. Sell-side filters can be used to enforce preexisting contracts like specially 
arranged pricing for preferred customers. This may be particularly important when responses 
to requests are submitted automatically by software coupled to the companies ERP software. 
Typically such software will appropriately construct a response but when no human is present 
to check responses sell-side filters provide a final sanity check. 

Other Applications of Constraints: An important consideration for any company in a 
supply chain is reliability. The company must be assured that its suppliers meet their 
obligations. In current supply chains this consideration is important enough that companies 
often limit themselves to a handful of trading partners with whom they maintain long term 



relationships. In order to fully reap the benefits of automated markets where many vendors 
may fulfill a need we need to resolve the issue of reliability. Filters are one way to address 
this issues. If for some historical reason a particular company doesn't like some suppliers the 
company may add these to be filtered. Thus the internal state f s, s) may include references to 
5 unwanted suppliers. Alternatively, reliability may be used to rank bids that pass all filters. 

6.4.6 Practical Concerns 

In this section we show how real world concerns of scalability, real-time performance, reliability, 
and security can be assured within the proposed e-commerce framework- 
Scalability The utility of an automated market is directly dependent on the liquidity 

10 (volume of trades) of the market. The greater the number of participants the greater the 
value to any individual participant. It is thus essential that the system scale gracefully to 
large numbers of market participants. Because the proposed e-commerce infrastructure is 
distributed scaling to large sizes is not a problem. 

The market is driven by consumer demand pulled through the supply chain. Interfaces to 

15 consumer interaction can run locally on the consumer's computer as Java language applets. 
The central server (or servers) for consumer interaction is then only responsible for relaying 
consumer information (expressed in compact XML documents) into the top level market and 
returning the list of consumer responses to the appropriate user computer. The workload 
is shared across multiple machines with more machines involved the greater the number of 

20 consumers. 

Similarly, the automated markets themselves serve mainly as relays of information and 
have little computational demands placed upon them. The markets send out requests, collect 
responses, and inform winning bidders of actual purchases. The bulk of the necessary computation 
is in the software which determines the parameters of the bids and responses. Since this 
25 software runs locally on each company's computer the system naturally exploits the parallelism 
inherent with increased vendor participation. 

Real-Time Response To be most useful to consumers it is important that consumer 
feedback is in real or near real time. This may be problematic on the rare occasions in which 
a query triggers a cascade down through submarkets. Before the top level market can clear 
30 {i.e. for responses to come back to the consumer) the tier one market must clear and before 
the tier one market can clear the tier two market must clear, etc. Fortunately, it is likely that 
all B2B markets will be monitored by software so that turnarounds on all markets will be fast. 
Market time-outs beyond which responses will not be accepted ensure timely responses from 
all submarkets. 

35 Nevertheless, in situations where the response time is slow we can offer the consumer an 

immediate but approximate response and warn the consumer that the response may not be 
accurate. This can be accomplished by caching previous responses and estimating the response 
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by interpolation of similar previous responses. Quick algorithms for the interpolation are easy 
to devise as long as the cached data can be accessed quickly as in any modern database. 



6.4.7 Reliability 

For any company to rely on the proposed infrastructure to run its supply chain it must 
5 be assured that the system is always up and running. For the system to be functioning 
for a company its own local server has to be up as well and all markets must be running. 
Each company is responsible for ensuring that its local servers are operating. The markets 
themselves will be served by a network of computers sharing the responsibility for managing 
transactions. The software which manages transactions can ensure that if a connection is 
10 broken during mid-transaction (e.g. if a computer goes down) it will be correctly reconstituted 
and resent. Such software is currently available from software companies. 

6.4.8 Security 

Whenever real dollars are being exchanged electronically there will always be important issues 
around security. We address each of these issues. 

15 Commitment Assurance: It is important the every participant is confident that contra 

parties honor their commitments whether in terms of payment or delivery. We have described 
one mechanism that allows companies to automatically filter out those companies with whom 
they have had bad past experiences. It might also be expected that word of any party 
defaulting on its commitments would rapidly spread to all other participants making any 

20 default very costly for the offender. 2 6 If these mechanisms are insufficient the infrastructure 
may require each participant to pay an up front fee from which offenders compensate their 
contra partners. 

Authentification: The system must ensure that all parties are who they say they are. This 
is an old problem faced by internet participants and has been solved satisfactorily by systems 
25 like PGP etc. 

Privacy: It is also important that companies private data remains private. If a competitor 
could intercept a companies responses to the market it would have a distinct and unfair 
advantage. All communications will be strongly encrypted to assure that this can not occur. 

Information Visibility: Similarly we may want different companies to have differing access 
30 to information. 27 For example data about consumer demand may be obtainable by a participant 
but only at a fee. This information may be served based on the authentification of the 
participant. 

6.4.9 Optimization of the Supply Chain 

The final ingredient in the infrastructure has been alluded to but not explicitly discussed. To 

35 achieve maximal benefit from the flexibility the proposed invention introduces at all levels of 

26 We might post defaults publicly so that other parties rapidly become aware of this knowledge and add the 
offender to their filters. 

^ Y This will certainly be the case if we wish to charge for access to some of this information. 
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the supply chain, real-time optimization systems can be used. This innovation is described in 
a patent application, titled, "An Adaptive and Reliable System and Method for Operations 
Management", Application Number 09/345,441 filed July 1, 1999 (the contents of which are 
herein incorporated by reference). 
5 Coupling to ERP Software: The optimization, systems monitor (in real-time) the state of 

the company and extrapolate from models of the company in its supply chain the optimal 
responses to requests that come in on the automated markets that the company supplies to. 
To achieve the best benefit, the servers or filters running at each company are coupled to ERP 
software to determine the best real-time response. The optimal response will be a function 

10 of the current instantaneous state of the company (as recorded in the ERP software) coupled 
with modeling and planning capabilities. A model of the likely future short term evolution 
of the state of the company may be used to best exploit the flexibility inherent in meeting 
the demand from the automated market. Such models might also take into account future 
expected demand and future responses from suppliers. 

15 Another important optimization that might yield significant payoffs is in the tradeoff 

between making a response that is very close to meeting a requester's demand versus the 
savings obtained by perhaps being further from the requesters true desires. 

6.4.10 Self- Organizing Supply Chains 

In the procedure outlined above it is the responsibility of each company to define its business 

20 rules as expressed though its operating constraints and ranking criteria. This is not a unreasonable 
request and is a formalization of the company's current operating procedures. However, the 
present system does suggest a novel extension for the supply chain of the future. We can view 
the activities along a supply chain resulting in a finished good or service as an exercise in 
theorem proving in linear logic. Each company in the supply chain expresses its capabilities 

25 as predicates and functions (to include constraint) in linear logic. The finished good or service 
is a logical expression of a set of predicates. A path or method of assembly of the final good 
can then be represented as a proof of the logical expression representing the finished good or 
service. If the good can not be built then no such proof will be found. The precise manner 
in which the product is assembled is read off the proof net describing all the necessary steps. 

30 The beauty of this approach is its simplicity and the availability of theorem provers for linear 
logic. Though theorem proving for linear logic can be NP complete proofs of supply chains 
may be much simpler. To summarize: the mechanism of self-organization in the supply chain 
is proving the final good or service. 

What are some of the axioms describing operations in the supply chain? There is shipping, 

35 manufacturing from components, etc. The state of the company includes inventory levels, lead 
times, etc. 

6.4.11 Learning in the Infrastructure 

There are many ways in which the data naturally recorded by the proposed infrastructure can 
be usefully used. 
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Prediction of customer demand: All the markets (both B2B and B2C) are rich sources of 
data for mining. Minimally, within any market we can use the data for both forecasting and 
inference of true customer desires (i.e. what they really want and now just what they are 
buying). The quantitative and nature of this market data should make mining it simpler than 
5 traditional data sources. 

Inference of consumer preferences: The present invention also improves upon the efficiency 
of the end consumers search by aiding the search with learning tools. Based on the consumers 
interaction with the top-level market software can extrapolate to guess preferences and speed 
the search. 

10 While the above invention has been described with reference to certain preferred embod- 

iments, the scope of the present invention is not limited to these embodiments. One skilled 
in the art may find variations of these preferred embodiments which, nevertheless, fall within 
the spirit of the present invention, whose scope is defined by the claims set forth below. 
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9 Abstract 



A system is described which allows buyers to define their preferences and sellers to define 
their capabilities, then determines which trading points maximize the utility of the buyer. The 
system suggests trades by exploiting the flexibilities and tradeoffs encoded by both parties, thus 
5 providing win-win trades. A second level of optimization ranks the trades with all suppliers, 
allowing the buyer to rapidly determine the best alternatives. The system allows for rich 
negotiation spaces and supports continuous, discrete, and range or interval decision factors. 
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A Distance Function for Range Variables 



The distance between buyer range variable r» = fe,rj) and seller range variable rj — (r^rj) 
is related to the overlap defined as 



i> r j) = / dx J\f(x;ri)J\f{x;rj) 



overlap (r j 



where M{x;ri) is a Gaussian in x centered at fa = (r^ + n)/2 with standard deviation oi = 
a(¥i — rj) and a is a tunable parameter. The range of integration is bounded by the seller's 
capabilities. Evaluating the integral, we find 



over lap (rj,rj) = 



2^27T(o?+<7?) 



'2(o? + aJ)J [ e 



Jj + (7jM — CTjUj + <?i 



By adjusting a we can control the hardness of the boundary at the ends of the range (smaller 
a corresponds to sharper dropoff outside the range). The maximum overlap will occur when 
fj,j = ^ and Gj = 0. Substituting those values in the above equation, we get 



maxOverlap = 



The distance d(ri,rj) is related to the overlap so that it varies from 0 (for identical range 
variables) to oo (for maximally distant range variables). Consequently, we define 



R(n,rj) = - In 



[" overlap (r^, rj)~\ 
I maxOverlap J 



B From Rankings to Numerical Values 

The numerical values in the tables Z{>ci,>cj) can be formed from the buyer's rankings in a 
20 number of ways. Since we are minimizing distance, we want preferred options to contribute 
little Z and progressively less favorable options to have progressively larger Z values. To 
approach this, we first consider the problem of defining a Z' value normalized to lie between 
0 and 1 and then take Z = - ln[l - Z'}. 

One possibility for Z' is a linear scaling. If di is the integer indicating the number of values 
25 in T>i, r\ a \i<j) indicates the buyer-specified rank (smaller rank means more preferable) for 
the option ^ = a € conditional on xj, then we might set Z'(^i = a, >tj) = (r\ a \vj) — 
l)/(di — 1). Alternatively, we might use a non-linear scaling and set Z(xi = a, vj) = (1 — 
^ (a) (^)-l)/(l _ ^-1) for some p € [0,1] (exponential) or set Z( Vi = a,Vj) = ((r[ a) {vj) - 



62 



type 




linear 
exponential 
algebraic 


ln L-r<°>(j 



Table 4: Possible numerical ranking functions. 



l)/(Ri — 1))" for some n > 0 (algebraic). The corresponding scaling functions for Z are given 
in Table 4. 

For those allowed values of >ti which are left unranked by the buyer, we can assign them 
a common numerical value greater than the lowest ranked value. For example if rj(xj) is the 
5 largest ranked value we might assign all remaining >Cj values to have Z{j<i, Kj) — ri(xj)/ (d* — 1) 
in the linear case or Z(x^) = (1— ^ r »C* f j))/(l-j9 d »- 1 ) or Zfaxj) = (l-l/(r i (x J )+l))/(l» 
1/di) in the non-linear cases. 



C Determination of Scaling Factors 

In this section we supply the details necessary to determine the scaling factors Q r and Qd- In 
10 the equations defining the average distances we have need of the integrals f v du exp[— (u — 
u)*C _1 (u - u)] and f v du (u - u)*C _1 ( u _ u ) exp[-(u - u)*C _1 (u - u)]. In general, these 
integrations are not analytically tractable; however we will assume that the region of integra- 
tion defined by x and x extends at least 2 standard deviations out in all directions so that the 
volume of integration is effectively infinite. In this case we can evaluate the integrals in closed 

15 form: 

f v du exp[-(u - u)'C-Hu - u)] „ 

and 

J du (u - u)*C -1 (u - u) exp[-(u - u)*C _1 (w ~ ")] ~ y 

To this point we have not been specific about what is meant by the average over the range 
20 variables. In fact, we will assume that, however we define this average, it is independent of 
what the buyer requested and thus can be removed from inside the sum over the discrete 
variables. In this case we will concentrate on setting the discrete and continuous components 
equal. We have (d c )/n c = {da) /rid where 

(d d ) = Q d Ev d d exp[-Q d rf d ] J v du exp[-rf c ] (dc> = 1 E v exp[-Q d d rf ] f v du d e exp[-rf c ] 

n>d ~ rid Ev exp[-Q d d d ] J v du exp[-d c ] n e ~ n c £ v exp[-Qd<f d ] fv du ex P[ _c? c] 
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and cancelling denominators Q d is determined by solving 

Q*Y d d exp[-Q d d d ] [ duexp[-d c ] = — y)exp[-Q d d d ] f dud c exp[-d c ]. 

nd ^ Jv n c y ./v r 

This can be simplified further given our assumptions about the Gaussian integral: 

Qd d d exp[-Q d d d ] _ 1 y> exp[—Qddg] 
n d y/detC- 1 ^) ~ 2 ^ y/detC-Hv)' 

5 In most cases there will also be at least one finite root, which is what we need. If there 

are a multiplicity of roots, we select that root having the smallest value to ensure that the 
distribution defined by the utility function is as flat as possible. If there are no real roots, 
then we select the Q d which minimizes the error. 

D Phase 1 Optimization 

10 In this section we detail the solution to the phase 1 optimization problem defined in Eq. (5). 

D.l General Formulation 

We begin by writing down the objective function: 

di(x,x) = (xW -pMYc- 1 ^ -**(*)) +R(t;x) 
where (see Eq. (15)) 

x« = (I-M«(*M {&W}))^ 

15 ' ' (15) 

as long as x^\b\ a l) < x^ < x^\b^l + 1). This must be optimized over x subject to 
the constraints Gi(x)x< s ' = gi(x) and G 2 (x)x( s ) = g 2 (x) for a given x. Since x is 
fixed throughout the optimization R{r;x) is just a constant which can be evaluated and 
then ignored as far as determination of x is concerned. We simplify our notation so that 

20 d x (u) = (x^ - /i) *C- X (xM R{r) . Similarly, 

xW = (I - MW({6SS»)- 1 c({6«}, {&$}) + (I - MW({6g}))- 1 MW({6g})x( 6 ) 
- ■({&$}. {$}) + S({&11 {^2})x^ 

where we have defined the n c -vector s = (1 - MW) _1 c and the n c x n' c matrix S = (I - 
25 M( s )) -1 M( & '. 28 n' c < n c is the number of continuous dimensions upon which the supplier 
response depends. Putting these pieces together we have 

di(x< 6 >) = x< 6 >'q({&$» {&$})x (b > + ^ x(&) + const (16) 

28 In the next section we consider S further, taking advantage of any redundant variables to simplify the QP 
task. 
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where Q is an n' c x n' c matrix given by 

q({*S5,{$!}) = s'^a^^^jc-^ciftW, {&{§» 

q is an n'-vector given by 

q(0>g, = 2S«({6S5, {6SS})C- 1 (8{6J2, {6iS» - 

and the constant term is 

const = {s{{b\%{bW}) _ -rfc-^sU^Ab^}) -n) + R(r) 

It is easily shown that the matrix Q is positive semidefinite. di(u^) must be optimized 
subject to the constraints 

GiS({&$, {&g>})x<» = gi - G l8 ({6g, {6iS». ( 17 ) 

G 2 S({^, {^2))x (6) < g 2 - G 2 s({6g, {&$}), (18) 

* W (0 < 8({&W {&$}) + S({6« < V i, (19) 

xW(0 <x< 6 > <#(i) Vi (20) 

where x^)' = [^(fcg) • • ■ x n (b%)] and x«(t)' = [a* (6$ + 1) • • ■ x n (b% + 1)] are 
the vectors of breakpoints for the f$ {k G [1, n c ]) functions. Similarly, for the buyer variables, 
xW(i)* = [xx^J) •-• and xW(i)* = [^(^ + 1) ••• * n (6g> + 1)] are the 

vectors of breakpoints for the {k G [l,n c ]) functions. 

di is optimized in two steps. First we assume a value for the integer variables {b\ s J.} and 
{bfl} fixing the constraints in Eqs. (17) and (18). We then minimize the quadratic di sub- 
ject to these constraints using a standard quadratic programming (QP) algorithm 29 to find 
xj p \({&g},{6g}). This result is then used to find x* (*:; {&$}, {&$}) = s({&$M&S» + 
S {©)^t({ 6 SS}> B y searchin § exhaustively over {&$} and {&$} ™ deter- 

mine 

x«(x) = argmin di(x*(x; {&$},{$}), x) 

Before turning to the manner in which this optimization is accomplished in section D.3, we 
discuss a few improvements to the input into the QP algorithm. 

D.2 The Quadratic Program 

In this section we discuss improvements in the specification of the quadratic programming 
task to make it better conditioned and more efficient. 

It is possible that the matrix Q used as input to the QP is singular or near singular, 

29 Numerous third parties supply quadratic programming packages, some based on interior-point methods 
which are appropriate, in the present context where matrix Q may be singular. 
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depending on the rank of Si. Consequently, it is important to use an interior point QP 
algorithm for which this is not a problem. 



D.2.1 Variable Elimination 

There may be rows or columns of S which are entirely zero. It is important to take advantage 
5 of this structure to eliminate unnecessary computation. 

If row i is zero, this indicates that supplier variable uf" 1 is independent of any buyer variable 
and is simply equal to s^. If there are any constraints involving only supplier variable i, then 
these can be checked before calling the QP solver to eliminate wasteful computation. 30 If the 
constraint is satisfied, we can eliminate it and then call the QP function. If column j of S 
10 is all zero, this means that there is no dependence of on and so there is no need to 
optimize over that buyer variable. We will, however, still have to check that the constraints on 
u( s ) are satisfied. This observation allows us to simplify the quadratic programming problem 
by eliminating this unnecessary variable. 



We formulate the problem to take advantage of both these observations. Assume there are 



15 m <n c non-zero rows of S and m! < n' c non-zero columns. Relabel the rows and columns so 



where Si is of size m x to' and, with a slight abuse of notation, the 0's indicate matrices of 
the appropriate size. In like fashion we write 



where x^ is an m vector, x^ is an n c — m vector, is an m' vector, and x^ is an n' c — m! 
vector, x^ indicates those components of x( s ) which vary with buyer variables while 
indicates those variables which are constant, x^ are those buyer variables which affect at 
least one supplier variable and x^ are those buyer variables which do not affect any supplier 
25 variables. If we also partition s as s* = s* s| then 



that 




20 




x[ s) = si + Six? 



:^ and x^ = s 2 . 



Finally, we will also have need of the partitioned covariance matrix 




where (C^.^)* = C 2 ^. w since C u ! is symmetric. 



'Additionally, the QP doesn't like to have a row of zeros in the constraint matrix. 
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Including the dependencies on /x, the quadratic form we must optimize is 

WsC^Wx W 2 C^W 2 ] [ s 2 -/x 2 J 
where we have made the obvious definitions 

-[T' 4.] -<•-£]• 

5 Expanding the quadratic, the constant term is 

( Sl -^)* W L C^ Wi ( Sl -Mi ) + (s 2 -M 2 )* W 2 ^^ 

The linear term is 

q = (( Sl - ^'WxC^WiS* + 2(s 2 - Ai 2 )*W 2 Cj ^WiS^x^ 

and the quadratic term is 
10 Q = (xf^SiWjCrjWiSix^. 

The constraints involving the optimization variables are 

Gi ; iSixS 6) = gl - Gi ;1 si - G 1;2 s 2 , 
G 2;1 SixS &) < g 2 - G 2 ,isi - G 2 ; 2 S 2 , 
x^W^sa+Sxxf^x^O") Vi, 
15 x^^^x^ <xj 6) (0 Vi 

where we have simplified notation by dropping the dependence on breakpoints and defined 

Gi=[G i5 i G i;2 ] 

i.e., Gj ; i is the columns of Gi involving the active buyer variables while Gj ;2 is formed from 
20 the inactive buyer variables. Similarly, the constraints which can be checked prior to even 
solving the QP are 

X 2 S) (i) <S 2 <X 2 S) (i) Vt. 

D.2.2 QP Normalization 

25 There are two types of normalizations we perform for increased numerical stability. If the scales 
of the continuous dimensions vary greatly (i.e., one dimension has an ideal value of 10~ 2 and 
another dimension has ideal value 10 6 ), then the matrix C _1 becomes badly conditioned, which 
makes S badly conditioned and thus makes the QP unstable. Similarly, if the breakpoints have 



[(Bx-^ + CxfVsl (S3 -/la)'] 
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radically, different lengths, the QP problem can also become unstable. We consider how to 
remedy each of these problems. 

Dimension normalization 

The quadratic form to be minimized is (x— /i)*C _1 (x— /*). We modify this problem by defining 
5 the dimensionless variables fa through xi = fain in terms of which 

(x - tfc- 1 (x - /*) = (0 - r?c-\ct> - 1) 

where C~ l = DC _1 D and D = diag(/i). When making this variable transformation we must 
assure that the constraints on x are adjusted so that GiX = gi and G2X < g2 become the 
constraints Q\<$> = gi and Qi<$> < g2 on <j) where = Gi j2 D. If there are global bounds of 
10 the form x < x < x, then these are translated to 4> < <j> < <f> where <f> = D _1 u and = D _1 u. 

In addition to these changes, we also need to modify the capabilities appropriately. Recall 
Eq. (2) which reads (after dropping the dependence on the discrete variables) 

k 

In terms of the dimensionless variables we have 

^ r 4MW) + 4MW) . {22) 

15 * y ^ K 

Thus, when we use the dimensionless variables, the function needs to be modified so that 
the y dimension is divided by ^ and the x dimension is multiplied by 

Once the optimization has determined an optimal <f>, we return to the original variables 
through x = D0. 

20 Breakpoint range normalization 

For larger problems it is sometimes important to normalize the input variables for the QP 
solver. This is carried out in order to improve the QP search performance in situations 
where there is large variation in the range of the segments (bk) < 4$ < (bk 4- 1) and 
the magnitude of their values. This can be accomplished by introducing the normalization 
25 variables 4>f = (^ - (b k )) / '(<^ 6) (b k + 1) - 4> k b) (b k )). Here ffl will have minimum value 
of zero and maximum of one. We can now express 4>^> as 

^) = A^ (6) +p (23) 

where A is the diagonal matrix A = diag(0? - <t> { ^\$ 2+1) - $ 2 \ ■■■ , cj>g m+l) - <f& m) ) 
and = [<$ 1 \$ 2 \--- ,4> { X m) }. 
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Expanding Eq.(16) using the definition of 4> in Eq.(23), we obtain 
rfi(w; {h, k }) =^ (6)f ADQ({6 i ,,})DA^ ) + q<({&;, fc })DA^ + p*DQ({6 i)fc })DA^ (b) + 
ip t DQ({6 i)fc })Dp + q t ({^, fc })Dp. 
Hence the quadratic programming problem to be solved in terms of <f> is 

QP(±W; {k, k }) = l^Qdb^})^ + q t ({bi, k })n (b ), (24) 

where Q = ADQ({6 ijfc })DA and q* - (q*({6i,fe}) + P i DQ({&i ! fc}))DA. Similarly, the equality 
and inequality constraints of Eqs. (17) and (18) expressed in terms of <p are 

6i({M)* W =gi and G 2 ({b i>k })j> {b) <g 2 , (25) 

where 

G!({^}) = G 1 S({6 i , fc })DA, 

G 2 ({b i>k }) = G 2 S({6 i)fc })DA, 
g! = gi - G 1 s({6 i>fc }) - G 1 S({6 iifc })D P , 
g 2 = g2 - G 2 s({& ilJb }) - GiS({& i)fc })Dp. 

D.3 Optimization Over The Breakpoint Variables 

We have detailed how to determine xg!t(x, {&$}) S iven a specification of all {&$} and 

{bfl} y but how do we determine the best feasible choice of {&$} and {&£>}? Before we answer 
that question, we have to examine the infeasibilities that might arise due to improper choices 
of{6$}and{&g}. 

D.3.1 Elimination of {b itk } Infeasibilities 

Most of the settings of {&$} or {&$} will not yield a feasible solution because Eq. (20) or 
Eq. (19) will be violated. The reason is that, for any k, it is highly likely to get two segments 
xf{bi, k ) < 4 6) < xf{bi, k + 1) and 4%, fe ) < < x^(b jJt + 1) that do not overlap (i.e., 
the optimizer will not be able to find a value that satisfies both bounds, and hence no 
feasible solution will be found). The same is also true for supplier response segment . 
To eliminate those infeasibilities, we construct a super set of breakpoints {va b k } such that 

4" = IM^S = - «). ■ ■ ■ .«f(»!S = 

i-l 

Hence each piecewise linear function f i)k will be described by a set of breakpoints w k = 
{x k {l), ■ ■ ■ ,x k {n k )} and function values / i>fc (a; fc (l)), ■ • • , fi, k { x ki K k))- These function values 
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at the new breakpoints are linearly interpolated from the old breakpoint-function pairs. After 
this reduction we have — 1 different € [l,«fc — 1] values, which eliminates the need for 
constraint Eq. (20). 

The continuous supplier variables and constraint Eq. (19) are more problematic. Sup- 
5 pressing any dependence on x, we know 

- «««&••• .e,M$>>+E^$Mr> +E«SM>*? ) - 

k<i k 

If {bfl) and {&$,••• ,6^_i} 

are known, then this expression allows us to bound the values 
x\ 8 ^ may take on, and thus bound the interval in which the feasible 6^ lie. If we define the 
functions 

{sm<Lx(x,y) ifs>0 I smin(a;, y) if s > 0 

and smin(s, y) — < 
srcim.{x,y) if s < 0 I s max(a;, y) if s < 0 

then the minimum and maximum values that a;^ may take on are: 

min(4«>) = . . . , jWj, {(,(«}) + £ S mm( m (J(6,W),4'' (6^) . 4"' + 0) 

+E^(^2).*? , ('drxf ) (»i2+i)). < 26 ) 

A; 

15 and 

max(^) = c ,({6g, • • - , frSjlj, + E smax K-;2( & S)' 4 5) (^2):4 S) (^ + 1)) 
+ E—^cS!).*?'^).-?'^ + «) < 27 > 

respectively. The above equations imply that we can fully translate a complete specification of 
buyer variable breakpoint intervals to a corresponding set of feasible supplier variable break- 
point intervals. We will search over the breakpoints using a branch-and-bound algorithm, 
and will need to be able to infer the restriction on the values of feasible supplier breakpoint 
intervals given a partial specification of buyer breakpoint interval variables. This inference, 
which narrows the domains of the supplier breakpoint variables, is called constraint propaga- 
tion. Before describing how the branch-and-bound algorithm search over breakpoint interval 
variables is accomplished, we show how the constraint propagation may be obtained. 

Narrowing of Supplier Breakpoint Domains Through Propagation 

Let jSj^ be a vector of integers specifying a set of breakpoint intervals for the kth buyer 
variable. \fik\ = 1 if the interval is uniquely determined, and \/3k\ = — 1 if no restrictions 
30 on the interval have been specified. We index the Ith. component of 0^ as I £ [1, - 1]. 
A set of breakpoints for all buyer dimensions is written as {^j^}. In analogous fashion, we 
define fij^J, I € [1, — 1], and {/3fc} for the supplier breakpoint interval variables. 



20 



25 
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Prom a specification {/3 k b) } we infer the feasible domains {p[ s) } using max(^ s) ) and 
min(a:^) where: 

max^) = £ ^{^({^ + !))} + 

g 6 J^J c ^^ + 1))} (28) 

and 

min(*W) =J2 min ■ {^({6«}) + smin^^^),^^), + !))} + 
fc b fc £ p* J 

^.ra J** ({6 * )}) + s ™< m ^\4 s \i\4 s \i s) + 1))} (29) 

k<i b k €01 } 

We begin with i = 1 (where no supplier variables contribute to x[ s) ) and determine the interval 
h = [min(4 fl) ),max(a;[ s) )]. We then form the vector p[ s) with components {p{'J\[x^ ) (fi li i),x[ s) (j3 li i+ 
1)] n Ji ^ 0}. Equations (28) and (29) can then be iterated to find the interval I 2 = 
[min(x^), max(4 s) )] using the previously calculated min(a;^ ) ) and max(4 s) ). With I 2 known, 
the components of fi^ are determined as above. Iterating this procedure, we can find the 
potentially narrowed domains of all {/3^}. 

We now return to the optimization problem: we need to search over all H k ^ x Hk ^ 
combinations of breakpoints to determine the combination that minimizes d\ (x opt (>*, {b^}, {b[ s) }), x). 
The simplest manner in which this accomplished is via exhaustive enumeration over all {b^} 
and {b ( k s) }. However, once all b[ b) have been specified, we need only search exhaustively over 
the corresponding {&(s)} inferred as above. There are, however, many ways to speed up the 
search over the buyer breakpoint intervals. We describe two methods. 

D.3.2 Local Search 

This method is extremely fast but not guaranteed to return the global optimum across all 
buyer breakpoint intervals. Here, initial random and feasible 31 {&jj. a) }, and {b^} sets are first 
created. After the continuous optimization is carried out using the initial {b^}, and {b^}, a 
new {b^}, and {bf^} set is generated from the old one by examining whether each generated 
4 6) ({4 S) }' solution is on the boundaries of the interval x[ b \b k ) < x (b) < x { k \b k + 1) 

or x ( k s) (b k ) < 4 S) < x k s) (h + 1). If 4 &) is on the lower bound of the interval, ~bf> = b {b) - 1; 
similarly, if it is on the upper bound, b k b) = b {b) + 1. If x { k b) is not on either end, then 
= 0f course this done while maintaining b k b) in the interval 1 < b[ b) < nf* - 1. 

Supplier breakpoint variables are updated in a similar manner. The procedure is repeated 
until no decrease in distance is attained. 

"Feasible in the sense that they do not violate the constraints in equations (17), (18), (19), and (20). 
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D.3.3 Branch-and-Bound Search 



Branch-and-bound is an efficient procedure used to exhaustively search a space and eliminate 
configurations that can be proven to be suboptimal. The theory behind branch-and-bound 
is well understood [5]. The idea is to structure the space as a tree of partially-completed 
5 possibilities and to explore the tree so that good complete solutions are constructed early on. 
Branches of the tree that can be shown to yield solutions which cannot be better than the 
currently best solution need not be explored. Heuristics indicate the order in which to explore 
the tree and provide a lower bound on the fitness of partially completed solutions for pruning 
of suboptimal branches.. 

10 In the present context, we carry out branch-and-bound search over the space of possible 

buyer breakpoint combinations. The tree is structured beginning from a node at which none 
of the have been assigned. From this node, the tree branches into the nf^ — 1 possible 
assignations for bf^ . From any of these nodes, the tree branches into the Kp possible values 
for b^P . The tree branches at each node to assign the remaining variables. At the leaves of 

15 the tree, all buyer breakpoint variables have been assigned. 

From any node, the search proceeds by choosing the most promising node to consider, 
as indicated by a variable ordering heuristic and considering the possible values for that 
variable in the order indicated by the value ordering heuristic. A complete solution and 
its accompanying distance is obtained when the search ends at a leaf of the tree. To avoid 

20 unnecessary exploration, a lower bound heuristic indicates when a partially completed solution 
(node) need not be expanded because all completions will have higher distance than the best 
located so far. We describe these three heuristics for search over buyer breakpoints. 
Lower Bound 

Given a set of buyer breakpoint interval vectors {ft^}, we can propagate to {jS^} and the 
25 intervals [min(a;^),max(x^)]. A lower bound on the value di can be obtained by ignoring 
the supplier capabilities and minimizing 

d 1 (x) = (x( s )- M )*C- 1 (x«-/x) (30) 

subject to 

Gix( s ) = gl (31) 

30 G 2 x^ < g 2 (32) 

mm{x\ s) ) < x[ s) < max(4 s) ), (33) 

The lower bound thus obtained is denoted by lb({fi^}). 
Variable and Value Ordering 
35 The lower bound \b({^}) can be used to define an ordering over variables and the values 
any variable may take on. We define the functions L k j which operate on the set {fi^} by 
locking the kth buyer breakpoint variable to At any node we can determine the lower 
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bounds lb(L kjl {pW}). We can select to explore the variables in order of increasing lb^ where 

Eli 1 1 

i.e., the average over all possible completions for b^P Within any given variable the (say 
k*) values can be explored in order of increasing lb(Lk* ; i({/3^})). If we are willing to do 
5 more computation, we can explore the branches in an interleaved order according to all the 
lb(L fcji ({^)})). 

Pruning and Search 

During the branch-and-bound search, branches can be pruned if the node has lb(l^ ; ;({jS( 6 )})) 

greater than the distance of the lowest examined leaf node. 
10 Search begins at the root node where each 0^ = [1, • • • — 1]. We applying the 

ordering heuristics above to move to a node where one breakpoint variable has been fixed. 

From that node, we move down the tree, fixing a value at a time until we reach a leaf. We 

then progressively backtrack up the tree, eliminating branches whose lower bounds axe greater 

than the best leaf examined. 
15 Once the optimal set of buyer breakpoints has been determined, we search exhaustively 

over the supplier breakpoints feasible with this choice according to Equations (26) and (27). 

Improvements 

We have not yet exploited the fact that, from any node having {fi^}, the constraints in 
Equation (19) may narrow the domains of the buyer breakpoint variables deeper in the tree. 
20 The computational cost of narrowing the breakpoint domains may be warranted if the pruning 
thus obtained is sufficient. 

E Phase 2 Optimization 

In this section we detail an efficient algorithm to optimize Eq. (5) : 

arg min d\ (x(x), x) + d 2 (x) subject to C (6) (x) A (jW (x). 
25 Our approach is a marriage of constraint programming and local search. 
E.l Local Search over ID-neighborhoods 

We begin with a definition. The (^-neighbors of a discrete configuration x are all those con- 
figurations, x*, which differ from x in d attributes. Attribute i differs in two configurations 
if Xf 7^ xj. Thus the distance varies from 0 if the configurations are identical, up to na if all 
30 attributes are different. The ^-neighborhood of x is the union of all d = 1 neighbors, d — 2 
32 We might also take the minimum over all possible completions lb fe ({/S w }) = min J6[1 lb(L k ,i({P (b) })). 
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Table 5: Configurations examined in the D = 3 neighborhood search around (a) the initial 
configuration 11111 and (b) the initial configuration 0+0+0 (duplicate search configurations 
are indicated in bold) . 



neighbors, up to d = D neighbors. The ^-neighborhood around x is indicated by N D {>c)- We 
note that, for any configuration, its n d neighborhood is the entire discrete space. 

Using the ^-neighborhood, we apply a standard hillclimbing algorithm which respects all 
the discrete constraints. D is specified before the algorithm begins and serves as an efficient 
5 manner to trade off running time with the quality of the final solution. The hillclimbing 
algorithm takes as input a starting configuration x t and finds the configuration x t+i e N D (x t ) 
which minimizes d 1 (x{x),x) + d 2 {x) subject to CW(x) AC< s >(x). The algorithm then iterates 
again starting from x t+l . The hillclimbing algorithm terminates when it finds a configuration 
x T which has a distance less than or equal to all of its D-neighbors. If D = n d , this is 
10 the global minimum. Given a function which exhaustively searches a ^-neighborhood for 
the lowest distance feasible configuration, a D-neighborhood hillclimbing algorithm is easily 
implemented. 

There is one noteworthy point about the implementation of an efficient ^-neighborhood 
hillclimber. Consider a D = 3 search over n d = 5 discrete variables. Without loss of gen- 

15 erality, we label the initial configuration as 11111. In the D = 3 neighborhood search, the 
configurations examined are given in Table 5(a) as indicated by the *s. Let us assume that the 
best neighbor in the D neighborhood of 11111 differs in components 2 and 4. We indicate this 
best neighbor as 1+1+1, which serves as the starting configuration for a new D = 3 search. 
The configurations examined in the second neighborhood search are given in Table 5(b). Prom 

20 Table 5(b) we see that the search branches 2, 7, and 10 are redundant, having already been 

performed in the D neighborhood search around 111 11. 33 Thus these branches may be pruned 

in the second D neighborhood search. A search branch covers a previous search and may be 

pruned if the vector of positions of t he wildcard characters covers the vector of positions at 

33 The duplication of search effort is most acute when D is large and the fittest neighbor is found very close 
by (small d). 
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which the previous best solution was found. The binary relation covers is defined as 



x covers x' x'- £ x V x| € x'. 

The covers relation is transitive but not reflexive. Thus, to increase the efficiency, we can keep 
a list of the previous branches searched and ignored search branches that have already been 
5 covered. 

E.l.l Improvements 

Throughout the present discussion we have assumed that D remains constant through the 
hillclirnbing process. There are good reasons why this should not be the case. With any 
reasonable landscape correlation structure, we know that lower distance neighbors are more 
10 likely to be found early in the hillclirnbing process. Thus a reasonable heuristic is to start 
with a larger D early on and gradually lower it with time. 

E.2 Exhaustive Search of a ^neighborhood 

To complete the discussion of the phase 2 optimization, we describe a computationally efficient 
way in which to find the lowest distance feasible configuration within a ^-neighborhood. We 
15 use a branch-and-bound method with heuristics tuned for the present problem. The factors 
tuned include the variable order in which we search, the order of the values any given variable 
may take, and a heuristic which says when certain branches of the search space may be safely 
ignored. 

We use heuristics which early on in the search generate low-distance configurations, and 
20 later prune parts of the search that we can prove must have higher distance than the best 
configuration generated so far. 

E.2.1 Heuristics And Lower Bounds 

We focus on obtaining good heuristics for the discrete contribution to the total distance, then 
later incorporate the continuous and range contributions since these are simpler. 
25 At any point in the search, some variables have been assigned and some have not. We 

divide up the sum over the Z^x^, Xj) into four contributions 

Q d Z w (*) = Q d ^ S u>iU)jZiA*i>*i) + Qd zZ WiWjZijfaKj) 

i locked j locked i locked j floating 

i floating j locked i floating j floating 

30 = Zl > 1 + Z ^ + Z U + Z U 

Since we are minimizing and all Zjj(xj, Xj), wi, and Qd are positive, we can prune any search 
node which has Z\ t \ > than the distance of the lowest distance configuration found thus far, 
since additional terms will only add to the sum. 
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We obtain a lower bound on the distance by lower-bounding each of the above four terms. 
Just taking the lowest value in the tables can yield differing values for any given variable, e.g., 
the values for a single variable xf. 

x* = arg^min Zijfa, xj), x* = arg min Z^{x h Xi), and x\ = arg min Zi^Xj) 

may all be different. To correct this we write 

QdZ w = z u + z iti + J2 z S^i) 

i floating 

where Z\^ and Zf >t are calculated as before and 
j locked 

The best choice for the floating variable x^ is 

= argjnin Z$(xi). (34) 

In this way all the possible choices for xi are at least consistent, and the true lower limit is 
still guaranteed to be higher than this (since the Z fjf contribution can only be higher). 

This approach also naturally suggests a variable ordering and a value ordering within any 
variable. We can try the variables in increasing order of their corresponding Z$(x%) so that 
variable i, which has the smallest is tried first. For any given variable x 4 we can also 

order the possible value assignations for «j according to the array of values Z$(xi) so that 
the value x { = xf is tried first and the value which maximizes Z^(xi) is tried last. 

With little extra work we can impose a consistency requirement raising the lower bound 
for Zfj. We decompose the sum as follows 

Zt,f = Q d WiWjZijixi, xj) =Q d w iWj {Zij^i, xj) + Zj^xj, Xi)) . 

i,j floating:.?^ jj floating: j <i 

Instead of taking the minimum component-wise as previously suggested, i. e instead of 

i,j floating:,^ " 3 

we take the minimum so that 



Zf,i = Qd J2 ^[ w i w 3( z iA*i,*j) + ZjAxj, **))]• (35) 

i,j floating:j'<i 

which enforces pairwise consistency between variables i and j. 
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The lower bound, lb, we have obtained is thus 

ib = z 1)1 + ^ f + £z»(^) 

where x? is determined from Eq. (34) and Z f * f is determined from Eq. (35). 
Integration Including Continuous and Range Contributions 

In this section we summarize the complete heuristic specification, including the effects of the 
continuous and range variables. As above, we use the definitions 

i,j locked 
j locked 

Z h = Qd ™m[w iWj (Zij («t, yc-) + Zj^Ki, x,))] 

i,j floating:j'<j 

so that a lower bound to the complete discrete contribution to the distance is Z^\ + Z\ { + 
J2i z \,f{Kt) where = arg min zj; % J(xi). 

We now extend this result to include the continuous and range variables. The continuous 
and range variables depend on the discrete variables through the mask. Since we expect there 
will be relatively few masks we exhaustively enumerate all the discrete possibilities. Conse- 
quently, we define d x (boating Nocked) as the distance contribution made by the continuous 
and range variables (the d x contribution to Eq. (5)). We define the continuous + range 
contribution for any given floating variable i as 

= flo mg^dl (^floating l^locked) 

where the min is taken over all floating variables other than i and where floating variable i is 
locked to value x;. 

Since each continuous/range/discrete distance contributes additively to the total distance, 
we can rank the possible values for any given floating variable i according to Z^fa) + d ( (\vi). 
The best (according to this heuristic) setting for floating variable i is 

x\ = arg min {z^fo) + <4° («»)}• 

The heuristic then suggests a variable ordering in order of increasing zf}{>cf) + d^(vf) so 
that we first select the floating variable with the lowest sum and last select the variable with 
the highest sum. 

Finally, we turn to the lower bound on the distance. Noting that the continuous and range 
distance contribution does not have the form of a sum over the dimensions (as the discrete 
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contribution does), we write the lower bound, lb, as 34 

lb = Z u + + ^df{x*)) + Zlz 

where i* = arg min,- + '(«?)) (i.e., i* labels the floating dimension with the least 

z \} + ^\ contribution to the distance). 

E.2.2 Improvements 

We point out that an extensive literature exists on parallel branch and bound algorithms. 
Many of the techniques developed are applicable to the present problem. 



is the Kronecker delta function which is 1 if and only if x = x' . 
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Figure 1: An architecture for the invention. 
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Figure 2: Schematic of a buyer-specific capability with examples indicating potential input. 
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Figure 3: Schematic of a supplier-specific preference with examples indicating potential input. 
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